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1. Historical introduction*. The theory of bending of elastic rods was initiated 
by JamEes BERNOULLI in 1691. In his work, as well as in a paper of VARIGNON 
from 1702, the basic principles of statics as applicable to a rod are distinguished 
from the particular hypothesis of elasticity, but this separation has often been 
neglected by later writers. An exception is furnished by a remarkable memoir 
of EuLEer [1771], where the fully general statical equations for a rod bent in its own 
plane are derived. While EuLER developed also a theory of bending of skew rods, 
all the researches in the seventeenth and eighteenth centuries remained primitive 
from lack of an adequate description of strain in a rod. 

Such a description was initiated by St. VENANT [1843] [1845], who introduced 
the concept of twist and the principal torsion-flexure axes (cf. also Binet [1844)). 
Other attempts to formulate a theory of strain in a rod were made by KIRCHHOFF 
[1859, § 2] [1876, Vorl. 28, §2] and Crepscu [1862, §§ 48—49, 55]. All this 
early work employs more or less hidden approximations and is difficult to follow 
with confidence. The first straightforward analysis is that of Love [1893, 2, 
§ 233]. While a great advance beyond its predecessors, his work is closely bound 
to the concepts of small deformation and linear elasticity. 

That six equations are needed to express equilibrium of a bent and twisted rod 
was first remarked by St. VENANT [1843], but he did not succeed in obtaining 
them without simplifying hypotheses. The exact general equations were given 
in principle, but obscurely, by KiRCHHOFF [1859, § 3], explicitly by CLEBSCH 
[1862, § 50]. These and other early treatments are difficult to follow, sometimes 
imparting the impression that some approximation is made. In fact, as was noted 
by Basset [1895, § 2], the statical equations referred to the actual position of 
the rod are exact; no question of approximation arises unless we attempt to refer 
the equations to a configuration assumed by the rod prior to its being subjected 
to the load under which it is in equilibrium. 

The theory of curved elastic shells is of more recent date. In the pioneer paper 
of Love [1888], neither the statical equations nor the description of strain is 
disentangled from the special elastic hypotheses and simplifying approximations. 
The exact statical equations were obtained by Lams [1890]. The descriptions 
of strain used by classical writers are rather intuitive than proven and egress 
little, if at all, from the domain of infinitesimals. 

The flood of papers on rods and shells brought on by the twentieth century 
has added little to principle. When not erroneous, in respect to the foundations 
they are for the most part repetitions or amplifications, sometimes more con- 
cisely expressed, of older ideas. Even works devoted expressly to finite deflections 
turn upon retaining additional terms in series expansions, etc. 

Exceptions to the foregoing summary are the paper of Hay [1942] on rods 
and that of SYNGE & CHIEN [1941] (cf. also CHIEN [1944]) on shells. These authors 
obtain correct and general descriptions of finite strain, independent of the elastic 
hypotheses adopted in the later parts of their papers. While they employ ten- 
sorial methods to describe the strain, they follow the classical custom of selecting 
special co-ordinates. Not only are these co-cordinates not necessarily the best 


aes full history of the early theories of rods will be included in The rational me- 
chanics of flexible ory elastic bodies, 1638—1788, Editor’s introduction to L. Euleri 
Opera omnia II 11, forthcoming. 
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suited to the statement of each special problem, but also in the theory of finite 
deformation they give rise to an unsymmetrical if not formidably elaborate set 
of equations. In their treatment of stress, however, SYNGE & CHIEN by using 
fully general co-ordinates obtain a most elegant and compact formal statement. 

A new idea, supple in application to a variety of mechanical theories and 
formalisms, was proposed by DuHEM [1893, 1, Ch. II]: A body is to be regarded 
as a collection not only of points but also of directions associated with the points. 
These vectors, which we shall call the directors of the body, are susceptible of 
rotations and stretches independent of the deformation of material elements. 
Such a model of an oriented body should include, among other possibilities, repre- 
sentation of a molecular concept in which the molecules have internal structure. 
But also, as E. & F. Cosserat [1907] [1908] [1909] remarked, in one and two 
dimensions it serves admirably to represent the twisting of rods and shells in 
addition to their bending. The CossErats constructed descriptions of this kind. 
These descriptions are sect in a single rectangular Cartesian co-ordinate system 
and are formally elaborate from lack of any direct notation; moreover, they are 
interjected in a long plea for a kind of generalized elastic strain energy which 
the CossERATs called ‘‘Euclidean action”. Except for an exposition by SUDRIA 
[1935]|*, this profound work of the CossERatTs has attracted no attention. Their 
considerations, while seeming to us incomplete as well as cumbrous, suggest the 
general theory of strain we now formulate. 


2. The nature and plan of this work. For the notable successes of the three- 
dimensional theory of finite elastic deformation, as well as more general theories 
of recent date, an exact description of stress and strain was a necessary prelimi- 
nary. Similarly precise and general theories of finite deformation of shells and 
rods might have been expected to follow upon these now well understood and 
developed branches of mechanics. That they have not, reflects, in our opinion, 
the lack of a precise and general description of stress and strain in rods and shells. 
Such a description, divorced from any constitutive assumption intended to de- 
scribe the elastic or plastic response of the material, we here construct. 

First we give some mathematical preliminaries and a sketch of those properties 
of an oriented body which are independent of the number of dimensions. There- 
after follow the theories of strain of rods and shells: 

4. Strain of position 

a. Intrinsic theory 
B. Imbedding theory 

2. Strain of orientation. 

The intrinsic theory, based on the first fundamental form, coincides with the 
classical theory of strain; the imbedding theory rests on the curvatures and is 
known from differential geometry. A brief restatement of these results is included 


for completeness. The bulk of the paper concerns the strain of orientation. Here 
again the analysis is divided into two parts: 1°, differential description of the 


* In his §9 SupRIA notes an error in the CossERATS’ argument and gives a dif- 
ferent proof of invariance. We do not investigate the details, since our treatment 
is simpler and more general from the start. We note that in his Ch. II Suprra analyses 
time rates, a topic not included in the present work. 


AM 
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undeformed oriented body, and 2°, differential description of the deformation. 
Proofs of invariance and completeness are given. Despite our acknowledged debt 
to the CossERATs, we consider this analysis new in the main; in particular, its 
approach differs fundamentally from those of Hay and of SYNGE & CHIEN, whose 
results we pause to include as special cases for specially selected co-ordinates. 

As is well known, it is far easier to obtain an exact description of stress than 
of finite strain, and even for rods and shells the only improvements to be expected 
here are in method. For rods we give essentially HEuN’s [1913, §19] compaction 
of the CossERATS’ argument. For shells we obtain the general equations of 
SyncE & CHIEN by a briefer method made possible by use of the concept of 
double tensor field and total covariant derivative, explained in § 3 and used in 
our theory of the strain of shells. While we regard our material on stress as 
essentially classical, we include it not only for completeness but also for method. 
As far as stress is concerned, our line of argument, essentially that used by HEUN 
for rods, is uniformly valid for a surface of # dimensions imbedded in a Euclidean 
space of #-+q dimensions. 

The reader is assumed familiar with the classical notions of stress and strain 
in three dimensions as explained, e.g., by TRUESDELL [1952], but no prior famili- 
arity with theories of rods and shells is needed. 


Part I. Preliminaries 


3. Some mathematical tools. We consider two spaces, one with real co-ordinates 
X*, K=1,2,...,.N, or simply X, and co-ordinate transformations AX (RDS 
and the second with real co-ordinates x*, k=1, 2,..., n, or a, and transformations 
x'=a2x'(x). A set of functions T§:p2:? of the pair of points X and @ is a double 


tensor field if it transforms according to the rule 


TK PhP TS US nth BOTS WOK eek ex 
Q...Rq...t V..Wv...w ax axu ax’0 tee axe 
' 3.1) 
S) Cee ax? ax” ox” ; 
ox* Ox" ax’ 4 Ox’? 


A familiar example of a double tensor field is the vector connecting two points 
in a Euclidean space: If #* and PX are the position vectors of the two points 
in the co-ordinates selected at a and X and if g¥ and gm are the shifters* that 
effect finite parallel displacement, then the vector connecting P to p is given by 


gk ph PK or pt—gh PK, 
according as the co-ordinates at X or at # are employed. 


Zs Let co-ordinates of X and # in a single rectangular Cartesian system be Z and z 
The shifters are then given by 


gf aah OE az gm = pt O" oZk 
OZ Og 


Alternatively, they may be defined as those double tensors of the indicated variance 


such € at W hen both Ss Vy S) tems oO SOR inates are th ame I ect an ular Car tesian 
h if co d ers 
system we have oM = Oo! &u = Onr S 
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Let both the space of a and the space of X be metric. Covariant differentiation 
with respect to x* and X* may be defined in the usual way, if we adjoin the con- 
vention that X 7s held constant when we differentiate with respect to x", and vice 
versa. These partial covariant derivatives, which we denote by the usual symbols 


ce »”» +” 


and “x’’, are double tensor fields of the type indicated by the number and 
positions of their indices, and the formal rules of ordinary covariant differentiation 
remain valid. 

Now. suppose 2 2 N and consider a mapping =a (X). The gradients dx*/ax* 
are double tensor fields, which we denote by 


as ax* 
aK = Kr (3.2) 
The total covariant derivative of a double tensor field 7: is defined by 
De Tea Ki (3.3) 


For a given T(x, X), the total covariant derivative T--. x has a value independent 
of whether or not « be eliminated through the functional relation #=a(X) in 
some or all of its occurrences in T(#, X). This property is not shared by the 
partial covariant derivative T-: x. These remarks are most easily illustrated 
for the case when 7 reduces to a scalar function T(x, X) of two real variables 
x and X; the partial covariant derivatives then reduce to 0T/0x% and 0T/éX, 
while the total covariant derivative reduces to 67/0X + (0T/0x) (dx/dX). 


The formalism just introduced we shall use in two senses. In the second, 
n=%3 and N=1 or 2, and the space of X is a curve or surface imbedded in 
Euclidean three-dimensional space. In the first, »=N=3, and we regard X 
and @ the co-ordinates of the same point before and after deformation of the 
material; in the now usual terminology of the subject, X is a material co-ordinate 
and a a spatial co-ordinate. The relation «=a (X) is assumed uniquely invertible 
to yield X = X (a), and we define 


ee a (3.4) 
TORN OE vig yo he Pw Gore (3.5) 

We have also the chain rule: 
ADEA SOUS gs A Ecareen wt Beare a Gare (3.6) 


In this memoir all considerations are local, and the functions occurring are 
assumed to have as many continuous derivatives as needed to justify the formal 
operations. 

4. Formulae from the general theory of strain of position. We consider two 
metric spaces of the same dimension m, with points X and a and with positive 
definite squared elements of arc: 


dS? = gx yd X* dX™, AS == 0,4 a xe. (4.1) 
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Under a deformation #=a(X), X=X(«), we have 
dS? = Gy dx dx", ds? = CgyaX* dX™, (4.2) 
Gi ey eee Crm = 8am “ik XM 
The tensors ¢ and C are the deformation tensors. Since they are positive definite, 
Sal of 
they have unique inverses and square roots, which we denote by C,C, ete. ihe 
principal axes of e and € are called the principal axes of strain at x and X. 

A non-singular real matrix a has the unique decompositions a=s-0—0- ss 
where 0 is an orthogonal matrix and where s and s* are symmetric and positive 
definite matrices. In a Euclidean space, application of this polar decomposition 
theorem shows that 


atte = Rhy CH =c"* R%, (4.3) 


where C and ¢ are defined above and where R is a rotation tensor*. This is 
a formal expression of the theorem that any local deformation may be resolved 
uniquely into a pure strain followed or preceded by a rotation. The rotation 
carries the principal axes of strain at X into the principal axes of strain at a. 


Part II. Theory of strain 


A. General theory of strain of orientation 


5. Differential description of the unstrained body. To the point X assign a 
set of p vectors D,(X), a=1, 2,...,p, the directors of the body at X. By a 
deformation we shall mean a transformation carrying X into x and the directors 
D, at X into directors d, at x. In equations, 


a= a(X), d, = d, (Dp) = dy (X); (5.1) 


in geometrical terms, a deformation consists in a displacement of the points and 
independent rotations and stretches of the directors. In the special case when 
di, = x* x DX, the directors are material elements, and their presence adds nothing 
to the usual description of strain**. In the special case when d* = gk DX, the 


* An orthogonal tensor o in a metric space is characterized by the property 
Sim nO" = Opa: When &% m= 4m, the matrix || o*,|| is an orthogonal matrix, but 
in general co-ordinates it is not. In a Euclidean space we may displace fields from 
one point to another by means of the shifters g& and gk, as mentioned in §3. A 
volation tensor such as that occurring in (4.3) is a shifted proper orthogonal tensor: 
Re, = oe Cee det ss = 1. 

| ** An exception should be made for the theory of anisotropic solids, where the 
directors D, are material elements selected on the basis of a priori knowledge con- 
cerning the nature of the undeformed body. The symmetries of the material. are 
stated in terms of invariance with respect to certain transformation of the Dg. Since 
di = ah D*, from (4.2) we obtain 
—1 
Ban = Cxm DED, Cem =8qyDeE Dy, 6 *™ = GP dh ay, 


etc. A formalism of this kind for the strain of anisotropic bodies has been constructed 
by EricksEn & RIvLIN [1954, § 2]. 
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directors are invariable elements and again add nothing®*. In general, the directors 
are neither material nor invariable; 1.e. 


dhtxteDE and dk-+egh DE. 


In an oriented body, strain and rotation are defined from (4.2) and (4.3) in 
the classical way. What is new is the relation between the directors. For the 
full range of interpretation mentioned at the end of §1, it should be possible 
to use an arbitrarily large number of directors. This degree of generality, 
however, we do not attempt, being content to consider in a space of 2 dimensions 
a set of m linearly independent directors D,. This number suffices for theories 
of rods and shells, which we regard as curves and surfaces embedded in a Euclidean 
space of three dimensions, the maximum number of linearly independent directors 
thus being, in both cases, three. 


Let D® be the reciprocals to the D,, so that 
DED YS Oy 1 Deny een DS PAGED ar, (5.2) 


where ** G,, = D* Di! ge y. The D* are the reciprocal directors. 


The director triads may be used to define anholonomic components. For 
example, if we set 
p= DEX, (5.3) 


then X*, = DX X$, and from (4.2), follows 


Cem = Gay XF v4 etc. (5.4) 
Now set i 
W*up = DE p Diy = — Dk M:P> (5.5) 


where for uniformity with later developments we use the notation (3.6) of the 
total covariant derivative, although of course DX p= D& p. From (5.5) follows 


Diu = Woy Dy. (5:6) 
From (5.5) and (5.2) we find that 
2Wx mp = Gav; p Dk Du- (5.7) 
Therefore G, = const. is equivalent to 


Weup = — Wuxr- (5.8) 


* Selecting three orthogonal and invariable directors, we may use them as a 
fixed anholonomic frame and so obtain invariant forms for the results given in many 
of the older treatments of finite strain, where all quantities are referred to a common 
rectangular Cartesian co-ordinate system. 


xx Note that 
det Gap = det (Dg Dy cos %q5) = (D1 De... Dy)? det cos Dqg, 


where Dg is the length of Dg and 8» is the angle between Dg and Dy .When the 
directors form an orthogonal unit set we thus obtain det Ggp= 1. 
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In this case, also, if we transform the D, at all points by the same orthogonal 
transformation, the components W*yp are invariant. From these results it 
follows that if the length of the directors and the angles between them are fixed, 
as is the case, for example, if the directors are chosen as an orthogonal unit 
triad, and if the point X is made to traverse the X” co-ordinate curve at unit 
speed, the quantities Wx y p are the components of angular velocity of the director 
frame D, carried by X. We do not need to use (5.8), and therefore we do not 
impose the restriction Gj, =const., but this special case serves to motivate our 
calling Wp the wryness of the director frame in the undeformed material. 


6. Differential description of deformation. Results dual to those in the previous 
section, obtained by systematic interchange of majuscule and minuscule letters, 
hold also for the deformed material, but the dual wryness tensor, since it refers 
only to the relative configurations of the director frames at different points in 
the deformed material, does not afford a comparison between the deformed and 
undeformed conditions. What we wish, in the kinematic terms used above, are 
generalizations of the angular velocities of the director frame at @ relative to 
those of the director frame at X when « traverses the curve into which the path 
of X is deformed. To this end, introduce the relative wryness at x: 


Fi p= ay, pdm = ee mare dy. R= Eg Roa (6.1) 


Here, however, we encounter the quantities dé. x, where di. x is the total covariant 
derivative defined by (3.3); these director gradients appear in the theory of 
deformation of oriented bodies along with the deformation gradients x* x as 
primary local variables. In general there are 27 director gradients; when the 


directors from an orthogonal unit triad, only 9 director gradients are independent. 
Set 


Ate = dt Dk, ak = DE at (6.2) 
then 
dc AD, Dy =e, OA (6.3) 
From (6.1) follows 
Fine = Alg, pa, + Ake Wigp aM, « (6.4) 


A deformation of an oriented body is rigid if not only C=1 but also the 
directors at a may be obtained from those at X by a uniform orthogonal trans- 
formation. In this case the tensor A defined by (6.2) is a covariantly constant 
orthogonal tensor, the first term on the right-hand side of (6.4) vanishes, and 
we see that F is orthogonally equivalent to W. 

To make use of these results, we consider the anholonomic components of W 
with respect to the directors at X, the anholonomic components of F with respect 
to the directors at x: 

Wave = Wer Dy Diy a 2S Dy x; p, 
Fue pe = Fi, mp dq dg = di dyn; Ps 


(6.5) 


so that by (5.6) and (6.1) we have 


Dyx;m = Waym Dk, Gyn: m = Foo dp. (6.6) 
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If we put ga5 = Sim4_ de, then (6.4) assumes the form* 


Popp = Are p Di dap + SacG' Way. (6.7) 
Also 2Fany p= Sab; P: 


Now A is an orthogonal tensor if and only if g3, =G,,. When A is a uniform 
orthogonal tensor, (6.7) reduces to 


Fup p= Wapr- (6.8) 
Conversely, suppose (6.8) holds. From (6.5) follows 
Sige di: p= &xm Da Dé p- (6.9) 
Since this holds for all choices of a and b, we deduce 


Sab; P = Gap; P- (6.10) 
Hence 


Sap = Gap + Kap, (6.11) 


where the Ky, are constants of integration representing constant differences of 
length and angle. Thus (6.11) asserts that the lengths and angles of the two 
sets of directors differ by constants, so that the tensor A is orthogonal everywhere 
if it is orthogonal at one point. From the foregoing analysis we conclude that 
necessary and sufficient conditions for rigid deformation of an oriented body are: 


1. Cem =&xm- 
2. Foox = Wane: (6.42) 


3. At some one point, A ts orthogonal. 


The tensor W thus appears as analogous to the metric tensor g, while F is analogous 
to GREEN’S deformation tensor, C. However, we must bear in mind that while 
(6.12), is a general tensorial condition, (6.12). is not, since the anholonomic 
components F and W are calculated with respect to different frames. Cf. also 
the footnote to this page. 

The relative wryness F is thus a measure of strain of orientation, as contrasted 
with the strain of position described by the classical theory of strain. Just as 


* Since the anholonomic components of W and F are defined with respect to 
different anholonomic frames, the usual rules for manipulating anholonomic components 
do not always apply. For our purposes the particular choices (6.5) for the anholonomic 
components are essential. For example if instead we set 


K M eee 
W", p=W up Dx Dy» F%, p= Fn pt G » 


we obtain from (6.4) 
k K 
Pig ea ad Dee Wee 


Therefore a necessary and sufficient condition for A to be covariantly constant is 
[eg a 
F’ypp=W pp. 


but this is not at all the same as (6.8). 
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the Cp are certain quadratic combinations of the deformation gradients ae K> 
the F,yy are certain linear combinations of the director gradients dj.x. In 
general, the numbers of independent quantities «x, Cu, G:x, and Fypy are, 
respectively, 9, 6,°27, 27. 

Just as there are many alternative and equally correct measures of strain 
of position, so are there also other possible measures of strain of orientation 
besides Fy, 4, but we do not take up the question of characterizing the class of 
strain measures. 

By systematic interchange of majuscules and minuscules we may obtain a 
dual description in which the deformed rather than the undeformed body is the _ 
standard of reference. 

Not pausing to explore the structure whose traits have just been presented, 
we note only a formula for the gradient of the relative wryness. By differentiating 
(6.1) we get 

F* ux = aux a, + damn. x, (6.13) 
— dé. ux a, — Fry Ft K- 


Hence 
Fupm;K a Fimu;K4adg + Fi, mu (da, x Ue J di, dg. x) , 


(6.14) 
= dyn. ux4at o° (Frem Funk + Room eax — Facu Fevx)- 


The foregoing analysis is freely adapted from that of E. & F. CossERaT [1909, 
§§ 48—50]. 


B. Rods 


7. Strain of position. Let a curve @ in Euclidean three-dimensional space be 
given by the parametric equations 


PG Snip Coal (9 ated <= (7.1) 


where it is convenient to regard S as arc length. In a deformation @ is mapped 
onto a curve ¢ given by 


Kean S\s Pa, 238 (7.2) 
where 
s =s(S). (7.3) 
The stretch 4 is given by 
ds 
A= FE (7.4) 


A necessary and sufficient condition that corresponding sections of ¢ and @ have 
the same lengths is A=1. This completes the intrinsic theory of strain of 
position. 

For the imbedding theory of strain of position, we remind the reader that 
if @ and ¢ have the same curvature and torsion as functions of S, then one may 
be brought into point by point coincidence with the other by means of suitable 
rigid motion. Thus a complete description of the strain of position, as far as the 


imbedding theory is concerned, is given by the curvatures and torsions of cand ¢ 
as functions of S. 
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8. Alternative interpretation. The foregoing approach is strictly one-dimen- 
sional, but for some purposes of interpretation it is preferable to consider the 
curve@ not as isolated but rather as imbedded in a three-dimensional material and 
deformed with it so as to assume the configuration ¢. This view is legitimate 
in all cases, since, given (7.4) and (7.2) and (7.3), we may always imagine many 
deformations of a three-dimensional body such as to carry the particular curve @ 
into the particular curve ¢. From the three-dimensional theory of strain, we may 
calculate the stretch A of the unit tangent T to @, obtaining 


9% S035 TET: (8.1) 


where Cis defined by (4.2), and this must agree with (7.4). In fact, a necessary 
and sufficient condition to be satisfied by the three-dimensional deformation is 
that one point of @ be carried into the corresponding point of ¢ and that (8.1) 
be compatible with (7.4). 

In what follows, we shall shift from one interpretation to the other as seems 
helpful. For example, for a three-dimensional double field A(X, a), evaluated 
at points on @ and ¢, we have A=A(X(S), x(s)), and we may set 


~ aoe 


Aw = At Ka: (8.2) 
From the definition (3.3) follows 
A ANS OAs Pas Ane eae ie: TOTES 


and in this formula appear only quantities which can be calculated from the 
equations of the curves @ and ¢ and the deformation s=s(S). For a strictly 
one-dimensional approach, (8.3) rather than the simpler and more understandable 


formula (8.2) should be taken as the definition of A; as far as results are con- 
cerned, it makes no difference. 


In conformity with this notation, set 


K k 
dx” ck ax ds (8.4) 


VK 
sites <9 SHCA Stas AS 


I 


where X* and x* are given parametrically by (7.1) and (7.2). Note also that 
for a scalar F(S) we have simply 
dF 


9. Differential description of the undeformed rod. A rod is a curve @ at each 
point of which are assigned three linearly independent directors D,. We set 


PED GS Ga» = gxu Dt Dy, = ADO (9.1) 
so that - ‘ 
a ee OW aus Gad te 1. (9.2) 


The scalars T°, anholonomic components of the unit tangent, are the lengths 
of the projections of the unit tangent upon the reciprocal directors; the Gay 
prescribe the lengths and mutual angles of the directors; the tensor W, whose 
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anholonomic components are the rates of change of the directors resolved along 
the directions of the reciprocal directors, is the wryness of the directors along @. 
By (9.2) we see that if the G,, are constant along @, and in that case only, W 
may be regarded as an angular velocity, as explained following (5.8). Sometimes 
it is preferable to use tensor components rather than anholonomic components: 


Wey ad De Du = Dp We Diy ib = Wey De W%, = Dr Wry Dy. (9.3) 


(If @ is thought of as imbedded in a three-dimensional oriented body, the wryness 
4, along @ is related to the wryness of the body, viz. Wp, as follows: 


W%, = DL Wk, DY, Wey = Wey pX?,) (9.4) 


Suppose the quantities T*, W%, and Gyy—=Gpq be given as any functions 
of S such that (9.2), and (9.2), are satisfied and such that G is positive definite. 
Equations (9.2),,. may then be regarded as a differential system for determining 
X*(S) and DX(S). From (9.2), any solution satisfies 


Gxu Di DY = 2Gx y DE Di W)- (9.5) 


We may write (9.2), in the form Gap= 2Geiq W%); comparing this result with 
(9.5) shows that Gg y DK Df! and Gay satisfy the same first order differential 
system. By the uniqueness of solutions of such systems, (9.1), will hold for all 
S if it holds for one value of S. Given a solution X* (S), D*(S) of the system 
(9.2)1,2, Suppose (9.1), is satisfied; by (9.2), follows 


REX, =Gap T°? T= 1, (9.6) 


so the solution represents an oriented rod with S as arc length. It is easily shown 
that two sets of initial data consistent with (9.1), are related by a rigid motion, 
perhaps combined with a reflection, and that any rigid motion combined with a 
reflection carries any solution of the system (9.2),,. into another. Uniqueness 
of solutions of this system thus implies that assignment of T*, W5, and a positive 
definite Gay satisfying the conditions (9.2)3,4 determines @ and its directors to 
within a rigid motion combined with a reflection. 


10. Special cases: curvature, torsion, and twist. In the special case when the 
directors are chosen as the unit tangent, principal normal and binormal to @, 
oriented so as to form a right-handed system, the Serret-Frenet formulae assert 
that 


One 0) 
|Wo||=||—* 0 I, (10.1) 
0. Se @ 


where x and t are the curvature and the torsion of @. No such result holds 
for a general triad of directors or would be of interest here, since we seek not 
the properties of the curve itself but rather of a set of axes attached to the curve 
but turning independently of it. 


For another special case, consider a single unit director D which is normal 
to the curve @, and let it subtend an angle (S) with the principal normal N, 
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measured positively toward the binormal, so that 
D* = N* cos + BX sin (10.2) 


where B is the binormal. The quantity ¢ is the twist of D with respect to @. 
If both the curve @ and the director D are subjected to the same orthogonal 
transformation, @ is unchanged. Since Dx NX=cos @, we have 


D¥ Ny + DgN* = —sing@: (10.3) 
by (9.3), and the Serret-Frenet formula for N follows 
Wy, DM Nx + Dg (—%T* +7B*) = —sing¢. (10.4) 
Substitution of (10.2) in this result yields 
Wem) N* N™ cos @ + (Wky N¥ BM +7+¢@)sing =0. (10.5) 


Thus far in this paragraph we have spoken of but a single director D, although 
for definition of the wryness W by (9.1), a set of three directors is required. 
Since we are here interested only in the single director D, and since this is a 
unit vector, there is no loss of generality in selecting the other two directors as 
of fixed length and subtending fixed angles with D and with each other. By 
(9.2)3 follows Wx y)=0, and thus, provided m=-0, from (10.5) we obtain 

Gia Wiyr IN = BM =a. (10.6) 


That is, the twist is the excess of Wy, x N* B™ over the torsion. The scalar Wy x N* B” 
is itself an anholonomic component of W with respect to the principal frame 
of @. If we let D* be a unit vector such that D*, D, and T form a right-handed 
unit triad, since D**——N¥*sing+ BXcosq@ it follows that Wyy BX N™ 
= D** Wy D™, so that by (9.3)3 and (10.6) follows 


@ = D* DE — x. (10.7) 


This formula is the basis of the classical attempts to construct a theory of strain 
for bent and twisted rods. 

11. Criticism of the classical description of strain of a vod. Since the rotation 
of D may be prescribed in any smooth way as we traverse @, there can be no 
general connection between the twists of two different directors. In the applica- 
tions of the theory to elasticity, it is customary to think of a rod as a line furnished 
with normal cross-sections; these are represented, for the purposes of the theory, 
only by their principal axes of geometrical inertia and perhaps one or two con- 
stants such as their geometrical moments of inertia about these axes. The twist 
of the unstrained rod is then defined as the twist of either of these axes relative 
to @. Since these axes are normal to one another, a unique twist is obtained. 

Now consider the deformation of @ into ¢, defined not only by (7.1) —(7.3) 
but also by a relation setting directors d, along ¢ into correspondence with the 


directors D, along @: qa dseiD De Dl (14.1) 


Pursuing the interpretation mentioned above, we might choose for one director d 
a unit vector along the direction into which one of the directors D, considered 
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as material, is deformed according to the three-dimensional theory. Such a 
director, dy, in general would not be orthogonal to ¢. In practice, d is selected 
as the normal to the unit tangent t which lies in the plane of dy and t. The 
triply orthogonal directions so determined are called the principal torsion-flexure 
axes in the deformed rod. The twist in the deformed rod is then defined as the 
twist of d relative to ¢. 

Obviously this classical procedure is motivated only by formal simplicity 
and furnishes an inadequate description of the strain of a rod. In the first place, 
twist is defined unsymmetrically with respect to the unstrained and strained 
rods. While the twist of the unstrained rod is uniquely determined, two different 
twists can be obtained for the strained rod, depending on which of the principal 
axes of inertia of the cross-section is selected for D in the unstrained rod, and 
by the above.remarks, these two twists are in general entirely unrelated to one 
another. Moreover, the insistence that d be normal to ¢ is merely artificial and 
does not represent any kinematical requirement. 


12. Exact differential description of the strain of a rod. For a more precise 
description, we need only adapt to ¢ what has already been done for @. Since 
the operation “‘~’’ defined by (8.3) is a derivative with respect to S, not s, the 
analysis, while parallel, is not merely dual to that given above. We set 


A 


iS = dix", Cab = Sim 4a b» P=; di. (12.1) 


The Cyy are the components of deformation; the F°,, those of the relative wryness 
of ¢. 

By analogy with what was said earlicr in connection with @ it follows that 
if ", Cyy, and F*, are given functions of S subject to the conditions that Cg, 
be symmetric and positive definite and that 

Gees — 2C ele Isl §f (12.2) 
the rod ¢ is determined to within a rigid motion combined with a reflection, its 
arc length s being obtained by integrating the equation 


ds? Pes 


B= = ht = Cop bh (12.3) 


where 4 is the stretch. What has been shown, then, is summarized in the follow- 
ing fundamental theorem on the strain of a rod: Given a rod © with arc length S 
and directors D,(S), prescription of the 18 scalars t*, Cyy, and F*,, as functions of S 
subject to the aforementioned conditions determines a second rod © with arc length s 
and directors d,(s), uniquely to within a rigid motion combined with a reflection. 
In other words, the quantities ¢, Cy, and F*, furnish a complete differential 
description of the strain of a rod. It suffices to specify only the 12 scalars ?°, Cam 
and Cy» ,), since one can then calculate J, by using (12.2). 


The apparatus constructed is very general, enough so to include what would 
be regarded physically as an anisotropic rod. To represent a physically isotropic 
rod, let D, be the unit tangent to @ and d, the unit tangent to c. In this case, 
Paf=1,Pr=F=Pfs f= 0. The remaining two directors, both along @ and 
along ¢, may be assigned arbitrarily. 
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13. Rotation.. Again we consider a general director frame. To determine the 
finite rotation of a rod, it is necessary to regard the rod from the Euclidean space 
in which it is imbedded. So as to motivate the definition shortly to be given, we 
first consider @ to be a material curve ina three-dimensional body and the directors 
D, material elements, this being the case that one-dimensional theories of rods 
are intended to idealize. Then, as remarked in § 5, in (6.3) we have A*,=x* x, 
a”,—X*,. By the classical theory mentioned in § 4, we can decompose x* x 
uniquely into a stretch and a rotation, and the rotation so determined is the 
rotation of the principal axes of strain at X into the principal axes of strain at a. 
In a purely one-dimensional theory, we do not have available the «*,, but from 
the two sets of directors we can define A’, by (6.2), and, since ‘the linear inde- 
pendence of the directors assures us that A*% is non-singular, from the polar 
decomposition theorem (§ 4) we can write 


Al = Rhy PK = pn R",, (13.1) 


where R is a rotation tensor and P and p are positive definite symmetric tensors. 
The tensor R then represents the local rotation of the rod. What we have proved 
regarding it is summarized as follows: Given two rods with directors D, and d,, 
a unique local rotation is defined by (13.1); the rotation is independent of the choice 
of the directors to this extent, that if we imagine and then leave fixed a deformation 
of a three-dimensional body in which © ts a material line and the given directors 
D, are carried materially into the given directors d,, then we may choose as directors 
any other set of linearly independent material vectors and obtain the same rotation. 


When the directors are chosen as above, (12.1), becomes 
Cm = Cra DD ici. (13.2) 


The quantities C,, are anholonomic components of the three-dimensional deforma- 
tion tensor C with respect to the directors of @. 


14. Special case: Hay’s formalism. The only previous analysis approaching 
the generality of what has just been presented is that of Hay [1942, §§ 2—3]. 
In effect, he takes the directors D, of the undeformed rod as any unit triad such 
that D, is the unit tangent to @, and he chooses co-ordinate systems such that 
not only D¥ = 6% but also dé = 06%. Hence (6.2) yields A*,= 6k, a%,= 6%, but 
it must be remembered that the co-ordinates at # are not generally orthogonal. 
In these co-ordinates we have from (12.1). 


Ca Cin nde — btm Ok 4s (14.1) 


thus the components of the metric tensor at # are numerically equal to the com- 
ponents of stretch, this being a generalization of Hay’s equation (3.9). From 
(12.1), we get in these co-ordinates 


dx? ds 


ad, k Mm. 
tea tne} % ds|dS’ 
= OE {hi} OF 4, (14.2) 


0g d 
__ sa kp ole 1 Spm 
Of Ob & ( ay” eon us 


PY = ay 
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where the rod ¢ is taken as the x!-curve with x1=s. The skew-symmetric part 
of this equation is essentially Hay’s equation (3.5); hence our relative wryness F 
includes and generalizes Hay’s rotation vector w. 


C. Shells 
15. Intrinsic theory of strain of position. Let the undeformed surface / be 


given by ye 
Ay Ras XK (V1, V2) Ses (Vy; (15.1) 


the deformed surface s by 
xt = x (yl, v2) = x*(v'), (15.2) 


where, as in the rest of this section, Latin indices run from 1 to 3, Greek indices 
from 1 to 2. The deformation is specified by the functional forms of the right- 
hand sides of (15.1) and (15.2), augmented by 


P= e(V); (15.3) 
For full generality, it is sufficient but not necessary to take 
vi = 62 V%. (15.4) 


The surface metric tensors are given by 


OS = AzcaV4 IVa =A gs dv dv, 


; , (15.5) 

ds* =a5,dv° dv* = aycdV4 dV", 

where E 
Age =AyaVis Vie, Gaz = MeV ais, (15.6) 
and where “;’’ denotes a partial derivative. If A(V) and a(wv) are known, and 


if the parameter transformation (15.3) is given, we may calculate for any material 
element the changes of length and angle occasioned by the deformation. Also, 
we may calculate the total curvatures R,.,. of Y and s. Corresponding to any 
assigned real symmetric tensor a(v), with A(V) and the relation (15.3) regarded 
as given, there exists a surface s into which is deformed, and any two such 
surfaces s, and s, are applicable. In principle, this completes the intrinsic theory 
of strain of position of a shell. 
In terms of the three-dimensional equations (15.1) and (15.2), we may cal- 
culate the induced surface matrices A and a, obtaining 
A y= 8xm <9 XY, Ags = Sim Xi wr (15.7) 


Voie 


16. Imbedding theory of strain of position. The second fundamental tensors 
B and b are given by 


Bsy= Nx Xe, Dse male He oak (16.1) 


where X¥ye= (X14), 2, "y¢=(x%),¢, the total covariant derivative “;” being 
defined by (3.3), and where N and n are the unit normals to “and s. We have 


ByzdV4dVi = Bye dv? dv’, bse dv° dv’ = bygdV4 ay. (16.2) 
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where 7s 
Bse= Baz V4 i baz = d5¢ wy Ug. (16.3) 


The tensors A and B are related by the equations of Gauss and MainarDI- 
CODAZZI: 


SAS Ee, Basie Bug = 0; (16.4) 


where K is the Gaussian curvature, d=detA,-, B=detB,-. There are dual 
relations connecting @ and b. When ¥ and a parameter mapping (15.3) are 
given, assignment of arbitrary symmetric tensors a and b satisfying the duals 
to (16.4) determines a surface s into which Y is deformed, uniquely to within 
a rigid motion *. 

The surface Y is obtained by solving the equations 


Xz = Bz, Nes Ni aA By Xf, (16.5) 


using appropriate assignments of X* and N* at one point, and s is obtained by 
solving the dual equations. In principle, this completes the imbedding theory 
of strain of position. 

The normal curvature K,y) of Y in the direction of a unit vector MN is cal- 
culated from 
Bz, N*=N4 ; 
AgyN?N* ’ 


the principal curvatures K, and Ky are the curvatures K,y) in the principal direc- 
tions of B, and the Gaussian curvature K and the mean curvature K are related 
to these principal curvatures through 


aha ee R= Roto A By. (16.7) 
In approximate theories of elastic shells it is customary to take the changes 
of principal curvature as measures of strain. To do this in the exact theory would 
be to introduce complications analogous to those resulting from using the dis- 
placement vector and the strain tensors in finite strain of three-dimensional 
bodies. 

17. Differential description of the undeformed shell. A shell is a surface S 
at each point of which are assigned three linearly independent directors D,. 
The theory of strain of orientation for a shell may be constructed in analogy 


to what was done for three-dimensional bodies in §5. Set 


Nee OI Gin == Gra DAD, = We == De DES (17.1) 
Then , , 
Xa De Aas Dia= Ds W{a. (17.2) 
Then by (15.7), follows 
Aje = Gigxixe, (17.3) 
while by (15.8), follows 
Baz = Nx DE (X4.2 + Wz X4). (17.4) 


* Cf., e.g., EISENHART [1940, § 39]. 
Arch, Rational Mech. Anal., Vol. 4 22 
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From the relation (17.2), 
Xe — DK (X8.,4+ Ws X%), (17.5) 
Disa Dé (W5.4 W e+ We,4)- (17.6) 
We also have the integrability conditions 


a2 X* K eDa 


——~ soo 2 -——s K — — . 
sr/tA 7 ela ce 2. ev4oev= Pa; az] = 0, pk 
whence follows axe 
{4 Wao X% =0, (17.8) 
ev) 
c b c b | OW aca 
Woe Wee =e eee el = 0. (17.9) 
Brome (7A)\ande (A722) or 
U7-4)e 17-2)s Gap, 4 = 2Weapya- (17.10) 


The equations (17.2) may be regarded as a differential system for the equation 
X (V) of the shell and for the assignment of directors upon it. When X%, Gay = Goa, 
and W®, , are prescribed as functions of V subject to the conditions (17.8), (17.9), 
and (17.10), the system (17.2) is completely integrable. Locally there will then 
exist a unique solution X(V) and D,(V) taking on prescribed values at one value 
V° of V. Provided that not all the quantities X%, X%, vanish, that G, be positive 
definite and that the prescription of D,(V°) is consistent with (17.1),, the solution 
represents an oriented shell and the relation (17.1), holds. By using the uniqueness 
of solution, it is a simple matter to show that X45, G,,, and W®,, determine S 
and its directors to within a rigid displacement combined with a reflection. 


18. Exact differential description of the strain of a shell. At the points x of the 
deformed shell s, assign the directors d,, and, as suggested by (12.1), put 
to a ay x gv", Con Camda ses (18.1) 
Tea: dy. 4 = dy a‘. 5 wns 
Pea hy adh V4, | dg = ORS V6. (18.2) 
The Cy) are the components of deformation; F%, 4 is the relative wryness of s; 
the x, are certain tangent vectors. Then 
Gate Can ead ea ae 
bog = mpdt [FT pa Vg} as Vento Vesa ou Wane 


Thus we see that knowledge of », x*,, C,,, F4,,, and d* as functions of V suffices 
to determine the first and second forms of the deformed shell s. 


By analogy with (17.8)—(17.10) we have 


omy 7b a 

aps apstay = 0, (18.4) 
cs oF* 
FoePaa—FoaFie+2 spa = 0, (18.5) 


Cap; 4 = 2 Cea Liat (18.6) 
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Further, when v, F4z, «4, and Cy, are given, the requirements (18.4) — (18.6), 
subject to |v z|=:0, and the conditions that C,, be symmetric and positive 
definite and that not all the quantities x{,x2, vanish, the differential system 
(18.2) determines an oriented shell satisfying (18.1), uniquely to within a rigid 
motion combined with a reflection. What has been shown, then, is summarized 
in the following fundamental theorem on the strain of a shell: Given a shell S 
with fundamental forms A(V) and B(V) and with directors D,(V), prescription 
of the 32 quantities v°, F%, 2, x*, and Cay as functions of V subject to the aforemen- 
tioned conditions determines a second shell s with equation «= a(v) and directors 
d,(v), uniquely to within a rigid displacement combined with a reflection. In other 
words, the quantities v°, F4z, x4, and C,, furnish a complete differential de- 
scription of the strain of a shell. 

A theory of rotation is easily constructed by analogy to what was done for 
rods in § 13. 


19. Resolution into normal and tangential components. Since X%,, X*,, and 
N* are three linearly independent vectors, we may write 


DS DEX = DNs, (19.1) 
where 
Dee AED Kes; D, = D* N,. (19.2) 
Then by (16.5) follows 
Df 4 = (Dg.4— Dy BG) X¥z + (DE Bay + Dy, 1) N*, (19.3) 
whence : y 
Nx DE =Di Baa t Daas (19.4) 
Xx.eDi 4 = Diz.a—D, Bey. (19.5) 


From the fact that B,; and A, , determine X(V) to within a rigid motion and 
that DX (V) is uniquely determined by DZ and D, when X(V) is known, it follows 
that B,z, Ayz, D4, and D, determine Y and its directors to within a rigid motion. 
Other than (16.4), there are no compatibility conditions to be satisfied by these 
quantities. 

It involves no restriction to require that s and one director not tangent to s, 
say d,, be material with respect to an unspecified three-dimensional deformation. 
That is, if x(v), d,(v), v(V), X(V), and D,(V) are given subject to the condition 
that d, be not tangent to s and D, not tangent to Y, there will exist infinitely 
many mappings «#(X) such that 


x(v) =a(X(V(v))), X(V) =X(x(v(V))). (19.6) 
nd 
: dj = *;x Df, (19.7) 


the quantities x*,« being uniquely determined as functions of V by (19.7) and 


Since the remaining directors can be assigned arbitrarily, they will not necessarily 
be material with respect to this deformation. 
22% 
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20. Special case: SYNGE & CHIEN’S formalism. The formalism just given is 
easily specialized so as to give results depending on particular choices of co- 
ordinates for the strain of position. We consider here only the most general 
form of the usual approach, that given by SyNGE & CHIEN [1941, 7]. The 
intended interpretation is that s and d, are material, so that (19.6)—(19.8) 
apply, with D, varying with the deformation #(X) in such a way that d;—n. 
Choose co-ordinates and parameters such that 


2 (X) =X. XV) (VT7, 0). ee) ae oF), 


(20.1) 
Ses = Ona- 
Then 
th = nt = m= bk DE = 08, Ny = by a(G*)-3, (20.2) 
and 
Azyg=Ggy 056%, 45 = Sim 6: OF. (20.3) 
From (16.1), (20.1),, (20.2)3, and (20.3), 
aax® 
2Bs4 = 2Ng( 2 + {lle} XM XT + (du) X44), 
Pie Gss {yx} of On, 
20.4) 
G oG ( 
— (G33 ~4( 2620 9b OGKs 5X 2 Wh ‘) 
( ) avi ae a av= 4 ox Oz é a 
+ (G38)- br 98 (“5 EA en cet 
Caan 
which we may write in the form 
ays 8804 Ga be + Gye 44 
(20.5) 


@Aga , 0OAga  OAzy4 
i] 


av4 ave ava 


ae (GAR) CK og | = 2(G33)-3 Bey. 


This is essentially SyNcr & CHIEN’s equation (69). Cf. also Cnren [1944, 
eq. (6.13) ]. 
From (19.2) and (20.2), 
De Gree. Dis = (Gay (20.6) 
Now, using (20.1),, (20.2),, and (20.6), we obtain 


2Xx.eDF. 4 = 267 Xs (22s. 


+ {fp} XM, Di), 
= 2G, x OL {ap} 04 03, (20.7) 


_ 8G7K J ok G73 I OGK3 6K, 
OXt Goeth Mav eat Ob Orem 


Rods and Shells 315 


From (19.5), (20.6), and (20.7) we have 


oG 

K K 

2Xx;( Dr, 4) = ae 6£ 6% = Dis. 4 Pits Drs), 
ODjz , D1, Zi 

— = 2D 4 —2D,Bz ) 20. 
Bee oG ~|0Ag 0A OAgs 

13 gf += “2 8) sia ve! Ate( ez 7 04a ye “A 3(G28)-1 Be 

ev4 pect hale: OV) ove 2 ey2 en ee: 


which is similar to, but not identical with (20.5). The apparent discrepancy is 
resolved by noting that the equations 


A SARA 01, 0 = 0, = GG, 0) = G'G,, OG bf Ass, (00) 


imply that 
(Git Ge 62 2 Gy 7 6), At. (20.10) 


For s, we have by analogy with (20.5) 


or, (20.11) 
SYNGE & CHIEN found it sufficient to introduce nine measures of strain. The 
set ayg, 0.8, Di4, and D, is equivalent to that which they used. These satisfy 
three compatibility conditions which may be taken as the dual* of (16.4). 
When Y is given, these quantities determine D, uniquely and s to-within a rigid 
motion. With the conventions adopted above, part of the problem involves 
determining a tensor g,,, consistent with the conditions (20.3),, (20.11), and the 
fact that the Riemann tensor based on g,,, must vanish; such complexity is the 
price one pays for eliminating the functions x(v) as unknowns. 


Part III. Theory of stress 


A. Rods 


21. Stress principle for rods. A rod will be regarded here simply as a curve ¢ 
which may be the seat of dynamical actions. The basic postulate is the stress 
principle: At each point on a rod, the action of the material to one side upon the 
material to the other is equipollent to that of a stress resultant vector S and a couple 
resultant M. If we denote by S, and M, the values of the vectors S and M 
appropriate to the action of the material to one side of a point s on that to the 
other, and by S_ and M_ the values of S and M when the roles of the two sides 
are interchanged, then the balance of forces and moments requires that 


S722 Sac Mi Ma (21.4) 


Thus, with an appropriate convention of sign, we may use the vectors S and M 
without further designation. 


* These equations are equivalent to the three given by SyNGE & CHIEN [1941, 
Eggs. (38), (39)]. 
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22. Differential equations of equilibrium in space. Let F and L, assumed to 
be integrable functions of s, denote the assigned force and couple per unit length 
acting upon the rod. Then for points where S and M are differentiable and 
F and L are bounded, a local statement of the principle of equilibrium is 

Y ith, alae: (22.1) 
M+pxS+L+pxFf=0, 


where p is the position vector to s from the origin of some inertial frame and 
where “‘~’’ denotes the intrinsic derivative, defined as in (8.3), except that s, 
the arc length in the actual configuration of the rod, is the independent variable*. 
By substituting (22.1), into (22.1), we may eliminate F and replace (22.1), by 
the equivalent condition 


M+txS+L=0, (22.2) 


where € is the unit tangent. At points where S and M are not differentiable, 
(22.1) and (22.2) need not hold, but application of the principle of equilibrium 
shows that S and M must be continuous if F and L are bounded**. 


23. Resolution into normal and tangential components. Since (22.1) and (22.2) 
are in vectorial form, they are valid in an arbitrary curvilinear co-ordinate 
system. It is customary, however, to refer them to a particular frame defined 
with respect to the rod c. Retaining full generality at the start, in the scheme 
of § 12 we assign three linearly independent directors d, and reciprocal directors 
d* to c. In general three-dimensional co-ordinates, let S* and F* be the contra- 
variant components of S and F, M, and L, the covariant components of M and 
L, and define corresponding anholonomic components: 


St=ds*, M,=d*M,, 


Fo=@F LL, =d*L,. (23.4) 
From (22.1), and the result dual*** to the reciprocal of (9.1); we have 
Be St aide eek 
Seah ie oS (232) 
where w is the wryness of the directors along ¢. Similarly 
ee aM. == (Cxpqt? S'+L,) + d@w*,, M,. (23.3) 


esis) theadualeo teste 


** If concentrated loads are present, they determine jumps in S and M: 
[S]=—,° [M]=—,. 


*** Note that w is not the relative wryness F of § 12 but rather the dual of the 
wryness W of § 9. 
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Hence the statical equations in anholonomic components are* 


oe + wi, 4+ F*=0, 
dM, z (23.4) 
ee gee a My + eapet? SS+ L,= 0. 
Thus far the director frame has been arbirary. We now require it to be a 
unit orthogonal triad such that d;=#, the unit tangent. By the dual of (9.2), 
the wryness w then satisfies w,, = — Wp, and may be interpreted as an angular 
velocity **. The component S‘, which is the projection of S onto the tangent to ¢, 
is called the specific tension of the rod; the components S? and S?, the specific 
shearing forces; M’, the specific twisting couple; M? and M®, the specific bending 
couples. This special choice of directors, while not simplifying (23.4),, implies 
that t’=1, #=??=0 and hence reduces the three components of (23.4), to the 
following explicit forms: 


dM, 

Ys ei My ar L, = 0; 

am — wy My — S?+ L,=0, (23.5) 
dM: 


2 — w®,M,+ S?+L;=0. 


ds 


B. Shells 


24. Stress principle for shells. A shell will be regarded here simply as a sur- 
face s which may be the seat of dynamical actions. The stress principle for 
shells asserts that the action of the part of the shell outside any imagined closed 
curve ¢ on the part inside is equipollent to a field of stress resultant vectors Sim) 
and couple resultant vectors Mq,) defined on c. The subscript *** n refers to the 


* By definition y 
C55 a Vesti aden 


where ¢;;; is the permutation symbol such that ¢.;= +1. Hence 


is le k aP 79 
Canc Tt VS Expq da dy de, 


=+ V2 £5, det dk : 


Now 
g(det d®)?= det (g,,, a d;’) = det go;, 
hence 


eabe Vdet g,¢ €abe? 


where the sign is so selected so as to agree with that of det di. The quantity det Set 
has been evaluated in the footnote on p. 301. In particular, for directors forming a 
right-handed orthogonal unit triad we thus obtain 


“ape — “abe: 

*x* The classical notation for the component w,, is tT or — t; the other two inde- 
pendent components are written as + and +x’. In the classical treatments it is 
not always made clear that these quantities refer to the loaded rod. 

*xkk This n is not to be confused with the unit normal to s, which was denoted by 
n in §§ 16—20, as it will again from § 26 onward. 
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unit outward normal to ¢; of course, n is a vector field defined intrinsically in 
the surface s, but Sj) and Mi) are fields defined in the three-dimensional 
Euclidean space in which s is embedded. A mathematical expression of the 
postulated principle, for the case of equilibrium, is 


c 


§ (Min) + PXSm) ds +f (L+pxF)da=0, 


¢ Sj 


(24.1) 


where s is arc length along ¢, where s, is the portion of the surface s inclosed 
by ¢, where p is the three-dimensional position vector to the running point of 
integration, and where F and FI are the three-dimensional fields of assigned force 
and assigned couple. The vector integrals in (24.1) are to be understood as 
abbreviations for the integrals of rectangular Cartesian components. 

A classical argument, not different in principle from its counterpart for three- 
dimensional bodies, yields not only 


Sa =—S-m, May = —M-w (24.2) 
pee Shy = Sn, Miy = Mn, a 


where, as in all that follows, Greek indices have the range 1, 2. The quantities 1, 
are the covariant components of the unit normal to ¢ in any curvilinear co-ordinate 
system v1,v? on s. By hypothesis, the quantities S,,) are absolute vectors 
and the Mi) are axial vectors; while to derive (24.3) rectangular Cartesian 
spatial co-ordinates were adopted, the results are double tensor equations in the 
sense defined in §3 and hence are valid in all spatial co-ordinates. The double 
tensors S*° and M*® are the fields of stress resultants and stress couples. There 
are six components S*® and six components M*°. 

In the classical treatments of the theory of shells, the vectors LZ and M are 
assumed tangent to s; in this case the number of independent components M*?® 
is reduced from six to four. (Roughly speaking, this assumption corresponds to 
the case when the assigned couples result only from forces acting a finite distance 
from s, not from internal couples such as occur in polarized media.) In mathe- 
matical form this assumption asserts the existence of surface tensors L® and 

£6 
M*° such that Poe, Min H Mes, (24.4) 
where x*,5=0x*/0v°, a=a(v) being an equation of s in general co-ordinates. 


25. Differential equations of equilibrium in space. Again we suppose the space 
co-ordinates rectangular Cartesian, we consider S and M as functions of w only, 
and we substitute (24.3) into (24.1). In regions where S§ and M are continuously 
differentiable, a classical argument based on GREEN’s transformation yields the 
differential equations 

q Ss’? FP = 0, 


Mires, 4 alt, Selo 4 Th? = 9, (25:4) 


where M*?° and L*? are absolute alternating tensors equivalent to the axial 
vectors M*° and L*, and where the subscript comma denotes the covariant 
derivative with respect to the surface metric a, except that z*5=02*/dv°, 
#—2(v) being a rectangular Cartesian equation of the surface s. 
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Now consider the equations 


where the space co-ordinates and the surface co-ordinates are arbitrary in- 
dependently selected general curvilinear systems, where «=a:(v) is an equation 
of s referred to these two systems, where x*,,;=0x*/dv°, and where other occur- 
rences of the semicolon denote the total covariant derivative (3.3). The equations 
(25.2) are in double tensor form; when the space co-ordinates. are rectangular 
Cartesian, (25.2) reduce to (25.1), already proved valid in such systems; therefore 
(25.2) are the general differential equations of equilibrium for shells. We may 
continue to regard S and M as functions of v only, or we may consider them 
to be functions of & also, as we please. 


The elegant simplicity of this derivation should not conceal the complexity 

of the result. When (25.2), is written out, it assumes the form 

Oss bY yp como E\ cho k 

~ 8 + tno} Ba S + {5} S¥? + F* =0, (25.3) 
where {,,} and {,<} are Christoffel symbols based upon the space metric g and 
the surface metric a, respectively, the two metrics being related as usual: 
Ass = 8pm Xi5 Xe. When (25.2). is written explicitly in terms of the axial vector 
M*°, it assumes the form 


amré 
av® 


+ {mp} ho M? + {ie} MPO + ef g 2%, SP? + LF —0, (25.4) 


where e¥,, , =" VE Esmp» Esmp being the permutation symbol such that ¢,93= +1. 
These formulae involve doubly contravariant tensors considered as functions of 
only. In terms of physical components or of components allowed to depend on «x 
as well, they would take on still more complicated forms. The results (25.3) and 
(25.4), specialized by means of the assumption (24.4), are the differential equations 
of equilibrium in the general form first derived by SYNGE & CHIEN [1941, 
pp. 104—111]. 


26. Resolution into normal and tangential components. Resolving the variables 
occurring in (25.2), we may write 


Be a Eg”, |B fern Bae os a ER 


206.1 
Sie et al Sloe ee te ae nk, cae 


where n is the unit normal to s. The following table connects the components 
occurring in (26.1) with the terms usually employed in shell theory: 


F, L = normal components of specific applied force and couple. 
F°, L° = specific applied force and couple tangent to the shell. 
S*® = cross force resultant. 
S’® — membrane stress resultant. 
M° = cross moment resultant. 
M°®’ = couple resultant. 
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In this notation, the classical assumption (24.4) takes the form 
D0 Ve ae (26.2) 


but this we do not at once adopt. The normal and shear components of S?? am 
an orthogonal co-ordinate system are called normal and shear membrane stress 
resultants; the normal and shear components of M’® in such a system are called 
twisting and bending couple resultants, respectively. 

To express the equations of equilibrium in terms of tangential and normal 
components*, we use the identities 


k (Me Aa) ee q 
C9 My = Cgnp Koy, Cap X90 =a” Cig My, 


(26.3) 
Ky 6 = by. n, Wy = — bya 


Go a 
where b is the second fundamental form of s. From (26.1) and these identities 


follows 
SiS = (S725, Nt Se | Wis, = [Sia + b,5 S?°) n*. (26.4) 


Substituting this result and (26.1), into (25.2), yields 
(SE nial boo SS = FS) xt y ae (S*.4 $3 b,,5 Se ot F) ni 0" (26.5) 


Taking the scalar product of this equation by n yields the condition for equilibrium 
of normal forces: 


35+ b,5S’°+F =0; (26.6) 
taking the vector product by n, the conditions for equilibrium of tangential forces: 
SY OO, S? ee, (26.7) 


Similar resolution of (25.2), yields the condition for equilibrium of twisting 
moments : 


M’., + b,5M+ 6, SL = 0 (26.8) 
and the conditions for equilibrium of bending moments: 
M”?.,— a”? b,5 M® 4+ aves, S°+ L’ =0. (26.9) 


In these formulae it is legitimate and natural to regard all fields as functions of 
the surface co-ordinates © only; in this interpretation, the double covariant 
derivatives “;’’ reduce to “‘,”, the usual covariant differentiation based in the 
surface metric a. The fully general equations (26.6)—(26.9), but in rectangular 
Cartesian spatial co-ordinates, were first given by E.& F. CossERAT [1909, 1, 
§§ 35—37], who derived forms in material co-ordinates as well. Cf. also HEuN 
[1913, § 20]. 

Under the classical assumptions (26.2) the total number of independent 
components of S and M is reduced from 12 to 10, and the equations (26.8) and 


* This resolution is effected by Syncr & CHIEN [1941, p. 109] by use of a special 
co-ordinate system. 


Rods and Shells 324 


(26.9) for equilibrium of moments reduce to 
bys Vie si Cy6 S’?= 0, 
Me a? oe, oS? eo TY 0: 


? 


(26.10) 


The first of these is a linear algebraic equation expressing the difference of shear 
resultants, S!#]) as a linear combination of the four couple resultants M”®. 

Further specializations appropriate to special co-ordinate systems on s or to 
special dynamical assumptions occupy much of the literature and are not dis- 
cussed here. 


27. Remarks on the relation of stress and couple resultants to stresses, and on 
the equations of motion. It is also legitimate to treat shells and rods as three- 
dimensional bodies and to regard the stress resultants and couple resultants as 
derivative quantities, defined in terms of the three-dimensional stress tensor as 
certain fields existing only on the line or surface but equipollent there to the field 
of internal stress in the whole shell or rod. For the theory of plates, this view 
goes back to Caucuy and Potsson, but for curved shells it was first carried 
through exactly and in some generality by NovozuiLov [1943, 1, §1] (cf. also 
NovozHILov & FINKELSTEIN [1943, 2, §§ 1, 4]), by TRUESDELL [1945, § 8], and 
by CHIEN [1948]*. We do not present a derivation by this method because such 
a derivation in general co-ordinates has been given ZERNA [1949, § 3] (cf. also 
GREEN & ZERNA [1950, §3] [1954, 2, §10.2]). The method consists in in- 
tegrating the three-dimensional equations across the shell; of course, no approxi- 
mation is involved. The result are formally identical with (26.6), (26.7), and 
(26.10). 

The interest in a derivation of this kind, beyond its great unifying value, is 
twofold. First, the applied force and couple F and L, which are necessarily 
quantities taken a priori in the direct theory given in this paper, appear as 
quantities defined in terms not only of the applied internal force but also the 
applied surface load, and in an exact theory the thickness and curvature of the 
shell have their influence on the net effect of these loads as it appears in the fields 
F and L defined on s. 

Second, this method enables us to obtain equations of motion, since all that 
is required is to replace the three-dimensional applied force f per unit mass by 
f—a, where a is the three-dimenstonal acceleration. No such simple device can 
be applied directly to the shell or rod; while, indeed, an element of the shell 
or rod is endowed with momentum, in an exact theory this is generally not 
proportional to the velocity of any point on the shell or rod. This distinction, 
if properly pondered, will convince the reader (if he has any doubts) that in 
truth, occasional conversation to the contrary, at bottom we always conceive 
the world as having three dimensions at least. Without such a conception, the 
distinction just made would not have meaning, but without this distinction we 
should be unable, in principle, to determine the exact motion of points on a shell 
from the exact dynamical equations. 


x All this work appears to date from the period 1942—1944 and to be done 
independently. The proper definitions of stress resultants and couple resultants in 
terms of the three-dimensional stress tensor are due to Love [1893, 2, § 339]; the 
essential idea was given by Lams [1890, 1, § 2] and BassET [1890, 2, §§ 5, 18]. 
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For rods, an exact derivation based on the three-dimensional equations of 
equilibrium has never been worked out*. It should be equally possible to base 
such a derivation on the exact equations for shells. 
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On the Riemann-Green Function 


E. T. COPSON 


Communicated by A. ERDELYI 


§ 1. Introduction 


The first general solution of the problem of Caucny for an extensive class 
of partial differential equations was given by RIEMANN almost a century ago 
in his well-known paper on the propagation of sound waves of finite amplitude *. 
Although stated only for certain special equations, it is applicable to any linear 
equation of hyperbolic type of the second order in two independent variables; 
it depends ultimately on finding a certain subsidiary function, often called the 
Riemann-Green function, which is the solution of a characteristic boundary value 
problem for the adjoint equation. RIEMANN gave explicit formulae for this 
subsidiary function in two cases of special importance in gas dynamics. 

Although a full account of the method in its general form was given by DaR- 
BOUX**, little progress has been made, and this seems to be due to the difficulty 
of finding the Riemann-Green function. The books often content themselves 
either with stating and verifying an expression for the function in one or two 
cases or with finding it in special cases by an inspired guess as to its form. It 
seems then to be worth while reviewing the present position in the hope of making 
RIEMANN’s method more generally useful. 


§ 2. A description of Riemann’s method 
A linear hyperbolic equation of the second order in two independent variables 
may be put into the form 


_8U  8U 


(DAG See ( oU 
On oy? 


Ox 
by choosing the coordinates x and y so that the two families of characteristics 
are x + y=constant; the coefficients a, b, c are functions of x and y alone. The 
problem of Cauchy is to find a solution U, given the values of U and its first 
derivatives on a certain curve C which has the property that no characteristic 
cuts it in more than one point. Under certain continuity conditions, the problem 
has a unique solution, given by 


OXY) =2lUV Jat e[UV a+ 
+2 SU(VU, — UV, + 20 UV) dx + (VU, — UV, + 2a UY) dy} 


2a 2b a cU=0 (2.1) 


(2.2) 


* RreEMANN: Abh. d. Kén. Ges. der Wiss. zu Gottingen 8 (1860), reprinted 
in Collected Works of Bernard Riemann, pp. 156—175. Dover Press 1953. 
** DarpBoux: Lecons sur la Théorie des Surfaces, vol. II, pp. 71—111. Paris 1915. 
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where A and B are the points in which C is cut by the characteristics through 
P(X,Y); V=V(x,y; X,Y), the Riemann-Green function, is the solution of 
the adjoint equation 


o2V eV a) a) 
LEV faa y3 2D = (aV) +2 os (6V)+cV=0O (2.3) 
such that 
OV OV 
a me ae on y—x=Y-—X, 
OV oV 
ee re on ytx=Y+X, (2.4) 
V=1 ate (XP ¥)% 


It is readily seen that 


[fis (VU, —UV, + 2aUV) — = VU, — UY, +26 UV)} dx dy 


vanishes for any region of integration. If we take the region to be that bounded 
by the characteristics PA, PB and the arc AB of the curve C, equation 
(2.2) follows at once by GREEN’s transformation. 
(see. Fig.) 

The proof of the existence of this Riemann- 
Green function V(x, y; X, Y) is omitted as we shall 
be concerned with its actual construction in certain 
special cases. For the moment, we note the follow- 
ing properties. 


(*) V does not depend on the choice of the 
curve C. 

(17) If we change the dependent variable in (2.1) from U to U, where U,= 
g(x, y) U, the Riemann-Green function for the transformed equation is 


Fig. 14. Riemann Problem 


x 
Vi (x, 9; X,Y) = ee 


(77) The function U(x, y; X,Y) =V(X,Y; x, y) is the Riemann-Green func- 
tion of the adjoint equation, that is gua function of (x, y), it satisfies the equation 
LU =0 and the conditions 


VAG VESEY a (2:5) 


aur 8U = Ae eee 
oa eae (a+b)U on y—*%=Y—X, 
0U oU 
= = 
aa By (a — b)U on y+%x +X, 
C= ate (erys)2 


Sometimes it is more convenient to use the characteristic variables 
(P= 00 == WV, Ve 5 GO 


The standard form of the equation is then 


tL c’u=0 (2.6) 
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where a’, b’, c’ are functions of y and s. In terms of these variables, RIEMANN’S 
formula for the solution of the problem of Cauchy is 
u(R,S) = 3 [wvja +3 [4 0]p — 


r , (277) 
—£ f{(vu,— wv, + 26 uv) dr—(vu,— uv, + 2a uv) ds} 
AB 


where A and B are the points in which C is cut by the characteristics y= R and 
s=S respectively through P(R,S); v=v(r,s; R,S) is the Riemann-Green func- 
tion, the solution of the adjoint equation 

o2u 4) ( ’ ) é (bv) eee (2.8) 


ry = 


ov Os Or Os 
such that ee 
pee aera} On hia 
Os 
Les on ‘s= 5S (2.9) 
or ? 
uA at (R, S) 


(See Fig? 2.) 
The two definitions are equivalent, in that V 
Fig. 2, Riemann Problem is transformed into v by the change to charac- 
teristic variables, and conversely. 
So far as I know, six ways have been used to find the Riemann-Green func- 
tion for particular types of hyperbolic equations. These will be discussed in the 
following order. 


(*) RIEMANN’S original method was based on the fact that the Riemann- 
Green function does not depend in any way on the curve carrying the Cauchy 
data. If it is possible to solve by some other means the Problem of Cauchy for 
a special curve C depending on one variable parameter, a comparison of the two 
solutions should give the Riemann-Green function. In the case of the two equa- 
tions considered by RIEMANN, it was possible to solve the Problem of Cauchy 
by a FourRIER cosine transform with Cauchy data on a straight line. 


(77) HADAMARD pointed out* that the coefficient of the logarithmic term in 
his elementary solution is the Riemann-Green function of the adjoint equation. 
It is possible to modify HADAMARD’s construction so as to give both functions at 
the same time. 


(zz) It is easy to construct an integral equation whose unique solution is the 
Riemann-Green function. 


(1v) CHAUNDY, in his work on partial differential equations of hypergeometric 
type, was able to construct the Riemann-Green function by the use of symbolic 
operators and power series. This work appears to be little known. 


(v) A. G.Macxre has constructed complex integral solutions of certain 
equations. Such a complex integral gives the Riemann-Green function for an 
appropriate choice of contour. To some extent, MACKIE was anticipated by 
CHAUNDY, whose approach was rather different. 


* HapaMARD: Lectures on Cauchy’s Problem in Linear Partial Differential 
Equations, p. 72. Dover Press 1952. 
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(vt) TITCHMARSH gave a direct solution of the characteristic boundary value 
problem defining V for the equation of damped waves by means of a complex 
Fourier integral. 

§ 3. The method used by Riemann 

Let us consider the Problem of Cauchy for the equation 


eU e@U ,. aU aU . 
Ae ay? r 2 a 2b By eu (3.4) 
under the conditions aU 
(oii), Te F(x) (3.2) 


when y= Yo, where yo is an arbitrary constant. The solution (2.2) then reduces to* 


Ea oS Vo 


U(X,Y)=2 f V(x, ¥9;X,Y) F(x) dx. (3.3) 
A oe Vo 
| 
8 
yt 
PI) P(X,Y) 
A 
A B 
ye L 
Fig. 3. Diagram for Eq. (3.3) Fig. 4. Diagram for Eq. (3.5) 


(See Fig. 3.) If we could solve this problem by some other method, a compari- 
son of the two solutions would give V(x, yy; X, Y) when x lies between X + (Y— yp) ; 
aS Yo is arbitrary, this would give V(x, y; X, Y) whenever X — x lies between 
+(Y—vy). Similarly if the data were 


aU 
U=0, $=) (3.4) 
x 
when % = %,, where %, is arbitrary, we should have 
Y+X—x%5 
U(X,Y)=3 J V(%o, ¥; X, Y) Gy) dy. (3.5) 
Y¥=—X+2, 


(See Fig. 4.) If another form of the solution were available, we should get in 
this way V(x,y; X,Y) when Y—y lies between +(X—-x). The reasons for 
this are (a) the solution of the Problem of Cauchy is unique, and (b) the Riemann- 
Green function does not depend on the nature of the curve carrying the Cauchy 
data. 
This is precisely what R1EMANN did in the special problem which interested 

him. The equation with which he was concerned was 

277 277 OTT 

ay re ea (6.6 


* Note that, if Y<y , the upper limit in this integral is less than the lower limit. 
Similarly in (3.5) if X< 4% . This is of importance in what follows. 
Arch. Rational Mech, Anal., Vol. 1 23 
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in our notation, where « is a constant. He solved the problem of Cauchy under 


the conditions = 
U=0, 3 =G() 


when x= 2% )>0 by means of a Fourier cosine transform with respect to y, and, 
by a comparison of the two forms of the solution, he obtained, in modern notation, 


le.) 
goth 


V(x, y; X, Y) = ———__; [ cos [A(Y — y)] x 


cos a7 xXet : 
(3.7) 
x ey (4 2) Iya X) — Jung (AX) yal 2] aa 
when Y — y lies between + (X — x). He then remarked that each Bessel function 
can be replaced by a definite integral, so that V is then expressed as a triple 


integral, which can be reduced to a hypergeometric function. He merely quoted 


the final result, which nowadays we should write as 
We r ges ed . ie 
V(x, y; X,Y) =() Patt 8) (3.8) 
where 
Fete Sad toed ee 9 
: 24 xX : 


and contented himself with verifying the result. WEBER’s editorial remarks 
in the ‘‘Collected Works”’ shed no light on how RIEMANN reached this remarkable 
result which antedated by twenty years the evaluation of similar integrals by 
SONINE®. 

RIEMANN was able to obtain a second solution of his equation by means of 
a Fourier cosine transform because the variables are separable. We show in the 
next section that the method can be extended to any equation with separable 
variables. 


§ 4. The generalisation of Riemann’s method 
When the variables are separable, equation (2.1) can be written in the form 


2U aU eu oU 


ay? | y 


where a, p are functions of x alone, b, g are functions of y alone. We consider, 
as is usual in such a case, the pair of equations 


ao dd 
aa +247 ++ A40=0 (4.2) 
1? p dy 
Be tog ++ 2) p=0 (4.3) 


where A” is the separation constant. Suppose that 0, (x, A), 0, (x, A) are linearly 
independent solutions of (4.2) whose Wronskian 


* SONINE: Math. Ann. 16, 46 (1880). See also Watson, G. N.: Theory of Bessel 
Functions, pp. 411—412. Cambridge 1922 
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we denote by W’(x, A). Similarly y,(y, 4), v2 (y, 4) are linearly independent solu- 
tions of (4.3) whose Wronskian is W’’(y, 4). 


We now try to construct a solution of the problem of Cauchy for equation (4.1) 
of the form 


= J {f(A) p1(y, 4) + fo(A) aly, A)} Or (x, a) dd 


where integration is over a fixed range which depends on the form of the equation. 
We wish to find U(X,Y), given that U=0, 0U/dy=F/(x) on y=yo, where yy 
is a fixed but arbitrary constant. This means that we have to find /,(A), fo(A), 
if possible, so that 


SAA (A) G1 (¥o, 4) + f(A) G2 (vo, A)} A(x, A) da = 
SAA pe cee ae 
where @;(y, A) denotes ae lf 
= f f(A) 9,(x, a) da, (4.4) 
these conditions are satisfied = 
f(A) Gi (Yo. 4) + fe(A) G2(¥o, 4) = 0 


f(A) Pi (Vo, 4) + f2(A) Po(¥o, A) = F(A), 
that is, by 
_— __ F(A) P2(Vos A) _ F(A) (do, A) 
h (A) me W’' (yo, 2) ) he (A) : ww’ (Vos 2) . 


Hence we get 


U(X, Y) = =[f0 ) {91 (0, 4) P2(Y, A) — (Y, A) G2(vo. A & 7 


But if the solution of the integral equation (4.4) is 
= J ,(x, 2) F(x) dx, (4.5) 


we then have 


U(X, Y) =[ fie ae z AG. 1) £6, (yy, A) P2(¥. A) — P1(Y, A) polo, A} dx da. 
Comparing this with 
X+Y—y% 
U(X,Y)=3 Le ree ) V(x, Yo; X, Y) dx, 
A 


we see that 


V(x, 9; X,¥)= £2 f ATONE fo (9, 2) gal VA —prl(¥, A) poly, A} ad (4.6) 


provided that x lies between X+-(Y — y), but that, when x is outside this range, 
the value of the integral is zero; the upper sign is taken if Y>vy, the lower if 
Ves 

Interchanging the roles of x and y, we see also that 


Vix,v;X,Y=+2 f BOR? EI 69, (x , 2) —0,(X, 2) Io(x, A} aa (4.7) 


23* 
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if the solution of the integral equation 
J g(A) ily, a) = G(y) (4.8) 
z (A) =SA(y, 4) Gl) dy, (4.9) 


provided that y lies between Y+ (X — x), but that, when y is outside this range, 
the values of the integral is zero; the upper sign is taken if X > x, the lower if 
NORGE 

Lastly, we observe that two other formulae for V can be obtained by inter- 
changing the suffixes 1 and 2 and changing the sign of W’ (or W”’). This comes to 
replacing #, and #, in (4.6) by J, and the corresponding #,; and similarly in (4.7). 

All this is entirely formal. To give a rigorous discussion would involve a 
careful consideration of the various cases which can arise, and the essential 
idea would be lost in a mass of detail. The important point is that, in many 
special problems, this technique gives the form of the Riemann-Green function; 
having got the form, one can usually justify the result quite easily. 


§ 5. Some particular examples 


As a first example, we consider the equation of damped waves 


aU e2U = 
Ox Oy? pli (5.1) 


The two associated equations are 


ao 
Sz t(2+1)8=0, 


Ce +o =0 
Hence 
’, = en ie 1. 3B, ies eit yard Ww = 21 (2 +1) 
and 
Opa eqe ee, Gee et™ We Oe A 


If we use formula (4.7), we are in fact solving the equation by means of 
complex Fourier transforms in y. For the range is — o<y<oo, and @G, is 
defined by the fact that the solution of 


fede?" da = Gly) 
is ae 
g(a) = J G(y) @(y, a) dy 
Hence 
Ply, A) re eb 


It follows from our formal analysis that 


CO 


4 OF ee oe SEI aoe x See 
bs | GOT ee ae 


is V(A? 41) 
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provided that y lies between Y+(X— x), and that otherwise the value of the 
integral is zero. We may evidently suppose X > x and take the upper sign; and 
since the integral can be written as 


co 


1 Z sin {VA? +1 (X — x)} 
P| cosy 


—_ —_— ? 


V2 +41) 


we may also suppose Y>y. To evaluate the integral when y lies between 
Y-+ (X — x), we put 


A— “42 Reoshasy Y= 7=—Rsinhe, <= sinh? 


and get 
1 i i R cosh (9—a —i Rcosh (+e 1 
Vaz | oR — cirri ta ag — © {HS(R) + HER} = FalR) 
Or 


V(x, y; X,Y) = Jol V(X — x)? — (Y— y)}]. (5.2) 


When Y —y does not lie between + (X — x), we may evidently consider only 
the case Y — y > X — x: if we make the substitutions 


X—*=RKsinha;, Y=y=Reosha, A=sinh? 


we find that the integral vanishes. This does not mean that the Riemann-Green 
function vanishes in this case, but only that to find it we must have recourse 
to formula (4.6). 

When X — x lies between +(Y— vy), we use A?—1 as separation constant 


instead of 42, and use Fourier transforms in x. The result is that in this case 


fo) 


n A eee sin {| Moe Ces AA V(x, y; x Y) 
Con, 


zs \(=1) 


where the upper or lower sign is taken as Y> or <y, and that otherwise the 
integral is zero. It turns out that 


V(x, y; X,Y) =1,(V(Y—y)? — (KX —)3]. (5.3) 


The equation (5.3) appears, very slightly disguised, in RIEMANN’s paper as the 
limiting form of the Riemann-Green function for the equation 


SU a Aka 
Ox? “Oy | <a Ox 


10 (5.4) 
where « is a constant *. 

Since (5.4) has a singular line x =0, we restrict consideration to one of the 
half-planes, say x 0, into which the plane is divided by the singular line: in 
what follows, x and X are positive. The associated equations are 


@9 | 20 dd 


PA RO 
Ga ee ee 
ap : 
ae 


* Put x=a+t¢ and make « tend to infinity. Then put U=e’U,. 
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74 


so that we may take 


Dane © I, ahve), Da ben” Joe aes 


Ww — 2608 7% 
Vary Leper 


FA 4 AY te oe 
Gie=e sae Yo @ W == AG 


if « —4 is not an integer, as we suppose to be the case. If « — 3 is an integer, 


we must take 


jad" T_{0), =2%y 40), Was 


yee? 

but the final result will be the same. If we use (4.7), which means that we are 
solving (5.4) by complex Fourier transforms in y, then y and 4 vary from — oo 
to + co and 9, =e'?”/2m. We then have, when y—Y lies between +(%—X), 


Doda sae pf AO Te-4A) -2 BX) — fry AX) Jy_a(A x) da 
360s aa Xo & ‘ 
By (5.5) 
4oyits 


= 25 [ 008A —¥) 49 aA X—L-g(AX) Iya} a 
cos ma X*~# 2 
0 
since the expression in brackets is an even function of 4; for other values of 
y—Y, the integral must be zero. This is the expression given by RIEMANN. 
The plus or minus sign is taken according as X> or <x. 


To evaluate Mile integral, we use the result* that, if a=b>0, c>0, 


Cc 


fooseaR(ba) I, (ai di= = @, (S55) (0<c<a—b) 


0 


=| (a—b<c<a+b) (5.6) 


=0 (c>a+b). 


We may evidently suppose that y>Y, x< X; but we must also have regard to 
the position of (x,y) in relation to the characteristics through (X,Y) and its 
image (— X,Y) in the singular line. These characteristics divide the half-plane 
into six regions as shown in Fig. 5. 


In III, y—Y >x-+X, and so the integral vanishes. In the part of JI where 
x<X, we have X—x<y—Y<X-4-4 and therefore 


+h om 
conan xed | SAY) 1(A) Fh a(AX) — Sea AX) Iyaldad} ad 
0 
eee ee ae pa P (ea eae 


2cos mwa X% \ ee 2% X 


fy pekinese only reference I can give is to OBERHETTINGER, F.: Tabellen zur Fourier 
Transformation, p. 65. Berlin: Springer 1957. That the integral vanishes when 
c>a-+b was proved by Battery: Proc. London Math. Soc. (2) 40, 37—48 (1936). 
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Lastly, in the part of J where y—Y>0, we have 0<x<X, 0<y—Y<X—x 
and so 


= anion Bee ay za) w+ X*— (y—Y)? 
ear gh) = mee <2 {Ox G Pe Cn )=O-« Dax I 


or 


; on x? —Y . 
Via yi XY) = 3 P(e tea), (5.7) 


It is now evident that the integral (5.5) vanishes except in the regions J and I’, 
and that, in these regions, the Riemann-Green function is given by (5.7). 


To get the Riemann-Green function in the rest of the half-plane, we use the 
formula (4.6). This means that we are solving the partial differential equation 
by means of a Hankel transform in x. The 
function #@,(x, 4) is to be chosen so that the 
solution of 


He, 
XX 
s 
It follows that, when « >0 
Oa Anes fe (Ada 


Hence we have the result that 


a2 ek a) nS 3 (A x) AEP (A X) sin A(Y std y) a). Fig. 5. Regions for calculation of V 


for Eq. (5.4) 


is the required Riemann-Green function when X — x lies between +(Y — y), the 
upper or lower sign being taken according as Y > or <y, and that otherwise the 
integral is zero. 

To evaluate this integral, we use a result due to H. M. Macponatp%, that if 
== 0s — 


Ce 


amici (bd) J, (aa)dA=0 (0<c<a—b) 
- tp ete ye : % 
 2V (ab) eA 2ab -} (4 ee a ) (5 ) 
COS V 7% c2 — a2 — b? 
a x V(ab) ( 2ab | (Gate 2) 


As before, we have to consider the six regions into which the half plane is divided 
by the characteristics through (X,Y) and (—X,Y); and we may evidently 
suppose that Y — y is positive. 


* Cf. Watson: Theory of Bessel Functions, pp. 411—412. Cambridge 1922. 
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In the part of J for which Y>y, we have O<Y —y<X —x; hence by the 
above formula with c= Y—y, b=x, a=X, the integral vanishes. Similarly in 
the part of I’ for which Y>y, we have 0< Y—y<x—X, and the integral 
again vanishes by taking c= Y—y, b=X, a=x. Next in the part of IJ’ for 
which x< X, we have X —x< Y—y<X-+4, so that there 


V(x, 9; X,Y) = 2, Pf (5.9) 


x? + X2 a (y _- ye 
es 


De 


since P_,—P,_,..And in the part of JI’ for which +> X, we have x—X = Vos 
y<x-+X, so that. the same formula holds. Lastly in 11’, Y—y>X-—-4)se 
that in this region 


AC oP. G be Z “di feces efioee zs 


sina — QO 
: a: 5 o— Da 


(5.10) 

7 
To sum up, the Riemann-Green function for equation (5.4) is given by (5.6) 
everywhere in the half-plane x >0, except the regions JJ and J/I’ where formula 
(5.10) holds. The unusual behaviour of V in the regions JJJ and III’ does not 
appear to have been noted before: yet it would be of importance if one tried to 
use RIEMANN’S method with data on the axes of coordinates and wished to find 
the solution when Y>X. If « is an integer, V is identically zero in [JI and JII’. © 


As a last example, let us consider the equation 


@U | 2a aU aU , 2p eU 


ae Seema oS ete 
One eb ay? | y ay 


(5.41) 


where «, 6 are constants. The associated equations are 


a9 2a dd Shaye 
dx x ax AU = 0 
ay 2B dy a Mee 
dy? ' y dy ep =0, 


so that we may take 


Oj a" fe (AR), ye OA 
5 2 Ti 
i- ie yer 2 

A=V? hb sy), ge— oP holy), We = 20828 | 
As there is symmetry in (x, y), we may use (4.6) or (4.7). If we use (4.6), we can 
deduce the formula corresponding to (4.7) by interchanging x and y and « and p. 

We take x, y, X, Y to be positive, which corresponds to solving the partial 
differential equation in the first quadrant of the (x, y) plane. As in the previous 
example, 


By(x, 2) =A y(ax) 
and so, if 6 — 3 is not an integer, the integral we have to consider is 


7m yttd yb tS 


BE cos x eho [AL sQ9 LAxx 
0 


x {Jp4(49) Jy_p (AY) — py (AY) Jy_pl(Ay)} aa 
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which should give V(x, y; X, Y) when X — x lies between +(Y—y), but should 
be zero otherwise; the upper or lower sign is taken according as Y> or <y. 
The final result will, however, hold if 6 — 4 is an integer. 


To evaluate this integral, we use the following result in the theory of Fourier 
transforms, that if 


ie, @) 


F(u) =f f(a) cosAuda 


G(u) =f g(A)sinduda, 


then 


the last integral being a Stieltjes integral, since G(w) in our problem turns out 
to be discontinuous. 


If we take 
f(A) = + {p44 9) Jy (AY) — yp (Ay) p-a(AY)}, 


where the upper or lower sign is taken according as Y> or <y, it follows from 
(5.6) that 
mi aetna y2+ Y2— w2 : . rae 
Fu) = |/ oy cos 7B Pa rer if 0<u<|Y— yl, 
=i) if u>|Y—y]|. 
If we take 


(A) = Sx-4(A2) ey (AA), 
it follows from (5.8) that if «>0 


G (u) =(0 Hii (0) Se 
1 : wet X2 — 2 
=p | 
] (20 x X) if Digg XK 


X— x| 


if |X—x«|<u<X+x 


sin a 7 u? — x2 — X? : 
Wana © if Se) i ee eee. 


where P,_, =P_,. It follows that 


oo |¥—y| 
fAfAg(Ada= J F'(u) dG (uv). 


We have to consider separately the cases X>Y and X< Y. In Figs. 6 and 7, 
the quadrant x20, y=0 is divided up into regions by characteristics 7 + y= 
constant. In the regions J and I’, |Y— y|<|X — x| for both cases; G (wu) vanishes 
everywhere in the range of integration, and so the integral (5.12) vanishes there 
as our argument led us to expect. This, of course, does not mean that V is identi- 
cally zero there, but only that our argument does not determine V. 
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Now let X¥>Y, so that the situation is as in Fig. 6. Then in II’ we have 
—(Y—y)<X—*<Y—y. Hence 


oo ; Y=—y 
f Arla) g(a) da=[Fm@) Geist + fF 
0 |X —2|+0 
Viens! é 
= cos af lp ye Y2— P etal apo a? +X? ve? 
~~ aVaxyV) -+| oe ee) a ( 7 ( 2xX 
m/2 
= 228 [Pali tm) + f Ppt + ncostd) dP_q(t + ésin®9) 
nV (aX y Y) 
where y' 
é (Xe 4 (Va ye 
; PEG IK. 3 
tl 6a) org tC Sa 
1] 2 Y ) 


(OX4¥) 


(0,X+V) 


(0Y-X) 


Sy 
Ss 


ORT SAY) 7¥*V,0) 7V-X,0) (X+V,0) & 


Fig. 6. Regions for calculation of V for Eq. (5.11) (X>Y) Fig. 7. Regions for calculation of V for Eq. (5.11) (Y>X) 
this follows by making the substitution 
u* = (X — x)? cos? + (Y— y)?sin?d. 


We thus have the result that, if X>Y, 


eae) 7/2 
V(x,y; X,Y) = ne P_s(1+m)+ [ P_p(1+n cos? d) dP_, (1+ &sin?9)| (5.13) 
‘s 0 
in II’. Since € and 7 vanish on the characteristics through (X,Y), this function 
obviously satisfies the prescribed boundary conditions. When I discovered this 
result in the summer of 1957, I thought it was new; but I have since then learned 
that CHAUNDY* gave it nearly twenty years ago in a paper which does not appear 
to be well known; his derivation is completely different from mine. (See § 8.) 
Similarly it can be shown that the formula (5.13) holds in JJ. 


* Cuaunpy, T. W.: Quarterly Journ. of Math., Oxford Ser. 2, 234—240 (1938). 
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The situation is rather more difficult in region III, for there we have 


[eee ee eee 
so that 
2 : ; X+x*—0 
Vey X,Y) = Fw COMA +S Flu) dG lu) + 
Xa-x| 20 
LNG 
+ [F(u) G(u)jxfi%9 + Sf Flu) dG(u). 


X+x+0 


The difficulty is that G(u) tends to infinite limits as w+ X + x +0, which makes 
the resulting formula more complicated. 

If, however, Y>X, as in Fig. 7, everything goes as before: the Riemann- 
Green function is given by (5.13) in JJ and JJ’, and more complicated formulae 
hold in IJ and III’. ‘ 

Lastly, to get V inJ and J’, we interchange the roles played by x and y, and 
by « and #, and make the corresponding changes in Figs. 6 and 7. 

A different form of the Riemann-Green function for (5.11) has recently been 
given by P. HENRICI*, w72. 


X \% 4) \P o- “: hie aaa? Viet ; 
eal lee sth a al re ele ae reel (5.14) 


Ri = (x —X)*— (y—¥)* 


where 


and F’, denotes APPELL’s hypergeometric function **. 
As a last example, I would like to mention briefly the equation 


GU I OU oh au . 
——— — }? — 3 al 
for which HeEnrici (loc. cit.) has given the Riemann-Green function 
8 2 
Vm R 4 
Fi aes 5 BY pee ieee » R2 ‘ 
V(x 95 X,Y) =e (6.4 rena rae Ie (5.16) 


This equation has its variables separated and could be discussed by the methods 
of this section. Another method would be to apply the limiting process suggested 
by RrEMANN*** to equation (5.11). If we put *=(«+x1)/x and make «+ oo, 


(5.11) becomes 
oo) 
or? 


gU  eU | 2p aU 


= “AL, 
ot ay? " y by’ (5.17) 


2% 
and the expression (5.13) for the Riemann-Green function then reduces to 


B 
Vii, y; i, X) ary ee) x aX 
7/2 A ee Tho (5 18) 
ee (1 +7’) + | JB a(t +. cos? P) d Jy (x sin ? Viv y)? — (T— i})| 
0 
* Henrici, P.; Z.A.M.P. 8, 169—203 (1957); in particular Table 2 on p. 180. 

xx See WHITTAKER & Watson: Modern Analysis, p. 300, Ex. 22. Cambridge 1920. 

*xkk See the footnote on p. 331. 
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when 7 —t lies between +(Y —y), and 


by A) ee 
/ DVey. ; 
It readily follows that the Riemann-Green function for 


2B eu 
y oy 


02 U 
o2U We ieee 
ot? Oy? 


is obtained by omitting the exponential term in (5.18). 


§ 6. Connexion with the work of Hadamard and Marcel Riesz 


Riesz’s theory of integrals of fractional order was originally devised to give 
a solution of the wave equation in any number of dimensions and was subse- 
quently extended by him to deal with the equation 4,U =—0 in a Riemannian 
space of normal hyperbolic type. The extension to a non-self-adjoint equation, 
the Euler-Poisson-Darboux equation, was worked out in detail by RutuH M. 
Davis*. 


In the case of two independent variables, say for the equation 


aU GAO! oU oU = 
L = ap a = Sa Xx - BACs i 
U AF Bye + a(x, y) a b(x, y) By + c(x,y)U OR 


this work depends on finding a function V%(x, y), depending on a parameter «, 
such that 


L* yet? — ya, 


(L* being the operator adjoint to L) and also such that, if « is large enough, 
V* vanishes on the characteristics through a general point (X, Y). This function 
V* is closely related to HADAMARD’s elementary solution. It has the form** 


where 


and V, is independent of «. 
In general, V** has a simple zero at «=0. The equation 


L* y2ar2 22 ye 


then implies that L* V2 = 0; and V2 turns out to be the required Riemann-Green 
function. If all the V; were zero, V?* would have a double zero at «=0, and in 
this case [0V**/éa],9 would also be a solution, actually a constant multiple 


of HApDAMARD’s elementary solution. This, however, only happens if a, b,c 
are all identically zero. 


x Davis, RutH M.: Annali di Mat. (IV) 42, 205—226 (1956). 
** It is simpler to consider V2% rather than V%: it avoids annoying fractions 4% in 
the denominator. 
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It is, however, possible to proceed in a slightly different way, analogous to 
the method of FRoBENtus for ordinary linear equations. Let us consider, for 
example, the function 

peti-t 


wee — BE Eger 9G 
4 HP (a+9) Pa +3) 
It is readily shown that 
G2W2e+2 ON oa Dyes Agere 
Ox? ay? DP(a+1) F(a +1) © 


Since the right-hand side has a double zero when «=O, the functions W? and 
[oW?**+?/dq],-9 are Solutions of 


aw ew 


Ee = ay? W = O. 


The former is the Riemann-Green function I (V0), the latter an elementary 
solution of the form 


Bie Pe eee 
I, (VL) log I 2 ee 


w(t) being the logarithmic derivative of /’(f). We note that the Riemann-Green 
function is the coefficient of log /°in the elementary solution. 


The difference between the procedure we suggest here and that used by Miss 
Davis, is that we do not require V; to be independent of «. We consider instead 
a function 

ee) peti} we 
= ip 


wee — Fe A | ne 
4 P(a+j) P(a+9) 27 
It easily follows that 
gy ilar Valk OF at ME 
pe Gay Lior} 
provided that 


2a+2 re 
ie ey vege 0 
GV 222 ‘ Wweute - 0% | We b| Veo ( ue ) L* Wat? — 9 
€ ds tai (x +4) 7 [( X)a l (Vv ) | 7 a7 ja 
where 
Wea ee eee A 
ae Crees Be ee oy 


If we require W?*t? to have the value unity at (X,Y), we find that 


(7, 9) 
Wet Sex Hf {(x —X)a+(y—Y) od} 


(X, Y) 


a 


S$ 


integration being along the straight line from (X,Y) to (x, y). Hence if (X, Y) 

is not a singularity of either coefficient a or b, Wo**? is regular in a neighbour- 

hood of (X, Y) and is independent of «. We denote it simply by W. Then since 
402 P*-1W, 


* T772a+2 __ tS 
pe T(e+1)0(«+1)’ 
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it follows that W2 and [éW?**?/da],-9 are solutions of L*W=0. The former 
will be the Riemann-Green function, the latter the elementary solution; the 
Riemann-Green function is the coefficient of log J’ in the elementary solution. 

We are, however, somewhat anticipating the result. We must first show that 
the remaining coefficients W?**? can be determined uniquely, that the resulting 
series for W?%'? is convergent and can be differentiated term by term. 


Since, 
| ree pa 
a(x X) ] b(y Y) A ds’ 

the equation for W?**? can be written as 

d 724+ 2 7 L* War 

Se Cam] ir | 4 (a +7) s%tI eae _=0, 
and so 

’ | 
ache e ee tag, (a +7) W% CASUT *Y72a+2 7. 
Me Sa eae i eee 
(CY) 


where A; is a constant of integration. If this is to be regular near (X,Y) no 
matter how large « is, A, is zero for all 7. Hence 


(x, ¥) é 
W=exp f {(x—Xat(x-V)}*, 
(X, Y) “ 
| ie (x, ¥) : 
2a+2 a + J) Pee 2a+2 
(X, Y) 


The convergence proof, when the coefficients are holomorphic, goes through 
as in HADAMARD’s book, and the rest follows. For 


Ww=W+ y fc “Ss 
Pan ivie 
gives 
ew? am 1 
ag ap te XW + OW), 
ow? am 1 : 
ay ay relma Lite kCAC IE 


Hence on the characteristic y—Y = x — X, we have W? = W, and 


ow? ow? =aW, ow, 
a. Bye ORE Bye ee oa 
and similarly on y—Y = — (x — X), 
ow? ow? ‘ 
i a = (a — b) W?. 


Thus W? is the required Riemann-Green function. In theory, we can find the 
successive coefficients W,**? by a simple recurrence process, and then obtain 
the Riemann-Green function and the elementary solution at the same time. 
In practice, this may prove difficult. 
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For instance, for the equation 


PUL OUTS) oS) ay oe 
ene ee ee 
it is easy to show that 
a ape v(4 —v) 7/ vee : a 
Wat, weet 29) F(2, 0445 0+ 2; =e 


but I cannot find any simple formula for the higher coefficients for general 
values of x. In the particular case «=0, we have 


we = (— 1 Ot DL» +3) 
(vr) P(A — y) (7X)! 
and hence 


eee y 6 Bes tage Dien. Mi IAS NY 2 
WwW = Flv, 4{— v5; 1; x) Bel 


§ 7. The integral equation method 


If we use characteristic variables (v, s), the equation takes the form 


_ @U aU eU a 
Be egos Pe ie a 
where A, B, C are functions of (7,s) alone. The Riemann-Green function 
V(r, s; R, S) satisfies the adjoint equation 

jc 8 (et ale as (a ye 


é 
ain a Vv = 
AE ates (BY) +-€V,=0 


Os 
under the conditions 


A Ay ON WAST: he 


SY ON Sea) 
os OY 


and 
Vers aes eae 


If V denotes V(7,s; R, S), we then have 


Jar fast, ai Wi abres Wen Wae tegraeti) 


Hence 
V (t,u;R, S)—V(t, S;R,S)—V(R,u; R, S)+V(R, S;R,S)—f Alé,s) V(i,s;R,S)ds+ 
S 
u t 
+ f A(R,s) V(R,s;R,S)ds —f Bir, u) V(r, u; R,S)dr+ 
S R 
t t uU 
6 EBS) V(7595 R, S)dr + fdr fds Cr, SAVAGES lees he 
R R S 


Now 
Vik aS Sie 43 
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Hence, by the condition on the characteristic s = S, 
t t 
fBir,S) V(r, S;R,S)dr=J V7, 5; 8, S) dr = VES; R,S)—1, 
R R 


and similarly 
[ARRS) VRS BiG Sld SV (CRS) —— 42 
s 


It follows that 


t 
Vie, ks, oA + [ Bir, u) V(r, u; R, S) dr + 
R 


uU t Uu 
tfA(ts) VE seR,S)as =~ (ariiahs C(tos\ 1 eo 
Ss 


R Ss 


Conversely, it can be shown that, under certain continuity conditions, this 
integral equation has a unique solution, which is the required Riemann-Green 
function. 


Unfortunately this is not a very good way of actually finding V. Let us con- 


sider two examples. 
For the equation 
aU 


TE eae oa 


where A is a constant, the integral equation reduces to 
r s 
Vir,s;R,S)=1+Afdtf[duVit,u; R,S) 
Ry Ss 


with a trivial change of notation. If we try to solve this by successive substitutions 
or (what comes to the same thing) by putting 


Vir, s: R&S) ==> AAG 5) 
0 
we get 


fo(7,s) =1 


is as s) = dt J du ae (é, u) . 
It follows by induction that 
fal. s) = HAN SV 


and hence that 


Vins; RS) Eee eon 
But if we apply method to 


@U  a(t—a) 
pa A a J= 

ov Os (y+ s)? v 
where «% is a constant, we have 


Vir,s; R,S)=1—a(4 — a) fat [au Al 
Ro US, 


t,u; R, S) 
(fxr 4)? 
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The known result is 


VQ,s:R,S)=P_.j1+2 ee, 


but solution by successive substitution gives 


V = Dial —1)"7,(7, 8) 


where 


AS 4 


f, (7%, S) = fai fentieses 
(¢+ u)? 


(R+s)(r+S) 
2, N10) =e 

A Ne Ge HRS) 

This is, in fact, the correct first term in a rather unorthodox expansion of the 
Legendre function. It does not seem possible to get any further in this way. 
HENRICI (loc. cit.) attributes the integral equation of this section to I. N. VEKUA. 


In particular, 


§ 8. Chaundy’s papers 
In a series of papers written in the five years before the war, CHAUNDY* 
discussed in great detail what he called partial differential equations of ge- 
neralised hypergeometric type, namely equations for the form 


IT f,(6,) L= Ny XQ.-- %n TT 8, (0) Z 


where /,, g, are polynomials and 06, denotes the operator x,0/éx,; in particular, 
he found the Riemann-Green function for equations of this type when = 2 
and the order of the equation is 2. It would take too long to give a complete 
account of this work. Suffice to say that it involved a very skilful use of symbolic 
operators and great insight. It should be noted that in CHAUNDY’s work the 
term involving the highest derivative may be 


in our notation; this necessitates a slightly more general form of the adjugate 
equation and an alteration in the definition of the value of the Riemann-Green 
function at the vertex. It is, however, quite easy to transform his results into 
our notation. This difficulty does not arise in the example we now give. 


In the 1938 paper, CHAUNDY considered the self-adjoint equation ** 


au fm (m + 1) Mz (M +1) , Ms(mz3t+ 1) — mg(Mg+ Dh 
évés \ (r+s)? (y—s)? ' (1—~,7s)? (1+7s)2 


* CHAUNDY, T. W.: Quarterly Journ. of Math., Oxford Ser. 6, 288— 303 (1935); 7 
306—315 (1936); 8, 280—302 (1937); 9, 234—240 (1938); 10, 219-240 (1939); 11, 
101—110 (1940). 

xk There is a misprint in the enunciation at the beginning of this paper. The sign 
of the term involving mg is wrong. 
Arch, Rational Mech. Anal., Vol. 4 24 
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where m,, ij, ™g, ™, are constants. The Riemann-Green function satisfies 
this equation and is equal to unity on each of the characteristics y= R, s=S. 
He now makes the change of variables 


_ _ W—R)6~S) ey 
oS yee aa) a eer ery 
_ _ _—R)(s—S) v-  Y=Be=—S) | 
1 ate ee ee ee 


and shows that the Riemann-Green function satisfies the simultaneous equations 


6, (6, + 63 + 63+ 54) v = x, (6, — m,) (6, + m, + 1) v Cpe BORN ee 


and is equal to unity when x, = *,=%,=%,=0. Such a solution is 


oe) 4d , re 1 
ae DOA TE ~m) Caner) Gee een 


11, %0,%3,%4=0 \s=1 


where («),=a(«#+1)...(a+7—1). This he denotes by 


My Mz Mz Mg 
%y %_ %z Xq 


aK 


If any m, is zero, the corresponding terms in the multiple power series are re- 
placed by unity. 


In particular, if m,=m, m,=m,=m,=0, the Riemann-Green function for 


eu  m(m+1) 
dvds (r +s)? 


U 


is 
= A(T )= pa ae ts 4 = F(— m,m + 13 1; 4) = (1 — 2%) 
or 
ae (y — R)(s— S) 
v Belt 2 Cee 


which is equivalent to (5.9) above. 


Again, if m=m, m,=n, mz=m,=0, the Riemann-Green function for 


Ou = {ee pa} 
Ov Os (y +s)? (7s)? 
1S 


and this CHAUNDY recognises to be 


1 
v=F,(1— 2%) —2x.f P,(1—2%,+2%,1) P’ (1 — 2x, t) dt 
0 


which is, in fact, equivalent to (5.13). 
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§ 9. The use of contour integrals 


In his 1936 paper already cited, CHauNpDy obtained an expression for the 
Riemann-Green function of the equation 


Cu ou , OU 


ces) ov OS "Gy i as 


=0 


where m and n are positive integers; he used his method of factorising differential 
operators and found the Riemann-Green function as a repeated derivative of 
a rational function. He then observed that this derivative of a rational function 
could be expressed as a contour integral which had a meaning even when m 
and 7 were not positive integers. 

Twenty years later, A. G. MACKIE, 
quite unaware of CHAUNDY’s work, 


obtained contour integral formulae for 
the Riemann-Green function of 


¢’ 
au 1 Ou Ou ; 
ar | rare (2 or +B ae =(0 (9.1) Fig. 8. Contour for Eq. (9.2) 


where «, § are constants. He was led to this when extending his work* on the 


solution of certain gas dynamical problems by complex integrals: so far, MACKIE 
has not published this result. 


The Riemann-Green function associated with the point (R,S) satisfies the 
adjoint equation gua function of (7, s) and the original equation gua function 
of (R, S). It is easily verified that, if zis any complex parameter, 


(zg — R)~* (zg +S) 
is a particular solution of the original equation and that 
(+s) (273 (z+ 


is a particular solution of the adjoint equation. We are thus led to consider a 
function 


Ne ae esl 
v(r,8;R,S)= fa { i _ 2 * = f(z) dz. (9.2) 
Cc 


The integrand has branch points at 7, R, —s, -- S; and we can make the integrand 
one-valued by cutting the z plane along the real axis from —s to — S and from 
+yto +R. We choose the branch that is real for big enough real values of z. 


If v is the Riemann-Green function, it must satisfy the conditions 
Ga (KR s)he +S) =” son: 7 = K 
woe + S/P(R+AS)F on s=S. 


We take C to be the contour shown in Fig. 8, and see whether we can choose 
}(z) so as to satisfy these conditions. 


* Mackie, A. G.: Proc. Camb. Phil. Soc. 50, 131—138 (1954). — Journ. Rat. Mech. 
Anal. 4, 733—750 (1955). — Proc. Roy. Soc. A 236, 265—277 (1956). 
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Now when 7=R, we have 
R35 f (s/s 2) dz 
ant ey) gas) 2 


O(Ry.5; hp) 
where C, is now a simple closed contour surrounding R; the singularities —s, —S 
lie outside Cy. If f(z) is regular inside and on Cy, we have 


(Ry Ry) =f (R) 
which is what we need if f(Z) is unity. Thus 


Nak oo. 
v(,s;R,S) =e | te) Cari) dz (9.3) 


24 


satisfies the conditions on y= R. 
To deal with the condition on s =, we observe that 


<1 ESR, ap re 
v(r,s; BR, S) al GR es) ai 
2 


where C’ is as shown in Fig. 8; this follows from the fact that 


fa tee 


as 9-> oo. If we now put s=S, we find that the second characteristic condition 
is also satisfied. Hence the Riemann-Green function is, in fact, given by equa- 
tion (9.3). 


If we make the substitution 


equation (9.3) becomes 


a vts*@+spP% 1 fe ee is 
ICAP Su tey ea) (etsy axis! (= 4) eas 


where 
xy — (R= (S=s) 
(R + S) (7+) 


and /’is a simple closed contour containing the points 0 and 1 but not —1/X, 
From this follows that * 


v(7,s;R,5)= vs ee Cay vines ls 
(R +)? 
As CHAuNDy pointed out, the method is capable of generalisation, and in his 
1939 paper he obtained a contour integral representation of the Riemann-Green 
function for 


(R —1) eos) 


(R+S)v45) OA) 


@U | 2u+1 0U #U , 2v+1 aU 


‘ i ey | 
Ox? x Ox Oy? y oy 


* This can be expressed in a symmetrical form by using the formula connecting 
F(a, b; c; «) with the hypergeometric functions of argument 1— »%. 
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| where yw, v are constants, viz. 


(0+) 
1 22 NUS TENG ye 7 z 4 
2; () (+) | K, (= x1) K,(—yt) K,(X 1) K,(Yd tdt. 


| cI 


In the first instance, «— } and y—# are integers; but the result is extended to 
any mw and ». 


i § 10. Titchmarsh’s derivation of » (r,s; R, S) for the equation 


of damped waves * 
If we use the characteristic variables 7, s, the problem is to find the solution of 
Pu 
Sey) 
Ov OS 


which has the value unity on the characteristics y= R and s=S. If we put 
y—R=x,s—S=y, we have to find the solution of 


oie tae 
ax Oy 
which is equal to unity on the axes of coordinates. 
Let us write ee 
VC, y) Ss [° (x, y) e°* dx 
V(2) ‘ 
where € =£-+-2y and n>c>0. Then 
(oe) (ee) 
a ae fo.0) 2 
a papas oes dv ett dx = Ale ov =| oa =| - GO Ge dx 
oy = (an) J oY V(2a) ley 0 iC\(27) » dx OY 
0 
But dv/dy vanishes when x =0. Hence 
Cc co 
2 a oe 
ve Jif ev ft dx = —~—*__ [ uel*dx 
oy icV(22) « Ox OY iV (22) « 
so that 
ave Vv. 
Cyne TG 
This gives 


Now make y tend to zero. Then 


[e,0) in CO F ; 
A =f» G0) dy Can — 
It follows that 
tk+oo at 
sina— gh, f ember tetl — Kone 
ik—oo 


and hence that 


u(r, s; R, S) =1,(2V f(r —R) (s—S)}). 


* TircHMARSH, E. C.: Theory of Fourier Integrals, pp. 297—298. Oxford 1937. 
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Although the equation of damped waves is very special, this work of TITcH- 
MARSH’S is of interest in that it is a direct solution of the characteristic boundary 
value problem satisfied by the Riemann-Green function, and is the only one of 
which I am aware. 


§ 11. Conclusion 


In this report I have, I believe, covered all the known methods of finding 
the Riemann-Green function and have listed all the known cases, apart from 
trivial changes of the dependent or independent variables. I should like to take 
this opportunity of thanking Harvard University and the staff of the Harvard 
Department of Mathematics for their kind hospitality during the Fall Semester 
of 1957—58, which gave me the leisure to complete this work, and also Dr. A. G. 
Mackie of the University of St. Andrews for the many interesting discussions 
we had concerning the contour integral method of § 9. 
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On a Singular Boundary Value Problem 
for an Equation of Hyperbolic Type 


E. T. Copson 


Communicated by A. ERDELYI 


§ 1. Introduction 


One of the ways of solving the initial value problem for the equation of wave 
motions involves the use of spherical means. For instance, if u(x, y, z, ¢) satisfies 
V?u— 0? u/ot under the initial conditions u=f(x, y, z), du/dt=0, the spherical 
mean # (7, ¢) of w over a sphere of radius v and centre at an arbitrary fixed point 
(x, y, z) satisfies the equation 

Oru 2 Ou Cru 


ar? ' ~ Or at? 


under the initial conditions «@(7,0)=/(r), m%(7,0) =0, in an obvious notation. 
It readily follows that 


a(r,t) = ENG eV CI 


and hence that 
Paes a 
U(X, Y,2,t) = lim a (7, 0) == = {tf (d)} 
the well-known Poisson solution. The problem of solving the wave-equation is 
thus reduced to that of finding lim # (7, t) given u(r, 0). 


The same method can be applied to the special case of the Euler-Poisson- 
Darboux equation 


under the same initial conditions. In this case 


fu 20% Ou , 2 2% 


DVe ly 2OV ow Oe t, ot 


under the conditions #=f(r), #,=0 when t=0. The problem here has a sym- 
metry, since u, +0 as y->0. This suggests that 


ii(r, t) >i (t, 0) 


This research was supported by the United States Air Force through the Air 
Force Office of Scientific Research of the Research and Development Command, 
under Contract No. AF 18(600)— 1461. 


350 EE CORson: 
as y->0, and hence that u(x, y,z,t) =F(t), 


which is, ip fact, ASGEIRSSON’S solution. 
If one wished to extend this to the general Euler-Poisson-Darboux equation, v72. 


™m 


A24) a2, y ere 

Pu _ oe A ge (1.4) 
Oak. OR! Th et 

1 


one would to consider the equation 


eu , m—1 OU eu , Aa 
CO Gor ar Be 

and it is this equation which we consider here for general values of the constants A 

and m. It is convenient to change the notation slightly, and to discuss the equation 


eU | 20 0 aU , 28 au 
On? A dy ae Oy? Sy say 


where « and f are positive constants. Actually we shall assume that « and f are 
large enough to ensure the absolute convergence of all the integrals which occur. 
We can then deal with smaller values of « and f by analytical continuation, just as 
Diaz & WEINBERGER*® did in the case of the Euler-Poisson-Darboux equation. 

This partial differential equation looks rather special. The point about it 
is that it possesses two singular lines. Apart from the extensive work on the 
Euler-Poisson-Darboux equation, very little is known about partial differential 
equations with singular coefficients; a discussion of this special problem may 
serve to suggest the direction in which further progress may be made. 


§ 2. The problem 
The problem is 
To find a solution of 
#U | 20 OU 8U | 2p aU 
Ox? x Ox  ody2 ' y ay 


where «, B are positive constants, such that 
(i) U and its first derivatives are continuous in x20, y=O0, 
(ii) the second derivatives of U are continuous in x>0, y>0, 


(iti) U=}f (x) when y=0, x20; U=g(y) when x=0, y=0; / (0) =¢ (0). 


At first sight it looks as though there are insufficient data, since in a mixed 
problem (such as the vibrations of a semi-infinite string with a fixed end) one 
can assign two data on y=0, x20 and one datum on x=0, y=O. The fact 
that the data are given on singular lines of the differential equation upsets all 
this. By a simple extension of the arguments recently used by W. WALTER**, 
it can be shown that the problem has no solution unless U,=0 on y=0, x20 
and U.=0 on *=0, y20. Thus we have too much data, not too little. Moreover, 


if this problem has a solution, the energy-integral type of argument shows that 
it is unique. 


* Diaz & WEINBERGER: Proc. Amer. Math. Soc. 4, 703—715 (1953). 
** WaLTER, W.: Math. Z. 67, 361—376 (1957) 
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As things stand, the problem has no solution at all unless f(x) and g(y) are 
suitably related. This can be seen at once by considering domains of influence. 
For the data U=/(x), U.=0 on y=0, x>0, determine a solution in 0<ySx; 
the data U=g(y), U, =0 on *«=0, y>0, determine a solution in O<aSy. 
Continuity across the characteristic y = x leads to an integral equation connecting 
/(«) and g(y), and hence g(y) can be expressed in terms of f(x). This, of course, 
is precisely what one needs if one tries to solve equation (1.1) by spherical means. 


§ 3. Polynomial boundary data 
If /(x) and g(y) are suitably related polynomials, it does not follow that there 
is an analytic solution, still less a polynomial solution. Take the simplest case, 
when /(x) is a positive integral power of x, say 


te = reas os 


the coefficient being introduced to simplify matters later. The corresponding 
solution must be a oes function of degree N, and it is easily shown 
tober 


(«+ 1 Ae os aie 5 

memes cm aN F 5N,5-2N~«; B4 5354] (3.2) 

which, incidentally, automatically satisfies the condition U,=0 when y=0. 
There are two cases. If N is an even integer, (3.2) is a polynomial of degree 


N, so that this solution holds in the whole quadrant x20, y 20; and the value 


taken on x =0, y2Ois et , 
And we observe that, if >a, 
y 
y= Te rpresa | oer el aae | 64 
0 
an Bue ta P(b+s a 4a} 
t? g(t) = Fete! ee /(\t) (3.5) 


where J’~* denotes the Riemann-Liouville operator of integration of fractional 
order. 

But if N is odd, the hypergeometric series in (3.2) converges only if y?/x*<1, 
and so this solution is valid only in half the quadrant. It might be thought that 
we could continue the solution (3.2) analytically, regarding it as a function of the 
complex variable z= y?/x?. But the hypergeometric function has a branch 
point at z=1, and the analytical continuation (which is in general not real) 
depends on how one avoids this branch point—it is not a unique continuation. 
If, however, we start again in the case when N is odd, with g(y) given by (3.3), 
we get another solution: 


2) Se hlasert N Ss 5 ae ona 1 Re 
= ay (EN G-EN Bet 5 5) 69 


* The other solution of the hypergeometric equation, that of the type valid in 
|z|<1, 2° F(1+a—c,1+b—c; 2—c; 2), does not occur since it would make op 
infinite when y = 0. 
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valid in «2/y2<1. Although (3.6) is not the analytical continuation of (3.2), 
(3.2) and (3.6) together define a function which is continuous in the first quadrant ; 
and, moreover, the first and second derivatives are continuous across y= x for 
any integer N provided that «+f>1. Thus with polynomial or analytic data 
connected by (3.4) we can build up solutions with continuous second derivatives _ 
by using the polynomial solutions of even degree and the hypergeometric solutions 
of odd degree. 


§ 4. Continuously differentiable data 


If f(x) and g(y) are merely continuously differentiable functions, the previous 
method fails, and we have recourse to RIEMANN’s method. This gives* us 


U(X, Y) => [UVJat > [UV e+ 
+f f{vg,- U 4 +28 UV} dx + (4.1) 
AB 


i Laces UV, + 7% UV} dy 


DATA Bx Where V(x, y; X,Y) is the Riemann-Green function and 
Fig. 1, Diagram for Remann’s 4 and B are the points where the curve carrying the % 
method data is cut by the characteristics y—x=Y—X and 


ytx=—Y-+X respectively. 
Now if Y>y>0 and X—~x hes between +(Y—y), it is known** that, if 
6—¥ is not an integer, 


mo yfith Pts 


Vilaeye Xe Y) = 


cosaB  ya-byb-$ * 
x J A Ju—y(A*) Js (4X) {Ip-4 (4) Jy_p (AY) —Jp_a (AY) Ap (Ay) } aa. 


If we take the case 0< Y<X, the situation is as in Fig. 1, A being (X — Y, 0), 
B being (X + Y, 0). But as the datum f(x) is given on the singular line y=0, 
we have to consider the limiting forms of V and V, as y-+0. Since B>0, V tends 
to zero as y>0. If we assume that B>%, we find that, for all such Bp, 


(4.2) 


FN ie anaes 
y ai-B ed phd 


[ P*L104 40% baayyaa 


,) Piers (4.3) 
ln PG) xzytht" [#105 Bs Fe) 

where 

R = Y*— (X — x)* 


by a formula due to H. M. MacponaLp **%, This holds when x lies between X tY 
TO YS X; but if O0< XY, it holds only when x lies between Y+ X. Since 
V vanishes on y=0, (4.3) also gives the value of V, when y=0. 


x See, for example, eqn. (2.3) of the preceding paper by E. T. Corson. 
** ibid., eqn. (5.12). 
**kk See Watson, G.N.: Theory of Bessel Functions, p.411—412. Cambridge 1922. 
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It follows immediately that, if 28>1 ando<Y<X, 


MOI E 


Bee) a5 ae j [fo We F(a, ae hes we) d ween <(44) 


XY 


Similarly, if 2a>1 ando<XSY 


YerX 
e. D(%+) — 
err eaee | Ce ASB ea 4) 
Wave 


where R2= xX? (Y—y)2. 

We thus have one solution depending only on f(x) valid in 0<YSX, a 
second depending only on g(y) valid in 0< XY. We have now to see whether 
we can fit them together in such a way as to form one solution with continuous 
first and second derivatives valid in the whole quadrant. Actually, as we shall 
see later, it suffices to make U continuous. 


Now from (4.4) we have 


Ux, x)= TO 2-1 frexs ) wt8-1 (4 — e- F(a, VE ey rae 


1) 45%+26— 
_ L(p+ = 2 [rcxe Area a {2)8P—% Pl F(Z) dt 
wt 


and from (4.5) 


al 
1) ,2a+f8—-1 
LP(a+ = 4 { g (2X2) jxt+B-1 — f2)82—4 Pr) dt. 
Tt 


Ux, B= 


Thus / and g are connected by the equation * 


: 1 
2PI(B +4) S f(X 8 ete-1 (4 — Pye} PLPC dt 
A (4.6) 
= 2° F(a +4) fe (Xt) ere (yt Pt (2) at. 
0 


If «=, f(x) and g(x) are identical. When B>«, we should expect that 


y 


phd g(y) = Teh met yf PT) yt apm ads (4.7) 
0 


and this is, in fact, the case. It can be proved by substituting for g(X #) from 
(4.7) in the right-hand side of (4.6), inverting the order of integration and using 
the result ; 
ae 6(4 — #)ae-8 (A — 8-2-1 Pht) dt 
SP CNTR hohe? 2-1 (4 — 2 \hh 8 PEP (y),. 


* We have replaced 2X by X. 
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I can give no reference for this last identity. It can be proved by showing that, 
for all ea values of s, 


fics bs fp ACs 2yba— (¢2 — i) cai pad (2) dt 


UW 


1 
= Pep — a) f eet hte—3 (1 — a MPA PAP) di 


0 
and then using LERCH’s Theorem. We show in the next section that the matching 
of the two solutions along y= does suffice to imply the matching of the first 
and second derivatives. 
Since (4.7) can be written in the form 


10) tracers aioe pea fet 2)B— "se 1t(yt)d t 


the condition /(0)=g(0) of the problem is satisfied. Moreover*, if B—«21, 
this relation and the fact that f(x) is continuously differentiable imply that 
g(y) has continuous derivatives of order higher than unity, the actual order 
depending on the difference 6 —«. Conversely, if we have merely B —a >0, the 
continuous differentiability of g(y) does not imply the existence of a continuously 
differentiable {(x) satisfying (4.7). Thus if B>«, the continuously differentiable 
datum must be assigned on the axis of x. 


§ 5. The continuity of the derivatives 

In §4 we found two solutions, U, and U, say, of equation (2.1) valid respective- 
lyin 0<yS% and in 0<%S y, such that U, and -U, were identical on’ 4-7 
We now show that if «+ is big enough, the first and second derivatives of U, 
and U, are the same on y= x. 

If we introduce the characteristic variables r=x—y, s=x-+y, we have 
two solutions of 

(Uf | & : p Tea ( o : p : a 
ye (ots US sty tard . 

valid respectively on s>r=0 and on s> —rv=0O and identical on y=0. The last 
condition implies that 0U,/és and @U,/0s are equal on 7 = 0. 


Denote the values of 0U/ér and dU/és on r=0 by u(s) and v(s); then 
du | ba Dae a— Bp 


isan s s 


Un 
Hence, if ¢ is any positive constant, 
S$ 
u(s) st? — m4 (c) cP — (q — p) f +8-1 y (2) dit. 
Cc 


By hypothesis w(s) and v(s) are continuous in any finite interval O<s<S, and 
are therefore bounded. And as «+ 8>0, we can make c tend to zero, and obtain 


(9) tH — (6 —f) fev 9(9 at 


* I owe these remarks to Professor ErDELYI. 
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Hence Bois. oe 2 : 
[w(s) |i s7** = — (& — By ee v(t) |jdt =0. 


Hence 0U,/dév and GU,/d7 have the same values on y=0, so that the first deri- 
vatives of U, and U, are identical on y = x. 

Evidently the second derivatives 0?U,/dr és and ©2U,/ér és have the same 
values on y= 0 by the differential equation; so also have €2U,/0s? and €?U,/és?. 
To deal with 6?U/ér?, we have to assume the existence of third derivatives. 
It can be shown that if the value of d?U/é7? on y=0 is denoted by w(s), then 

dw | a+b, &-f+o-B , (@—ptatB, 


ei | 
ds Ss s2 s2 


Hence 
w (s) st*P — aw (0) oP + f PPR (a? f+ —B)u() + [2 —f)? +4-+flo@}al. 


If «+f>1, we can make c tend to zero provided that we assume w(c) bounded 
as c->0. The equality of the values of the second derivatives on 7 =0 will then 
follow provided we pay the price of further assumptions on the second derivatives. 
Fortunately, as we show in the next section, it is not necessary to do this. 


§ 6. An alternative method 


An alternative method of dealing with the case Y>X of §4 is suggested by 
that used for solving the wave equation 


Oru anu 
Ox? ay? 
when the data are a 
| —— 4 a 2° — x = 
Up Ke), By g (x) Vie A0 
u = F(y) Sh (0) eee Os 
We suppose that the data imply that 
Ou Wik Ss 


where G(y) is an unknown function. For such a solution 
fu,dx +u,dy =0 
net 
for any closed curve I’. If we take J” to be OBPAO of Fig. 2, we readily find that 


YoRX Ya 


u(X,Y)=2/(Y+X)+3F(Y—4)4 eS BME G(y) dy 


which involves the unknown function G(y). But if we take I" to be the triangle 
OCA, we get 


Y=X Y- 
=3/(¥Y—X)—sF(Y—X)4 a J pe des G(y) dy. 


And, if we subtract, we find that, when Y>X>0, 
Yon 
u (X,Y) =2 (V4 X)—F/(V— XV 4 FV — X44 S eae. 
V=¥ 
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For the equation (2.1), everything turns out to be even simpler. We assume 
that Y>X>0 and apply formula (4.1) to the contour OBPAO of Figure 2. 
By (4.2), V behaves like a multiple of x** when x is small; hence if we 

assume that «>, the integral along AO and 
: P(X, Y) the point A make no contribution, and we 
are left with 


(YHA u(X,Y) => [UV]a+ 
Ea i Erica 2p \ 
+4 f{rg UY, + =F UV hae. 
OB 
0 C1V-X,0/) QGUXO x : 
Fig. 2, Contour for the second method As before, V tends to zero as y tends to zero, and 
V 


a te Sc a ae aan eed ae 


if B>}%. When Y—X<x<Y-+X, the value of this integral is given by (4.3); 
but when 0<%< Y — 4, 


eae 2 (6 —4) aed Oia air tail eB 
CES Nie yo 
a—-B+2 a—B+1.. 1. 4x2 X? 
x F( ae ey oe eh eet eae 


It follows then, that, if a and # are greater than $, the value of U(X, Y) 
when Y= X>Ois 


yen 
U(X, Y) = 2r'(6 +3) wl ) x22 X2_ y2)P-o-1y 
ae P(a+4)P(B—«) VP | ae a 
0 
Eee eee eee 
x F| : : Fate Sg ogee Tee (6.1) 


| C(p+3 2) x x 2p—-1 on RS 
Sern ee [ fe RiP F(a, 1— a; Bs Gte| dx. 
Y-—X 


This is the required solution when Y= X, and should be compared with (4.4) 
Gig SO. 


Lastly, when B>«, we may make X tend to zero in (6.1). It is readily seen that 
the second line tends to zero, and we are left with 


EU rae, alee, APT fr ) 2% (V2 — x2)P-a-d gy 


which is, in fact, precisely equation (4.7). 
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Abstract. This paper! develops, with an eye on the numerical applications, an 
analogue of the classical Euler-Cauchy polygon method (which is used in the solution 
of the ordinary differential equation 


for the solution of the following characteristic boundary value problem for a hyper- 
bolic partial differential equation 


Ugy = [(4, 9, U, Ux, Vy) , 
u(¥, Vo) =9(*), 
u(%,¥) = ty), 


where o(%)) =T(¥Vo). The method presented here, which may be roughly described 
as a process of bilinear interpolation, has the advantage over previously proposed 
methods that only the tabulated values of the given functions o(¥) and t(y) are 
required for its numerical application. Particular attention is devoted to the proof 
that a certain sequence of approximating functions, constructed in a specified way, 
actually converges to a solution of the boundary value problem under consideration. 
Known existence theorems are thus proved by a process which can actually be em- 
ployed in numerical computation. 


1 This paper was issued on 16 January 1957 as NAVORD Report 4451, U.S. 
Naval Ordnance Laboratory, White Oak, Maryland, and was presented to the Ameri- 
can Mathematical Society in October 1956. 
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§ 1. Introduction 


The classical initial value problem for the ordinary differential equation 


dy _ 
Vira =aif (ey y)s 


(where the real valued continuous function {(x, y) is defined for XpSxs Xp ta 
and — c<y<+ 00) consists in the determination of a real valued function y(x), 
defined on x)<x< % +4, which satisfies the given ordinary differential equation 
on this interval, and also satisfies the initial condition 


¥(%) = Yo» 
where yp is a given real number. 

Among the many methods which have been employed for proving the existence 
of a solution y(x) to this problem, mention will be made here only of PIcarD’s 
method of successive approximations (see e.g., G. SANSONE [2/, vol. I, pp. 9—14], 
E. L. Ince [12, pp. 6365], E. A. Coppineton & N. Levinson [28, p. 11—13], 
or E. Kamke [16, pp. 54—56]); of L. ToNELLI’s method (see, ¢.g., L. TONELLI [13], 
G. SANSONE [2/, vol. I, pp. 45—48]); and of the Euler-Cauchy polygon method 
(see, ¢.g., G. SANSONE [21, vol. I, pp. 36—45, vol. II, pp. 208—283], E. L. INcE 
[12, pp. 75—81], E. A. Copprncton & N. Levinson [28, pp. 3—7], E. KAMKE 7 
[16, pp. 62—64], or G. A. Buiss [9, pp. 86—92)). 

For the numerical purpose of the actual construction of a solution the Euler- 
Cauchy polygon method is usually the most advantageous. The construction of 
the Euler-Cauchy polygons may be described as follows. For each positive in- 
teger m, let 

xo = Xom—< Xm Xo om TOES None Xm,m = *o =F a, 


be a subdivision of the interval x»<*<x,+a into m closed subintervals 
Ky mee XS X44 m, Where k=0, 1.0, m—4.)\ On each such subintervaly the 
ordinary differential equation is, so to speak, replaced by one whose right-hand 
side is a (suitably chosen) constant, so that the corresponding function ap- 
proximating a solution turns out to be a linear function on each subinterval. 
More precisely put, the polygonal function y,,, which is an approximation to a 
solution, is defined recurrently by the equations 


dy : 
a (x) = T\Xorns Yo): Vn (Xo, =) == Yoo on Xo, n= Ss x = Xt m> 
dy 
a te) = f(%4, m3 V1), Vin (Sa) SiG on I pe SS Xe es 
dy i : : 
Pe (x )= f (%., m? Vis Vn (Xb on) = Veo on Xk m = x S Ne+4,m , 
for k=0,1,...,m—414. Notice that, for simplicity in writing these equations, 


the symbol y, is used to denote the value of the function y,, (x) at x, ,,, a value 
which is obtained from the definition of y,, as a linear function on the preceding 
subinterval %;_4 << %,»,,, and which is used as an initial value for the function 
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¥(x) for the “miniature”’ initial value problem (of the same kind as the original 
one, but whose differential equation has a constant right-hand side): 


wai (x) =1(%p ms Ve) Vn (Xk, m) Sy he 


on the next subinterval x, ,<*< %,11,,,- For each positive integer m, the func- 
tion y,,(”) is continuous on the interval x»<*<x,)+a, but its derivative will, 
in general, not exist throughout the interval, since it may jump at the subdivision 
numbers Xj in- 


Under the sole additional hypothesis that the function f(x, y) is bounded in 
absolute value on %<%S%+a, —co<y<+ oo, it follows that the sequence 
of functions {Vm (x)} is equibounded in absolute value and equicontinuous on the 
interval x» x<%)+a, and hence, by AscoLi’s theorem [1] (see also TONELLI 
[11, p. 70—86]) there is a subsequence of the sequence {5m (x)} which converges 
uniformly to a continuous limit function on x, x<x,+a. If, further, it is 
supposed that the maximum length of the subintervals of the subdivision of 
Ko %S%y+ a approaches zero, 1.¢. 

PO MAX PHN ory — 0, 


m—>co +k=0,1,...,.m—1 


then every such continuous limit function is a solution of the original initial 
value problem, whose solution need not be unique. (It should be noticed that 
the condition on the maximum length of the subintervals is automatically satis- | 
fied in the most common case when the m" subdivision consists of m subintervals 
of equal length, namely a/m.) If, besides this, the function /(x, y) satisfies a 
Lipschitz condition with respect to y, 7.e. there is a number L=O such that 


| f(x, V1) — }(%*; V2) =1\y,— y» 


whenever %S*X% +a, then the whole sequence {Vn (x)} converges uniformly 
on %SxS% +a to the (known to be unique) solution of the original initial 
value problem. 


? 


The purpose of the present paper is to develop, with an eye on the numerical 
applications, an analogue of the Euler-Cauchy polygon method for the solution 
of the characteristic boundary value problem for the hyperbolic partial dif- 
ferential equation 


Ux y = I(x, y, u, Ux, Uy) ; 


(where the real-valued continuous function f(x, y,%,~,q) is defined for all 
(x,¥,2,~,q) satisfying 
atau ta, WSVsyoto, and —~cx72,p,9— + ~). 


The problem in question consists in the determination of a real-valued function 
u(x, y) which satisfies the given partial differential equation on the rectangle 
Xp SxXK Xp ta, V—ySVSv+4, and also satisfies the conditions 


u(x, Yo) = 6 (2) for %»S%x*S%+4, 


U(X», Y) = T(y) for Vy SV SV 94-18, 
Arch. Rational Mech, Anal., Vol. 1 se 
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where o(%))=T(Vo) and a(x) and t(y) are given continuously differentiable 
functions on the characteristics y=) and x=%9 of the given hyperbolic 
equation. (The treatment of this boundary value problem by successive approxi- 
mations goes back to E. Picarp [5] and has been considered by various other 
methods by many writers since that time.) For each pair of positive integers m 
and », consider the following subdivisions of the intervals 


KoSxSx%y +a and yySySyt+, 
- — | 
xo = Xo,m < *1,m =< Xom ae < Xm —-1,m <*n,mn = xy 7 4, 
— | 
Vo = Vom Vin = a ee Vi eae Vn,n = Yor b, 


which produce a subdivision of the rectangle x» xS%)+a and yy ySyo+ 0. 
The miniature problem in the present method (see Section 3 for details) is as 


follows: 
Um n 


ax ay (x, ¥) = Agi, FOL 4, Se Nee, Ores 


Umy(% V1) = Dei t Baile —%), for %, SxS x41, 


Unn (Xn ¥) = Dait Caily — i), for y;SySV41; 


where A,,, By;, Cx, and D,,; are suitable constants, depending on the sub- 
rectangle (for simplicity in writing, x, has been written for x, ,, and y, for y; , 
in the formulation of the boundary value problem for the subrectangle). This 
means that on each subrectangle, the approximating function w,,,, is bilinear in 
(x, y), te. it is a hyperbolic paraboloid: 


mn 


Un (X>V) = Api e— 2) (Via 9) oe Bale = e—) Ca ea 


The process just described reduces in the special case of the equation u,,—=/ (x, y, u) 
and equal subdivisions of the intervals x»SxSxjt+a, yoSvSvot+b to the 
process given by G. ZWIRNER [24, pp. 222—223], who did not consider the more 
general equation treated here. Similar methods, analogous to the one described 
above, have been employed to prove existence theorems for the same boundary 
value problem by P. HARTMAN & A. WINTNER [26], R.H. Moore [29] and 
R. Contr [27], but they do not appear to be as convenient for numerical pur- 
poses as the one described above, which requires knowledge only of the tabulated 
values of the given functions o() and t(y) (from which the difference quotients 
needed may easily be calculated) and does not require the tabulated values of 
the first derivatives o’(x) and t’(y). Mention is also made of a different, but 
closely related, method, also analogous to the Euler-Cauchy polygon method, 
given by H. Lewy [/4] (see also H. BEcKERT [22]) for the solution of the initial 
value problem for second order quasilinear partial differential equations in two 
independent variables, which appears to require more differentiability assump- 
tions than the present method. 


The statement of the known main results and their connection with the existing 
literature is given in Section 2. Section 3 contains the precise description of the 
analogue of the Euler-Cauchy polygon method and the construction of the double 
sequence of functions {w,,,,(x, y)} approximating a solution. Each function w, , 
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is continuous, but not necessarily differentiable with respect to + and y on the 
rectangle %S%5%+4, YySySyot+b. Section 4 contains an inequality, termed 
the convergence inequality, which is used, together with a theorem of C. ARZELA 
[7, pp. 119—125] on the convergence of certain not necessarily continuous 
functions to continuous limit functions, in order to complete the proof of the 
existence of a solution in Sections 5 and 6. 


§ 2. Statement of known results 
Theorem 1. // 


(1) the real-valued function f(x, y, 2) is defined for all (x,y, 2) such that 
MySXS%+A MZSVSM+t), -—w<z<to, 


where Xo, Vo, 4, 6 are real numbers, and a=0, b=0, and if f(x, y, 2) ts continuous 
and bounded in absolute value, so that for a certain non-negative constant M one has 


|7(*,y,2)| SM 
for all these (x, y,z); 


(2) the real-valued function a(x) ts defined for all x such that xySxSxy+a 
and possesses a continuous first derivative o'(x) for all these x, while the real-valued 
function t(y) ts defined on the set VyoSySyo+b and possesses a continuous first 
derivative t'(y) for all these y (it being understood, of course, that o'(X»), for example, 
denotes the right-hand derivative of o at Xo, etc.); then 


(3) there is at least one real-valued function u(x, y) defined on the rectangle 
Ri %SxXS%+a, YSySyot+4, 


which is continuous, together with its partial derivatives Oulox, Ouloy, Puldox oy 
(= @ul/dv 0x) on R, satisfies the partial differential equation 


nae (x, 9) =t(x,9,4(%,9)) for (x,y) in R 
and the characteristic conditions 
UK No) = 0 (%)- ‘fOr xB K SX 4, 
U(%,V)=T(y) for YWSVvSyot d. 


It is to be noticed that this theorem asserts the existence of at least one 
solution to the characteristic initial value problem under consideration, but that 
the uniqueness of the solution is not asserted, and is, in fact, in general not true. 
(See P. MonteEx [8, pp. 279—283].) One need only consider the following simple 
example of a characteristic problem (cf. P. HARTMAN & A. WINTNER [26, p. 84] 
and P. LEEHEY [23, p. 23]) consisting of the partial differential equation 


2 
Ot 8. | yf for O<x<a, 0XySXb, 


U(X, Oj 2 fOr 0S % = 4; 


MONV)— 01), S100 0 SHS D, 


362 eS wOTAZ: 


a 


which has as solutions both 
i ke) = 05 
and 
Uy (x,y) = te x? y?, 


on the rectangle 0S *Sa, OS ySb. 


Theorem 2. [/ 
(1) the real-valued function f(x, y,2, p,q) ts defined for all (x, y, 2, p, q) such 
that 
HyeSKS MEA, WSYS Mtb, —OK<2z,p,q9<+0, 


and is continuous and bounded in absolute value, so that for a certain non-negative 
constant M one has 


l{(~y2p9| SM 


for all these (x, y, z, p,q), and if f satisfies a Lipschitz condition in the three argu- 
ments z,p,q (that ts, there is a constant LZO such that one has 


f(4,9.2,0.9) —Kea Pr. Gs lke—a_|+Llp—fp)4+Llla—qa; 


for any (z, p,q) and (2%, p,, %), whenever (x, y) lies in the rectangle 


Ri S25 %6454, Y Ves VSVo Os 


(2) the real-valued function o(x) ts defined for all x such that x»xSxSx%y+a 
and possesses a continuous first derivative o'(x) for all these x, while the real-valued 
function t(y) ts defined for all y such that Vy vSyo+b and possesses a continuous 
first derivative for all these y; then 


(3) there 1s one and only one real-valued function u(x, vy) defined on the rec- 
tangle R, which ts continuous together with its partial derivatives 


Ou OU Cru i= C7 u 


ax’ ay’ ax dy et, 


~ Oy Ox, 
satisfies the partial differential equation 


Oru 


ou ou ; : 
ax Oy (x, ¥) = Hx, 9, u(x, y), : (x, ¥), (x, y)| for (x,y) in R, 


Ox 


and the characteristic conditions 
U(X, Vo) == Ox) > VfOr x, See, ae ee 
u(%,v)=T(y) for ySyvySyot+. 


This second theorem does not contain the first theorem as a special case, 
since the function /(x, y, z) of Theorem 1 is not assumed to satisfy a Lipschitz 
condition in the argument z. However, if in Theorem 2 the function P(X AVSS, PEG) 
does not depend on f and g, then Theorem 2 yields the additional information 
that if f(x, y, z) of Theorem 1 does satisfy a Lipschitz condition in the argument z, 
then the solution whose existence is assured by Theorem 1 is indeed unique. 
Theorem 2 is the classical theorem of Picarp [5] mentioned in the Introduction. 


The Polygon Method 303 


Theorem 3. I} 


(1) the real-valued function f(x, y, 2, p,q) is defined for all (x,y,z, p,q) such 
that 


Koh S% +a, VYSVSYytl, —~<2,h,9< 0, 


and 1s continuous and bounded in absolute value, so that for a certain non-negative 
constant M one has 


lf ,2,0,9)| <M 


for all these (x, y, 2, p,q), and if f satisfies a Lipschitz condition in the two argu- 
ments p,q (that is, there is a constant L=O such that one has 


[f(%9.2,0,9) —1(% 92%, b.n)|SL|p—A| +Llq—a| 


for any (p,q) and (p,, q,) whenever (x, y) lies in the rectangle 


Ri MSxX*S%ta, YSySyt+), 
and z 1s any real number); 


(2) the real-valued function o(x) is defined for all x such that 
ho % SX + a, 


and possesses a continuous first derivative o'(x) for all these x, while the real-valued 
function t(y) ts defined for all y such that 


Y=VSY¥ot 6 
and possesses a continuous first derivative for all these y; then 
(3) there is at least one real-valued function u(x, y) defined on the rectangle R 
which ts continuous together with its partial derivatives Ouldx, Ouldy, CulOx doy 
(=@uldy 0x) on R, satisfies the partial differential equation 
4 (x, y)| for (x,y) in R, 


Cu ou ra) 
x oy 


an oy (x,9) = Hx, 9, u(x, 9), ey, (x, 9), 


and the characteristic conditions 
UNL, Vo) =O (2) Je JOT Noa XS Kg a, 


U(%,y)=t(y) for YoS VS Yo + b. 


This third theorem contains the first theorem as a special case (and the 
same example used there is applicable here). The hypotheses made in the third 
theorem are such that the part of the second theorem concerning the existence 
of a solution follows, while the second theorem yields the additional information 
that if the function f(x, y, z, p,q) satisfies a Lipschitz condition in (z, A, q) 
together, rather than just in (f, g), the solution u(x, y) whose existence is asserted 
by the third theorem is indeed unique. Theorem 3 was first proved by P. LEEHEY 

23] and P. Hartman & A. WINTNER [26]. For more general theorems see 
R. Conti [27] and A. ALEXIEwIcz & W. ORLIcz [30]. 
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§ 3. The double sequence of functions approximating a solution 


Let m and n be positive integers and consider the corresponding subdivisions 
of the intervals x,<xSx)+a and ySySyo+8), as follows: 


Xo = Xo, < 1m < %a,90 <9 °* S Sn —tym = Fmm = A071 4 
Vo = Von <Viyn <V2,n <0 Vpn S Van = Yo + b. 


These subdivisions of the intervals x,< xx +a and yy ySyo+0 produce a 
subdivision of the closed rectangle R into m-n closed subrectangles Rj", where 
k=0,1,..-,m—1 and /=0,1,...,%—1. The closed subrectangle Rj" consists 


in all (x, vy) of R which satisfy the inequalities 
Ry KS Xia m Vi gd San 


Given the functions o (x) and t (vy), defined on the closed intervals %jS7S%9+ a 
and yy<ySy_o+b respectively, a continuous function w,,,,(x, y) will be defined 
on the rectangle R by a recurrent process, consisting in solving, on each sub- 
rectangle Ry", a boundary value problem of the form 6? w,,,,/0x 0y = constant, 
with assigned (linear) values for w,,,, on the two rectilinear closed intervals of 
the boundary of R;'” which intersect at its lower left hand vertex (%, ,,, Vi, »)- 
Of course, the constant involved in the partial differential equation, and also 
the linear boundary values, both depend on k and / (and on m and n). The fact 
that two adjacent rectangles, say Ry" and Rj", ; for instance, have a common 
boundary interval (since they are both closed subrectangles) will create no dif- 
ficulty concerning the definition of the function w,,,, for poits lying on the 
common boundary intervals, since the specific process employed in defining u,,,,, 
will be such that the values assigned to w,,,, will coincide in this situation. 

Suppose, for the moment, that w,,,, has already been defined on the sub- 
rectangle of R with lower left vertex (x9, yo) and upper right vertex (%, ,,, Vim) 
7.e., the subrectangle defined by the inequalities 


XS x Ss Xm Vos My = y 


Ln? 
where 15k<m—1 and 1</<n—1. Then the definition of the function Voy, 
will be extended to the slightly larger subrectangle defined by the inequalities 


XS XS Xp i1m Vos V Sa per ses 


by first defining it on the closed subrectangles 


Ren mn mn 
k, 0» Ryda sey k,l-1 


in numerical succession (7.e., passing from Rz'} to Ry1, and so on); then defining 
it on the closed subrectangles 


mW mn mn 
0,1 > i Pee Aan k—-1,l 


in numerical succession (i.e., passing from Ro; to Ri'7, and so on); and finally 
defining it on the remaining closed subrectangle R%” in order to complete the 
definition of w,,, on the rectangle 


KS VS ts) Voss Vi+4,m- 


f 
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(A simply drawn figure will readily make the process intuitive to the reader.) 
Alternatively, the function w,,,,(%, y) may first be determined on the m sub- 
rectangles in a row: 


nn mn mn 
0,17 3 IATA es OSS m—1,1> 
for the rows /=0,1,2,...,—1 in succession. There remains only to make 


precise just exactly what boundary value problem, 7.e., what partial differential 
equation and what boundary conditions, is to be solved on each subrectangle 
Ry. This will be done by showing how the process is started in the initial sub- 
rectangle Roy and how the step-by-step scheme indicated above can then be 
carried out, using the given data, the given functions o(x) and t(y). The final 
result will be an explicit formula for w,,,,(”, vy) at any point (x, y) of a typical 
subrectangle R;””. 


On the rectangle Rog the function w,,,, is required to satisfy the partial 
differential equation (with constant right-hand side) 


e Uma 


(¥ m — O\Xom Vin) — T\Von . mn 
“Ox Oy (x,y) =f (Xo, Vo» (ep) was im) 8 (tom) 71049) Un) | for (x,y) in Roo, 


, 
4m — om Vin — Von 


subject to the boundary conditions 


O.(Bn i) 0 (Xo) 
Um (% Vo) Aes 6 (Xp) al er eerie ie (x —_ Xo) for XS %S Mm 
nN 


(V,n) — 1(Ie) 
vA Vo 


Umn (Xo. ¥) = T (Yo) + (yv—%) for ySySy 


n° 


Roughly speaking, what is done in defining w,,,, on Rog is to take as boundary 
conditions along its left boundary edge and its lower boundary edge certain 
linear functions derived in a natural manner from the given functions t(y) and 
o(x), and to use the value of u,,,, at (%9, Yo) and the slopes of these linear func- 
tions in determining the constant value to be assigned to é?w,,,/0x dy on Roy. 
It is clear that, the boundary value problem for u,,,, on Roo being explicitly 
solvable, 


Umn(% y) >r Unn (%o, Vo) ae 


Umn (4%, m’ Vo) — “mn (%0 ¥o) (x %o) ++ 
Xm — Xo,m 

Umn (qe Y1,n) —Umn (%, Vo) (y Vo) | 
Vi,n — Vo,n 


Umn in (4%, mr Mesa Umn (Xo, Yo) 


te i(r. Yoo Umn (Xo. Yo) oi i a , 


%1,m — Lege 


Umn (Xo; V1, n) iz eet (x 
Vin — Vo,n 


Xo) (Y — Yo) 


for (x, y) in RQ’, where, for uniformity in the writing of formulas to appear 
later, u,,,(%9» Yo) has been written instead of o(%9) or T (Yo) eéc. It is to be noticed 
that w,,, is bilinear in (x, y) on Roo’, @.e., it is linear in x for each fixed y and 
linear in y for each fixed x. (From ae point of view the process of defining 

u,,,, being described may be thought of as a process of bilinear interpolation, 
so to speak.) 
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wo ae 
Consider now the definition of w,,,, on the rectangle R;z'", it being assumed 


that u,,, is already known as a linear function on the left boundary edge, where 
x = X,.,, and on the lower boundary edge, where y= y,,,,, of the closed rectangle 
Tah, Then u,,,, on the rectangle Rj", is required to satisfy the partial dif- 
ferential equation with known constant right-hand side 


mn 


, 


F Umn(Xh+i,ms Vin) — Umn (hm: V1,n) " 


(x, y) = fm’ Vin» Mea dito Nie) ) 


Xet+1,m— “kim 


PUmn 
Ox OY 


Umn (Xk, mi Vita,n) — Un (%k, ms 7a!) for (x,y) in Ry” 
EEE se £ fi 
Viti,n — Vin 


and to coincide with an already known linear function of y on the left boundary 
edge, where x= ~,,,,, and with another already known linear function of x on 
the lower boundary edge, where y=y,,,. It is clear that w,,,, is bilinear in (x, y) 
on Ry” and that 


Umn(%, y) = Pen Cp Vin) a 


Umn (X41, m3 Via Umn (Xk, ms Vin) (x eae ) ue 

4k+1,m — Xk om bY 
Unn (¥p,m3 Vi-+1,n) = Umn (%e,mi Vi,n) ( 
Wein ype Gy tert 


Vea V1 n) ot 


‘ : : . Unn (¥pit mi Y1,n) = Umn (Xk, m3 Vin) “ 
ole He m' Vin» Umn (Xk, m els } ~ ee } 


/ 


? 


Xk+A,m Xk, m 
Mm n (%b, mi Vt-+3,n) — Monn (%e, mi Vt,n) (x—2» m)(Y¥—%n), for (x,y) in Rey". 
Viti,n — Vin 

This last formula does not exhibit the explicit dependence of the function 
Um»(X, Y) On the given functions o(x), t(y). In order to obtain a formula which 
makes evident this explicit dependence on o and 1, which will be essential in the 
convergence proofs to follow, it is convenient to use an abbreviated notation 
yielding more manageable formulas. For example, when considering the function 
Umn, With m and m regarded as fixed throughout the discussion, a functional 

value such as 


Uy, n (Xp, m ; Me by) 


will be denoted simply by w,,, and a functional value such as 


Umn (Xe+1,m3 Vin) Te Umn (Xeni Vin) \ 


> 


H(a,m' Vin» Unn (Apes Virgil F 
Vk+1,m — Xm 


Umn (k, mi Vi+-1,n) —Umn (%k, ms Vin) 


Viti,n — Vin 


will be denoted merely by f,;. Further, x, and y, will replace x,,, and ¥,,,, 
respectively. 


In this notation, the above formula for w,,,,(x, y), for (x, y) in Rp", may be 
rewritten 
ae 


a. Un+at Un Iti — Up 
Mm n (%, y) = Uni + pana (%— 4%) + eS — (y— 1) +her (x — x) (¥ — 1). 
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Using this abbreviated notation, one has the following formulas for w,,,, on each 
of the rectangles Roo’, Rio, Roy’, and Ry’, which are all special cases of the 
last formula just written for Rj7”, for (k,/)=(0,0), (1,0), (0,4) and (1, 1), 
respectively. In the first place 


iyo fey) ie OS Mater yi) tg (x — a) A) 
x1 — Xo Vasa 0 


for (x, y) in Roo, that is, when x» xXx, and yySySy,. In the second place 


eet 


eg ee (fg ee (ye at file — XN) vale 


ve — X%y Sir — Se 


for (x, y) in Ryo’, that is, when x,Sx<x, and y<ySy,. In the third place 


U4, — U 
Bee Mel ee Toul = to) (at) 


LOF\V1y) eno that is, when 7, = 7=%, and ¥,=y=¥,5- In the fourth place 


yy (yn) ae en) 


X_— Ay Vos Vy: 
fori(4.y) ini, , that is; when 7, =+= 7% and 4,2 y=. 
The formulas for Roy’, Ryo and Ry" will now be rewritten so as to reveal 


the exact influence of the given functions o(x) and t(y). From the formula for 
(x, v) in Rog it follows that 


M11 = U9 + M91 — Moo + foo: (%1 — Xo) (Yi — Yo) - 


i a mn 


Substituting this expression for “,, into the formulas for (x, y) in and Roy’, 


one obtains 


Uog — U u U 
ogy (2, 9) = thay + M9 M9 (x — xy) + Hoe Mee (y — 9) 4 
2 1 1 0 


+ foo: (%1 — %o) (Y — Yo) + fro: (% — 1) (Y — Yo); 


when (x, y) is in RY’, and that 


tao ton y\ 1 Moe M1 
eu (% — %) 4 ya, (y — 4) 4 


Unn (x, y) = a Uo1 aie 


+ foo: (% — %o) (¥1 — Yo) + for (% — %o) (Y — M4)- 


when (x, y) is in Roy’. 


mn 


Now from these last two formulas for (x, y) in R75’ and Foy’ one obtains 


Ug = Ui + U29 — Moot foo (%1 — Xo) (V1 — Yo) + ho (%2 — %1) (Yi — Yo) » 
and 


Uy 9 = U9 + U2 — M044 foo (%1 — *o) (V1 — Yo) + for (41 — %0) (Y2 — 1); 


these, together with the already known equation 


Uy = M9 + Mp1 — Moot foo (%; — %) (Vi — Vo)» 
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can be used to rewrite the formula for (x, y) in Ry’ as follows: 


Uso0— “0 , Uno — Yor a 
Umy (X%,V) = M19 + Mor — “007 <7 (% — %) + Aree ase (3a) ae 


+ foo: (%1 — %o)(¥1 — Yo) + ho: (% — %) (1 — Yo) + 
+ for: (%1 — %) (Y — Vy) + hile —A) (—W"), 
for (x,y) in RY. 


From the preceding considerations, the following general formula may be 
obtained by a process of mathematical induction: 


Upn+19— Uk, 0 Uo 141 — U,1 
Umn (X,Y) = Up,9 + Uo,1— Yoo ep ae (% — x,) + ars aan (ia 
Ea ed) l 
: ie 1,j 1 (%; — %;—-1) (Y =, 44) ay A —1(% =p) (Vg eg) ee 
4=1 7=1 j=1 


a >; fi-a,o(% = 4%.) (9 —F) a 


for (x, y) in Ry", that is, when x,<*<4%,., and ySySy4,, where k=O, 1, 
.,.m—1.and 1=0,1,...,%—41. It is readily seen that by putting (, /) equal 

t0;5(0;,0) 7043 0); (0, 1), (14, 1) in turn one obtains the formulas given above for 
60> 10> or > 11 > respectively, as special cases. 

For each pair of positive integers m and n, there has been defined a sub- 
division of the rectangle R into m-™m closed subrectangles, and there has also 
been defined on the rectangle R a real valued continuous function 4,,,,,(%, ). 
This double sequence of continuous functions {nalts y)} is equibounded in 
absolute value on R. For let A, B, C, D denote non-negative real constants such 


that 
lo(x)| SA, |o(x) —o(«*)|<Clx — x*], 
whenever %)S*3%)+4, and %Sx*S%,+a; and 
lé(y| SB, |r(y)—t(y*)|SDly —y*|, 


whenever yop ySyo+b, and yy<y*<y +b. (The existence of these constants 
A, B,C, D follows from the assumptions made about the functions o(x) and t(y) 
in any of the three theorems of Section 2.) Then, given (x, y) in R, one has 
XpSXS%,4, and ySySy,,, for some suitable pair of integers k and 1, with 
OSkSm—i1 and 0S</<nu—1. Hence 
mn #901 S|, +| toa +] Mae| + eee A fa — ag +) OAS yy 
| +t cd | 
l 


k l 
pa a (%; — %j—1) (Vj; — Ha) + 2 M(% — %) (vj — 9-1) + 
i=1 j= j=1 
k 
a 
=1 


M(x — x1) (Yn) + M(x —m) (yy), 


where M=0 is an upper bound for the absolute value of the function f (see the 
hypotheses of Theorems 1 to 3). Thus, by use of the definitions of the constants 
A, B,C, D just given, it follows that 


|%mn(¥,¥)|S24+B+Ca+Db+Mab, 
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where the numerical constant on the right-hand side does not depend on the 
point (x, y) of R or on the pair of positive integers (m,n). This proves that the 
double sequence of continuous functions {nn (x; y)} is equibounded in absolute 
value on R. 


Now, for each pair of positive integers (m,n), let the positive numbers « 
and fB,, be defined by 


“mn 


a» — max Ae Ny. 
Hk cee ook ®) 


and 
B,, 


Bi ie ees. ial V1) , 


so that the product «,, - ,, is certainly not less than the area of the largest sub- 
rectangle of the subdivision of R corresponding to the pair of positive integers 
(m, n). Under the additional restriction that 

lin-oe = 0 Granda slim ps0 


m— Co 


(which implies, but is not implied by, the fact that the maximum area of the 
largest subrectangle of the (m, )'* subdivision of R approaches zero) it will be 
shown that the double sequence of continuous functions PAGS, y)} is an equi- 
continuous double sequence of functions on R. By this is meant that if {thing ng(%s y)} 
is any singly infinite sequence of functions (with bm mM, = lm Ng oo) extracted 


from the double sequence {w,,,,(x, y)}, then the set of all functions u,,,,,, where 
y=1, 2,3,..., is an equicontinuous set of functions. 

In order to show this, one has to find an upper bound for the absolute value 
of the difference w,,,,(%,V) —Umn(*, y), where (¥,¥y) and (x,y) are points of R. 
There are really four cases to consider, depending on the relative positions of 
the points (x,y) and (x, y) with respect to each other namely; *<¥% and ySy; 
%Sxandysy;*s%*% andysy; *Sx and yxy. The first two cases are essen- 
tially the same by symmetry, 7.e. by interchanging the roles of (x, y) and (%,¥), 
and a similar remark applies to the last two cases. Only the first case mentioned 
will be considered here, since the treatment in the third case is exactly analogous 
to it. In the first case one has x=%, ySy and %,S%*%S%,44, Vi SVSYi41, and 
XpSXS R44, VISVSy41 for suitable pairs of integers (Fk, /) and (k, 7). Further, 
%ySxR, WSyy and %yiy SX, Wi SV1- 

From the definition of the function w,,,, it follows that 


Un, AC y) — Unyn (%, y) = Uz 9 — Uz, 0 Uo 7 Uo 1 


Uk+1,0— “Yk,o ;= , “o,f+1 > Y40,t (= ¢ 
To _ j XR 7 ——(y—y 
Xki1— ¥k ( ®) ye ieee). i ( y 1) 
Ui, 01a 4k, 0 (x t x;) Ley CAN) @ on y)) a 
*pti— *k as We = yr i 
k I k l 
An D>; ya =e Ds [fs j= 1 (%j — %- Oe) a) er 
4=1 7=1 4=1 7=1 
if 
+ D fej We AY, = ya) + Dhaai A) (yy 
q=1 
l 
= Lbs, i ees: Xp) oe Maps 4) ms fa % ig Xj 4) (v “+i Vi) Sie 
j= Oa 


~. 


 fey(% — xg) (¥ — v1) — feil® — *) (¥ — )- 
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Hence 
| nn (x, y) a Unmn(%, y)| a C( x5, ae Xp) a D (yt41 er y)) ois 20 Cpt pb M os oF 
4M -[ (x — x9) (vr— Yo) — (44 — %0) (Yi — Ma) | + 2M Oa, + 2M 0B, + 2M om Bu 


in terms of the constants A, B, C, D, M which were introduced earlier. However, 


gpa — Hy S (Et Og) — (4 = Hp) = F — 2) + Lot, 
and similarly 


pag me Vee (Qh Dx) eeepc i=) 42 2B, 
while 
(xg — %o) (¥r— Yo) — (%e — 0) (Yi — Yo) 
= (a5 %) VB) — ie Bn — Yo) 
= (% — 2%) WY = V0) = (% — %) (Y — Yo) 
+ Om l(Y¥ — Yo) + (¥ — ¥o)] + Bn L(® — %o) + (* — %)] 
<= (% — %) (V — Vo) — (x — %o) (Y — Yo) + 204, + 248, 
so that finally 
| mn (8, ¥) — Un (%, Y)| S4(C + MD) ty + 4(D + Ma) B, + 
+ C(x —«)+D(y—y)+2Ma,,B, + M [(% — x9) (¥Y — Yo) — (% — Xo) (vy — Yo) ]- 


Suppose ¢>0 is given. Since 


/\ 


lim «,, = lim B, = 0, 


m—> CO N—> CO 


there are positive integers m, and n, such that 
4(C + Mb) d+ 4(D + Ma)p, + 2MamB,<te 


whenever m>m, and n>n,. Further, in view of the continuity of the functions 
involved, there is a number 6,>0, which does not depend on m and and is 
such that 


C(x — x) + DY — y) + M[(® — x9) (Y — Yo) — (% — %) (¥ — vo) ] <2 


whenever |x —x|<6, and |y—y|<0,. Thus, whenever m>m, and n>m, and 


[Ae ap ces eee Vier eee 
one has Segre 
| nn (%, ¥) -. Umn(%, y)| <6. 
Now, let {tn,n,(%, y)}, where lim m,= lim n,= o, be a singly infinite sub- 
%¥—> CO %1—> CO 


sequence extracted from the double sequence {u,,,,(x, y)}. Given e>0, one 


certainly has m,>m, and n,>n, for all but a finite number of positive integers 7, 
and hence 


| emrernn (%> 51) = Una ay y)| <¢ whenever both | x —x| <6, and ly =| 220.8 


Since only a finite number of values of y are excluded and the corresponding 
finite number of excluded functions Um,n, are continuous (hence uniformly 
continuous) on the rectangle R, it easily follows that the set of functions wu 
where 7=1, 2, 3,..., is an equicontinuous set of functions as desired. 


My Ny? 
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It follows then from Ascoxt’s theorem (see Asco1t [7] or TONELLI [11]) that 
there is a subsequence of {Wn,n,(%, y)} which converges uniformly on R to a 
continuous limit function. This information is all that is really needed to com- 
plete the proof of Theorem 1 of Section 2 (where f depends only on (x, y, 2), 
as can be easily seen by particularizing the considerations of the following sections, 
and for this reason the proof will not be carried out in detail here. 


The formula for u,,,,(x, y) given above was derived by carrying out a step- 
by-step process such as would take place in an actual numerical solution. An 
alternative derivation of the formula for w,,,,(x, y) will now be given. This 
second derivation seems to have the advantage of leading more quickly than the 
step-by-step method to a formula of the desired kind for other boundary value 
problems as well as for the present one. 


First, it will be recalled that if the function F(x, y) is continuous for (x, y) 
in k, and the functions G(x) and A(y) are continuously differentiable on 
XypSxSxyta and yo ySy+, respectively (and G(x») = H(yo)), then there 
is one and only one function w(x, y) which is continuous in R&, together with 
Ow/dx, Ow/dy, and w/dx dy (=Cw/dy Ox) and satisfies the boundary value 
problem 


Pw ieee f 2 
Peon Cae) SC) We tore Cay een ses 
DieVyy = Ge) 10F Aj Ves as 


w(%, y) = H(y) for ySyvSyt+). 


The function w(x, y) is given by the formula 
“u Y 
w (x,y) = G(x) + H(y) — w(x, ¥0) + J SJ F(E,n) dé dy, 
where w(x, Yq) =G(%9) =H(y)- 
Consider the subdivisions 
Hq Hy << Mg Ny Xt 4, 
Vo Vr<Vo< 0 SIM <VMo + 8, 


which were employed in the step-by-step process leading to the equation for 
Um, (x, Vy). By use of this subdivision of the rectangle R, the formula for w(x, y) 
may be rewritten as follows: 


xy 


Rk l 
w(x) =G(x) + Hy) —w(ro yo) +X DY J fF (Gn) ab dn + 


*je-1 Vj—1 


y ey 


+. y [ PFE Mazant ¥ i f Fn) d&dn +S f Fn) dé an. 
J=1 ke Vj—-1 u i-1 VI XK VI 

This rewriting of the equation for w(x, y) makes no difference under the 
assumptions made about the functions F(x, y), G(x), and H(y). But it makes a 
difference when the differentiability and continuity requirements concerning 
F(x, y), G(x) and H(y) are relaxed slightly. Specifically, suppose that F(x, y) 
is bounded in absolute value throughout R and continuous at all imterior points 
of each subrectangle R7)”, with possible discontinuities allowed on the boundary 
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of any such subrectangle. Suppose also that G(x) is continuous throughout 
Xp<x<x)+ a and continuously differentiable for each x interior to a subinterval 
(i.e., such that x,<*<%,,, for some k) but that the derivative of G(x) need 
not exist for the subdivision numbers x,. Similarly, suppose also that H(y) is 
continuous throughout yy< y< yv)+ 6 and continuously differentiable for each y 
interior to a subinterval (7.e., such that y,<y<y,,, for some /) but that the 
derivative of H(y) need not exist for the subdivision numbers y,;. The require- 
ment that G(x») = H (yo) is still retained. Under these relaxed assumptions, the 
rewritten formula for w(x, y) shows immediately that w(x, y) is continuous on 
R and satisfies the partial differential equation 


o2 ) 
(*,9) = gray (9) = Flay) 


o2w 
ox Oy 


whenever (x, y) is mterior to a subrectangle R;’/'. Further 
Ww (4,Vopa= Ge)" fer a, Ss WSS Aut 4, 
w(%,¥) =H(y) for ySySyo+ b. 


This last observation and the rewritten formula for w(x, y) furnish im- 
mediately the desired formula for u,,,,(x, y) upon taking F(x, y), G(x), and H(y) 
to be certain suitably chosen functions. One need only take for F(x, y) the — 
following (piecewise constant) function defined on x x5 x%)+4, vySySvo+d, 
b 
F(x,y) =fei for %SxS%q_i, and ySySyiy, 
while for x=x,+a and y=y,+) 


F(%y+4,¥) =fn-14 for w=¥< Has, 


Fx Yo b) = Te n-1 for XS X< Xp, 1> 


F(% + 4, V9 + 8) = bynes 
where R= On Ay css iia A Ol Oh eo ee 


while for G(x) and H(y), respectively, one takes the polygonal functions (compare 
the description of the Euler-Cauchy polygon method in the introduction): 


o(x,) + ee 7 o(%%) (x — x,) 
Rader 

G(x) = fot: x, 4 ey Mand - FO) teow 

a(x» + 4) for *¥=%+tad=Xp, 
and 

(yy) op Tia) 20) fy — 9) 

H i Jt+1— Vb 

ies for y;Sy<Syjegs and L=O0)1,.,%—4 


T (Yo + 0) for y=y+b=y,. 


In verifying this remark, it must be remembered that, in the abbreviated notation 
one has for example 


> 


Ugo = U(%q, Vo), T(Xp_) = Mo, T(¥)) = Mo,1- 
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§ 4. The convergence inequality 


In order to complete the proof of the theorems of Section 2, one has to con- 
sider the two double sequences of “partial derivatives” with respect to x and y 
of the double sequence of approximating functions of the last section. The 
quotation marks enclosing the phrase “partial derivatives’ are a reminder that 
these functions must be precisely defined on R, especially along the boundaries 
of the subrectangles of R, where jumps may occur. The exact definition of 
what is meant by “‘partial derivatives’’ will be taken up in Section 5. Since the 
“partial derivatives” in question are not necessarily continuous functions on R, 
in considering their convergence one cannot make use of the theorem of Ascott [7] 
on equibounded, equicontinuous sequences of functions employed in Section 3 
above. Instead, appeal will be made to a theorem of ARzELA [7, pp. 119—125] 
asserting the convergence of certain sequences of not necessarily continuous 
functions to continuous limit functions. The lemma of the present section 
furnishes an inequality concerning finite sums which serves as a basis for the 
application of ARZELA’s theorem in Section 5. The result of the lemma is termed 
here the “‘convergence inequality’’ because of the central role it plays in the 
convergence proof of Section 5. It is remarked that in the theory of the ordinary 
differential equation dy/d x = f(x, y), an entirely similar réle is played by another 
convergence inequality (see, for example, Biss [9, pp. 88—89]). The proof of 
the inequality of the lemma below resembles that given by M. Bretor [J8, 
pp. 31—32] for an inequality occurring in the theory of the ordinary differential 
equation dy/dx = f(x, y). Compare also the inequality employed by H. BECKERT 
[22, p. 13]. 

Lemma. // 


(1) tis a pees integer, fy, fx, ---, fy iS a sequence of t+ 1 non-negative numbers, 
and %,2,,---, 2,18 anon-decreasing sequence of t+-1 real numbers (so that z;—z;-420 
ROL ft A oy 8b), 

(2) the numbers L=0 and e=0 ave such that the inequality 


l 
iS oor Ld h- es) 
i 
is valid for 1=1, 2,...,t; then 
(3) the inequality 
t 
psf T+ Le—aa}{et Lhe —a)} 
11 “ 
OLAS HOF 20 —— WE ony te 
Proof. It will be shown by mathematical induction that 


e+ LD hal 41)S {T+ 26-29) I} {e+ L fol —z)} 


for 1=1,2,...,¢, which implies the desired conclusion of the lemma, since 
ee Di e e= 1 for t= 1,2,...,0 and 2-227, (27— 2) = 0- 

For /=1 the asserted inequality follows from hypothesis (2) and the fact 
that 1+ L(z,— 2%) 21, because 


e+Lfy(a—%) = {1 a Lay a) }{e + Lfyla— 2o)} 
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Now for the inductive step. Suppose that the inequality to be shown holds 
for a positive integer /<¢—1, then it will be shown to hold also for the integer 
141 in place of J. This is readily seen, because then, by the inductive hypothesis 
and hypothesis (2) of the lemma, one has 


pet LD ha—4-)S{ UT U+Le—aal}{e + Lhola—al} 


which, together with the equality 
141 { l | 
E+ LX h- 1(% — 4-1) = \é ne Eda oS 3-1) a dpe Le pig =e ap 
i= i= 
implies that 
141 141 \ 
e+ Ld halj—Aa) S| M+ Le —aaJ]p{et+ Lila —a)}, 
fe 


and the proof is complete. 

It is of some interest, although it is not needed in the considerations that 
follow, to point out that the inequality contained in the lemma just proved is a 
finite difference analogue of an inequality due to T. H. GRonwa zt [10, p. 293], 
in the continuous case. (See also G. SANSONE [2/, vol. I, pp. 30—31].) Making 
suitable changes (in order to conform with the present notation) in the statement 
of GRONWALL’S inequality, as given by SANSONE, one obtains the following result: 

If f(z) is a non-negative continuous function defined on the interval z3Sz7< 
%j+a and there exist numbers e20 and L=O such that 


O<f()Se+SL fat 
for 2%SzSz+a, then \ 
OS f(z\ee iors, Sar =, ee 


In order to compare this last inequality with the inequality of the lemma 
proved here, for each positive integer ¢ consider the following subdivision of the 
interval zySzS2,+ a: 


2o =%o1 S41 S%q1S-:: SY 1 ee tot as 


and suppose that the hypothesis (2) of the lemma holds, with z, and f; being 


replaced, respectively, by z;, and f(z;,). Then the conclusion re the lemma 
proved reads 


a ess Le i—Anad]} fe + L fo (21,1 = 20%) 


i=1 
on 
Il ok mana {e |. Dae, — z)}; 
that is, 
f(z) Se” {e+ Lfy(i—%)}, 


whose relationship in the limit to the inequality of GRonwatt cited above, 


iG) see" 
is clear. 
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§5. The double sequence of functions approximating the partial derivatives 
of a solution 


Consider the double sequence of approximating continuous functions u,,,, 
defined in Section 3. It has been pointed out at the beginning of Section 4 that 
the partial derivative with respect to x of w,,,, exists in the usual sense and is 
finite on R save possibly when x is equal to one of the finite set of numbers 
(recall the abbreviated notation introduced at the end of Section 3) 


Hy X's Seats, 


where jumps may occur. (Of course, it is understood that when x= x, and 
X—=X +a, by the “partial derivative with respect to x’’ of the function u 
are meant the one-sided limits 


mn 


lim “mn (%, ¥) = Un (%o Y) 
x 


X—> Xo — Xo 
%> Xp 
and 
lim Umn (%, y) mi Umn (%o + 4, y) 
E> ata % — (%) +4) i 
K<xH+a 


respectively.) A similar statement applies to the derivative with respect to y 
of the function w,,,,, the possible jumps now occurring when y is equal to one 
of the finite set of numbers 

Ve Se nao 


a corresponding agreement being made about the “‘partial derivatives with 
respect to y’” of the function u,,,, when y=y, and y= )+0. For reasons of 
symmetry, it is clear that one may restrict attention to the x derivative, similar 
considerations being applicable in the case of the y derivative. Intuitively 
speaking, it will now be shown, using the lemma of Section 4, that the magnitude 
of the jumps in 6u,,,,/0x can be made arbitrarily small by choosing both m and n 
sufficiently large. 

First, consider the function w,,,, on the closed subrectangle R’”, where 
k=0,1,...,.m—1andl/=0,1,...,2—1. By its very construction, the function 
Um, 1S bilinear in x and y on the subrectangle Rj". In view of the formula for 
Um, given in Section 3, when (x, y) is a point of the rectangle Rj" which is not 
on its closed left and right-hand rectilinear boundary intervals (7.e., when the 
point (x, y) satisfies the inequalities x,<x<%,,, and ySyS Viex), then 


l 
U —uU 
eas @, 9) =e + > faj—a — 3-1) + farlv — »). 
ja 


Ox : Yeti — *k 


On the other hand, when the point (x, y) is on the closed left-hand rectilinear 


boundary interval 
K=%,, MSVSVi41, 


then the right-hand x derivative 


’ 


i: : U X,v)—u Bee 
FUmn Cae y) — lim wince y) all k y) 
Ox : E> XR ¥ — Xp, 
Arch, Rational Mech. Anal., Vol. 4 26 
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where %,<¥<%,4, and y,SySvj4,1, exists and is a linear function of y. Simi- 
larly, when the point (x, y) is on the closed right-hand rectilinear boundary 


interval 
R= hp, WS Sie 


then the left-hand x derivative 


OU an ( Umn(¥, ¥) — Umn (XR +19) 


- = = > 
OX E> Xk 41 Da ip 


where %,<%¥<%,4, and y,Sy<y,,, exists and is a linear function of y. It is 
to be noticed that the “partial derivative with respect to x’ is constant on Rj” 
for each fixed y; that is, for each y such that y,<y<vy,,, one has 

OUmn el OUm n —_ OUmn 

a (%_,¥) = oe (%,y) = Fe (%p+1,9) 
for all x satisfying %,<*< %p,44. 

The maximum absolute value of the difference between the values of the 
“partial derivative with respect to x’ of u,,,, on two subrectangles Rz'”" and Rzy" 
at the same y level will now be estimated by use of the lemma of Section 4. 
Suppose, for definiteness, that k=k. For the rectangle R7/" there are formulas — 
for 0u,,,,/0%, etc., similar to those just derived for Rj", which need not be recorded 
here explicitly. One also has that for each y such that y,S ySy;., the equality 

Ons One Oya 
ox ~ (%5, ¥) = ee (%, 9) =- rr (Xp+1 ¥) 
holds for all x satisfying x,<*<x;,,. Consequently, the problem of estimating 
the maximum absolute value of the difference between the values of the “partial 
derivative with respect to x’ of w,,,, on the two subrectangles Ry” and R*," 
reduces simply to the estimation of maximum absolute value of the difference 
of the two functions of y, 
mn 


O* Um, nw ( , 
Ox 


oO an and XE 

(Xx, ¥) Ag ts Ys 

which are linear functions of y on the interval y,Sy<y,,,. In view of the 
linearity of the two functions involved, the desired maximum absolute value 
of their difference, 

| OU Mn Ot Uy n 


se 
max x aoe eT 
EE ee | Ox (Xx, ¥) ax (5 9)| 


WSIS V4 
is just equal to the maximum of the four numbers 


at 
a) z u mn or u mn 


AK (%p. Yrya) — aa (%R, Yi41)), 
Tet (mv) — Se (a, 9), 
we (%p Vira) — 5 (Xp, yi), 
Tem (ory, 93) = Tm Cg, yeaa). 
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It is to be noticed that only the estimation of the first two of these numbers 
requires special attention, since it will turn out that the last two can be made 
arbitrarily small whenever the difference y,,,— y, is chosen sufficiently small, 
the reason for this being the continuity of é*w,,,,/0x with respect to y for each 
fixed x. For example, 


Ot. otu 
|: oe ~ (Xm Vi4a) — ay ae (xg, 91) 


+ 
a) Umn ot Umn 


Oru otu 
= ee a oa ke (%5, Via) — ge (xe, Yi). 


Aya (Xp, V4) = ee (XR, yee 1) 


Further, since the first number is obtainable from the second merely by replacing 
1 by /-+-1, all that remains is to estimate, for each pair of fixed integers k>k, 
the +141 numbers 


GMa 


Oo Um n 
pms (%,. Yi) — cae (XE, |, 
where /=0,1,...,”. For 1=0 this absolute value can be made arbitrarily 
small, and the lemma of Section 4 will now be used in showing that the absolute 
values for J=1,..., can also be made arbitrarily small. 


or Um nv 


Now, from the definition of ae (x,, ¥;) and 


or =e (Ae, y)) (recall, for 


example, that wu,,,,(%, y), for x, *S%,4,, is a linear function of the single 
variable o) together with the previous formula for 0*w,,,,/0x obtained in this 
section, it follows that (recall that, for example, u,,,,(%,41, Vz) = Up 41,1) 


oUmn WA, ae Lae UR+4, irs UR, u 


a ae (%_,%1) — tee (519) == 


es Ve+1 — *R 


Up+1,0 — Uk, 0 UR, 0 4k, 0, 

— eae (fn,7-1 — fe, j-1) Oy — Yj-1) for J=1,...,% 

‘ 4 i 7 
Vkt+i— ¥k Vht1— +k Ta 


The proof that the absolute value of this last difference can be made arbitrarily 
small provided that m and m are chosen sufficiently large will now be completed, 
at first under the hypothesis required of the function in Theorem 2, 1.e., that f 
satisfies a Lipschitz condition in all three of its last arguments z, p, g (see hypo- 
thesis (1) of Theorem 2). The argument will be carried out first in this case 
because it is somewhat simpler than the corresponding argument when f satisfies 
a Lipschitz condition only in its last two arguments , ¢ (see hypothesis (1) of 
Theorem 3). It will also be supposed at first, again for the sake of simplicity in 
writing, that the function f(x, y, 2, p,q) does not depend on x and y, that is 
f=1@.P.9)- 

Accordingly, under the hypothesis (1) of Theorem 2, one has that (recall the 
description of the abbreviated notation /, ;, introduced in Section 3): 
Up+i,j—-1— YR,j-1 Uk+1,j—1— Vk, j-1 


| 
I 


T 


ees = fe, j-1| aL i a up ;—1| al 


Nht1— “k Mk+1— Xk 


}, 


26* 


Uh Chg | Say Cael 
Vji— Vi-1 Vi — Vj-1 
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and consequently 


| Unti,1— Un,1 UR+1,1 — UR,I ; 
Arti vk Xh+i— *k : 


Un+1,0— %k,0  “E+1,0—— UR, 0 
An+1— “hk AR+1— Vk 


l 
ES eae sal 
7=1 


nie Cia ie Sis 
Vj; — Vji-1 Vj; — Ve—-1 


I 


Upet ee Up ae’ Seah 
Ve+1— %k AR+1— Vk | 


| (oy — ya): 


The term in the last inequality which involves the difference quotients with 
respect to y requires special attention. Consider the function w,,,, on the rectangle 
Ry; -1. From the formula for w,,,, given in Section 3 it follows that 


=o 
Uy, 7-1 = Up,0 + Uo,j-1 — Yoo + bor J- 1( i—1) (vy = V7-1) , 
Ps 
and a 
Un, 7 = Up,o + %0,7 — Yoo 4 ashy: ys fi-1 a= 1(%; — %j-1) (v= Va) 
Fae 
k 
ud 
+ 2 hi-ai-1 (Spe a) IR HAN 
4 
hence j 
Up; — UR, j—-1 Uo,; — Uo, j-1 
= hea (GI — Hale 
Vj Vira Vi Via af 17-1 
Similarly P 
UR; — URj—-1 40,7 — Yo, 7-1 
= — a = aL = E 
y, = Vi Vj Vi— ah Ph s 1) 
and thus 3 
Up 7 Up G=1 UR, j Uf j-1 
Sy ee e. 
Vi Vina Vi Vjna co maja a) 


Further, in view of this 


UR, 7 = Up j—1 Uk, 7 aa UR, j—-1 


Fay te at Wy Aa Vj—1) 
x l nay 
<M| > (%— i-)| oS a =M| »s %— %-2)| -b=Mb(x,—~,). 
t=k+1 j=1 EO 


The inequality for 


may now be rewritten in the form 


| “ae+1,t— “et Vata MR | — |Matio— Uh,o — UR+1,0— YE, 0 
Vk+1— ¥k Hh4+1— 4k | | Hay ¥R “R+1— *E T 
+L > | 4, ;—1 — Uy, ;—-1| (Yj — Yj—a) + Mb (xg — %) + 
j=1 


Un+1, 7-1 — UR jg—1 UR+1, 7-1 — VR 7-1 
Ahi — WR VR+1— VE 


l 
nes! 
j=1 


— Y;—-1) for Vee ted ene nN. 
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This last inequality is precisely of the same type as that of hypothesis (2) of 
the lemma in Section 4, upon identifying, in particular, ¢ with n, the jf; and 2; 
occurring there (for 7=0,1,...,) with the present 


it 3 ii eed 8 Py a 


Veti— %R *ht1— Xk 
and y;, respectively, and the ¢ of the lemma with 


Mn ty0r 7 40 YE+I,07 “Ro! , 
aecb A 


max 
Osk<kSm-l1 


4k+1— Xp 4 +1— VR 


1 
oe LL 2 ppp Up, ;—1| (9, — 92a MO (eR x} ; 
re 


which, as will now be shown, can be made arbitrarily small merely by choosing 
mand n sufficiently large and |x,— x,| sufficiently small (in view of the assumed 
continuity of the derivative o’(x) and the equicontinuity of any subsequence 
{tm,n,(%, y)}, with im m, = lim n,= 0c, which was shown in Section 3). In 


%1—> © 
verifying this, one can use the mean value theorem of the differential calculus, 
since for k=1,...,m—1 


Un+i,o— Uk, o UR+1,0 — “ko 1% - 
SL oat) —a3p), 
Mati “bk Ahi “Fk 
where 2%,< %, <Njiy and x,< 4 <%p14. Let s>0 be given, then: there exist 
(see Section 3) positive integers m, and 2, and a positive number 6, such that 
whenever m,>m,, n,>n, and |x,—x,;|<06, one has 
\ 
max 1 aed War FS ye ea Ys pe ea) ae INO te ee es 

ees, 2 kyj-1 Ri 1| (seta) (x; a) 2€, 
where w(x, y) is written for U,n,(%, y). Also, in view of the unzform continuity 
of the function o’(x) on the interval x *< x +4, it follows that 


“E1,00, “ko | “YRta,0~ “Eo 
Veta Wk *k+1_ “k 


<-— 6. 


2 


Consequently, from the conclusion of the lemma of Section 4 it follows that 


Up+3,1— URI UR +1,1— YR, 1 


Me+1_ %k VR ae OR 


& { Ul [1 +L (y,-9:-)]} {e+£ 6(1—Yo)} S {I oral fet €(¥1—Y0)} 


=e "te + Ley %)} fOreti—t2- sae 


(It should be noticed that the last inequality also holds for /=0.) 


The last inequality has been obtained under the two assumptions that the 
function f satisfies a Lipschitz condition in all three variables z, , g (hypothesis 
(1) of Theorem 2) and that / does not depend explicitly on x and y; that is, 
f=/(z, p,q). The derivation of a similar inequality, in the case when /= 
(x, ¥, 2, B, q) satisfies only a Lipschitz condition in the two variables £, ¢ (hypo- 
thesis (1) of Theorem 3) will now be sketched. As in the previous case, the initial 
step, where the Lipschitz condition is applied, is in estimating the absolute value 
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of the difference f, ;-1—/g,;-1- This can now be done as follows by adding and 


subtracting the number 


Un+1,j--1— “a j-1 , esse 


x Aa ; 
i Rk» Me Sy) eR a D> Xp ) yj;— 


Xe+1— Vk G1 


to the difference in question. One obtains 


i } = flee et jet CR ga eR i=) 
j—1_ /k,j-1 ke Yj—-1s “k,j—-1) ax ; 
k,j—1 X, 7 ) 7 ) 7 Hera — ¥p , yj; — Hin 
/ XEs ay) ‘ Un+1,j-1 UR G—1 FE UR, 7 — Up, fat) ae 
SRS Ds i eee ee 
ae , *kt+1i— *¥k ‘ Vj— Vij 
oe Loe . Mp+ijj—1— Yejj—-1 | UR “RG -1 
+E %Es Vj-as Ua — gies 
Xk+1— Xk bere Sei 


Ui Pasgent Uh gets ee oe 


= |e Means gas 
| : a Akt1— Wk ses | hac Se) 


Using this and the Lipschitz condition with respect to # and g, one has 


Unti1,l — URI UR +1,1— YE,1 


Vet+1— Yk Vk+1— Vk 
Vert Omens CRs ee OO 
sg |e o +¥U -al ,—y--1) 
= 5 zs Rips ae hg 2 et. 
Ve+1— Vk ae | : ‘ | ; 
< Untio— Uk,o0  Uk+1,0— “%k,o| | 
= ic = | 
He+1— *k Xk+1i— *k 
i 
os fx Eas _ Met1,7-1 — UR,j—-1 | UR, j ai | 
k i—1 Ril? oe a : re 
ER pe ee ae Xpt1— Xk a ee 
Upiy j—1— Up mat Up 7 Up peat ; 
= flee; Yeu Mey see tan cee saga hee tat) 0) 
| : 4 tii — Hp” Vg Vad Ak as 


URI, 7 sil < ee jal UR4+1,7 = fir UR, j—-1 


+13] 


Xeta— %p Vh+i— *k 


! (¥; Vj—1)- 


The term in the last summation involving the explicit difference quotients with 
respect to y may be handled exactly as before, yielding the same result: 


1 
Fe gee em 


ae Peg reel 


(¥;—= ¥21) SO (a — 4,2 


j=1 
Thus, one has finally 
Weta — “Rt UR+1,1— UR, 1 
Xet1— ¥p Xk+1— VR 


Urn+1,0 — Uk, 0 Url, ee a 


¥et+1— Xp VR+1— Vk 


Vj D> Up 1) 


Kad Ta OY et SU Al Bees MPa. a pis 
Xp+i— Xp ts poD pee se 


es : . Yneraj—1— Ya j—a , Ukyp — Mega 
lac Dal 3 Ss veer a Viewer pa T) id 


Vk+1— Vp Sale Vi-1 
idan cee TA iam ae bi res Ree ea 
Vkt1— Xp Vk+1— 


OF eal) 


l 
+ MO(%—m) FES Mee 
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This inequality is again precisely of the same type as that of hypothesis (2) 
of the lemma in Section 4, upon identifying, in particular, ¢ with 1, the j; and z, 


occurring there (for 7=0,1,...,) with the present 
Led Rd 7 eee Se ed | 
| *kt1— 4K *h+i— *k 


and y,, respectively, and the « of the lemma with 


u — Up Pate eed) 3 
max k+1,0 OE Ae eeG oR, 0 
0Ssk<kSm—1 


Nei1— Xp Vk+1— *k 
1 
ere ; ee ge le ge ied 

aE PE Gel rae, j= tos Xhiy— Hb 7 yy, 

j=1 \ ri 7 Fiza 

Up ee oe Up We ean \ 

ss a ’ c s é spat Py Nat | bya 2) Pies | pte 

i[ves DO as tha — Xp ray Sean ) (Y¥j—¥j-1) + Mb (xg—%) 

am ; a 


which (as will now be indicated, without entering into the detailed argument) 
can be made arbitrarily small (7.e., less than any positive number given in ad- 
vance) merely by choosing m and » sufficiently large and |x,— x,;| sufficiently 
small. In showing this, use is made of the assumed continuity of the derivative 
o'(x); of the equicontinuity of any subsequence {hsric (5 y)} with him mM, = 
lim 1, = 0c, which was shown in Section 3; and of the wniform continuity of the 


function /(%, y,z,,q) on any closed and bounded set of points (x, y, z, p, 9) 
satisfying 

Ca) i Le pee SP a ee) eae aa), 
with Z, P, Q positive numbers. Notice that it can readily be seen, from the 


definition of w,,,, and of the difference quotients involved, that there exist 
positive numbers Z, P, Q such that 


|v,1| SZ, 
Rees NP, a ior b= 0,4, Bak = 0 and: £0, (eae, 
Xe+1— *k 
nl Nera ama aE 0 ator ON on andes l= 0541, 2... 4 
Viti — Vj 


and for any pair of positive integers m and m, where one uses the abbreviated 
notation, %,;=Umn(X,, vi), etc. In particular, since 


Up+y,1— Up, 1 Unt, 0m “4k, 0 


Npiy— Xp NRiy— Vk 


Pie Ci b 


1 
+ fea = Vi24)> 
j=1 


one may choose 


in terms of the constants C, M and b of Sections 2 and 3. This being granted, 
one obtains exactly as before, by an application of the lemma of Section 4, that 
if e>0 is given, then there exist positive integers m, and m, and a positive 
number 6, such that whenever m,>m, and n,>n,, and | %,— x,|<6, then 


Until — “k,1 UR+1,1 — YVR, 1 


ae se Le(y,—= yor) for L=0,1,2,...,%. 


Vpiit— Wk Mk+i— *k 


382 \PBSDrAz: 


It is now time to define the double sequence of functions {Pn (x, y)} cor- 
responding to the double sequence oer 9: (x, y)} of Section 3. The double sequence 
{p,,,(x, ¥)} will be, roughly speaking, a sequence of functions approximating 
the partial derivative with respect to x of a solution. In view of the possibility 
of jumps in 6u,,,,/0x, the function #,,,,(x, y) has to be defined carefully in R, 
to make sure it is single-valued. For each pair of positive integers m and 1, 
the function #,,,, is defined as follows, for points (%, y) in the closed rectangle R: 


enn (x,y) whenever %,<%*<%414, for some k=0,1,...,m—4, 
Pinn (X,Y) = is oe (%,y) whenever x=x, for some k=0,1,...,m—1, 

z sae (x,y) whenever x =%,,=% +0. 
The function #,,,, possibly has jump discontinuities only when x= %},..., %,—-4 


and is continuous in the two independent variables x and y at all other points 
of R. 

This double sequence of functions {Pmn he y)}, as may be readily seen from 
the formulas given for 0u,,,,/0x, O*u,,,/0x and d-u,,,/0x given earlier in this 
section, is equibounded in absolute value on Rk. That is to say, there is a positive 
number P, which is independent of m, and of (x, y), such that 


| Peel y)| = Pe. 


for any positive integers m and n, and any point (x, y) of R. 


Let {2,,, »,} denote any singly infinite subsequence of functions (with 
lim m,= lim n, =o) extracted from the double sequence {Pan (3 yr. Let s=0; 


s—> oo S—> co 
From the preceding considerations it follows that there exist positive integers m, 
and n, and a number 0,>0 such that whenever (%,V) and (x, y) are points of 
R satisfying 

lA e050, 


and m,>m,, n,>n,, then 
| Dieesetshs) AST Pig al y)| < é. 


(In ARzELA’s terminology [7, p. 119], the subsequence of {Pac ie, y)} for which 
m,>m, and n,>m, is equioscillating by less than ¢. This can be proved by an 
argument similar to that used in Section 3 in showing that the sequence {tna ee, y)} 
is equicontinuous. There are again four cases to consider, depending on the relative 
positions of the points (¥, vy) and (x, y) with respect to each other. As in Section 3, 
only the case when x<¥x% and yS<jyj need be considered in detail. Here one has 
ga N19 MISS Migas and x, SXS p44; YySVSyp4, for suitable pairs of 
integers (k,l) and (k,/). Further x,<x,; y,<y, and Kye Ss Mpa Vee a Va 
The inequalities deduced earlier in this section for 


UR+1,1— URI Ghar t, be eRe 
Vet+1— *k vh+1— *h 
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may then readily be employed to obtain the desired result, the details being 
as follows. Now 


Pn, y) ive Phene y) aa [Pn (%¥) a. Pug (x, y) | a5 Pag (x, y) rin Dine y)I, 


where the point (x,y) is in the subrectangle Rz;" because xg XX xg, and 
¥iSyvSyi41. This, together with the definition of the function 4#,,,,, implies the 
inequality 


| Pn (X,¥) = Pina he y)| =| Pi eV) — Pant x »¥)| + Doma = Paatts y)| 


{eee | (7 bi Mi 


11S MRT appa ron 
BGS — Oe PNY yr) f a | | ¥hi4 — XR =i (y 1) far} 
UR + — Us nee eye 
eS + (y — 91) fra} — | + Cy — 9) fa} 


Nkt1— Wk Xe+1— Xk 


| 
rat 


| erg ORY a tig ee UR+1,1 — URI Sid Jae pind A 

— | ¥e41— HR ¥h4+1— ¥F XR ia— *R thty— Xp | 
1h) — ~ te z — ie 
zt Ms v1) tet| + 2|(y — i) fer] + 1(y — 0) fail 

= UR+1,1— URI Un+1,1— Up, 1 

= 5 ie 1( 5 Vi eis = »| + RS — oe x = 4M B,,, 
j= i+-41 k+1 di 3 AR+1 k 

in case %,5%*%S%,., for some integer R=0,1,...,m—1 (this only excludes 


x= 4% +a, which will be treated separately below) and the integer & is chosen 
(if possible) so that x, <*%<4%,,,, with R=0,1,...,m—1 (if ¥=x)+4, which 
is seemingly excluded at first, the inequality just ten still continues to hold, 
but with k replaced by m—1). If x =x +a, a case definitely excluded above, 
then one must necessarily have ¥ =x +a (=x,,) too, and then 


Panne — Prat (x, y)| 
Um I— Um—-1,0 = i i if Um, i 'm—1, ee 
ea Tea id 6 et Ee Wp 


ae Xm—-1 Xm — i rf 


ze 


+2M6,. 


fn—ri-al = wie) 


These leet Hi now readily furnish the desired “‘equioscillation”’ property of 
the singly infinite sequence {P,,,,,}- 

Since the sequence of functions {Pm, xf is equibounded in absolute value, and 
since for each «>0 there are positive integers m, and », and a number 6,>0 
such that for all points (¥,y) and (x, y) of R satisfying |x —x|<06,, |y—y|<6, 
and for all m, and n, satisfying m,>m,, n,>n, one has 


| Pngn, (%) V) a Pina (%; y)| <r 


it follows from a theorem of ARZELA [7, pp. 119—125] that there is a continuous 
function p(x, y) defined on R and a subsequence of the sequence {/,,,,,} which 
converges uniformly to the continuous function £(x, y) on R. For a proof of 
this particular result needed here, carried out under the equivalent hypothesis 
that the given sequence of functions has zero “‘Grenzschwankung”’ (see CARA- 
THEOpORY [17, p.3] for the definition of this term), reference is made to 
H. BECKERT [22, pp. 24—27]. 
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For reasons of symmetry, without further discussion it is clear how the 
double sequence {9,,,,(«, y)}, which approximates the y derivative of a solution 
is defined. It is also clear that there is a positive number Q, which is independent 
of m,n and of (x, y) such that 


l9mn(%, ¥)| SO 


for any positive integers m and n and any point (x,y) of R. Let {4m,n,¢ denote 
any singly infinite subsequence of functions with iim my = lim N,=oo extracted 
from the double sequence of functions {On ae y)}. Again, by ARZELA’s theorem, 
one concludes that there is a continuous function q(x, y) defined on R and a 
subsequence of the sequence Cie (4%, y)} which converges uniformly to q(x, ¥) 
on R. 


§ 6. The existence of a solution 


Consider the double sequences of functions {tn n (x, y)}, {Pian (%: y)}, and 
{mn(x, y)}. In Section 5 it was pointed out that there exist positive numbers 
Z, P, and Q such that for any positive integers m, m and any (x, y) in R, one has 


lee alte Vy le Se 5 © i ')) je ma nn %, Y= @. 


It is remarked, since use will be made of this fact immediately, that the continuous 
function f(x, y, 2, p, g) is uniformly continuous in (4%, y, z, #, g) on the closed and 
bounded five dimensional set of points defined by 


MySKS%+a, MWSyvySt+4, |z| SZ, lel=P; lq] <Q. 
That is, given ¢>0 there is 6,>0 (which may be chosen to be less than ¢, for 


later convenience) such that whenever (%,, V1, 2, 21,91) and (%», Ve, 2s, Be, Je) 
satisfy the inequalities 


kos hiss tg A; Voss V; Vo 4, \24\eeZ. \P VSP, 1G) =O. fierce =e 


and 


(ee Males, [V1 — Val <6,, |2, = 39| <6, |: — pal <6, 19. — 9s| < 6,, 
then 
|#(%1, 71,71, Pi. U1) —F (%e, Vo, Ze, Bo, Ge)| <e. 


Let {Uny Paden y)}, where lim m, = lim n, =o, be a singly infinite subsequence 
v—> CO v—>O 


of the double sequence {thin n (%, y)}, and suppose further that (see Section 3) 


Dm tm, n, (%, ¥) = 4 (%, 9), 

where the convergence to the continuous function w(x, y) holds over the rectangle 
R. From Section 5, it follows that the corresponding subsequence {int nhs; y)} 
itself contains a subsequence which converges uniformly on R to a continuous 
function p(x, y). For simplicity, suppose the subscripts have been chosen so 
that the subsequence {,,, »,} itself converges uniformly on R to p(x, y). Making 
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a similar agreement about subscripts, it may also be supposed that the cor- 
responding subsequence {4n,n,(X, y)} itself converges uniformly on R to a con- 
tinuous function g(x, y). Summarizing, one concludes that 


1M yng n, (%s ¥) = U(x, 9), 
EM. Pm, n, (x, y) a p (x, y) 2 


BM Ym, n, (x, y) = q(x, y) , 


the convergence to the continuous functions u, p, g being uniform on R. It will 
now be shown that the function w(x, y) is a solution of the boundary value 
problem under study. 


In view of the above mentioned uniform continuity of f on a certain closed 
and bounded five-dimensional set of points, it follows that 


Jim 7 (x, y; Um, n,(%, y), Pinan y), Cpa es v)) = f(x, y, u(x, y), p(x, y), q(x, y)), 


the convergence to the continuous limit function being again uniform on R. 
Furthermore, since the limit function 


i(x,y, u(x, 9), P(%,¥), 9(%,¥)) 
is continuous on &, the following Riemann integrals exist for all (x, y) in R: 
J St(E.n, (En), b(n), (En) dé an, 


Si(E.y. wey), BEY), 9.9) a, 


fi (4m), B (on). 9m) dn: 


the order of integration with respect to € and 7 may be interchanged in the 
double integral without altering its value. All this information will now be used 
in order to show that the function # is precisely the x derivative of the function 
uw and that the function q is precisely the y derivative of the function w. 


Let ¢>0, and let 6,>0 be such that s>6,>0 (the restriction ¢> 6, is made 
for later convenience) and that also 


[f(®1, Va, 2, Pr, M1) — f(¥2» Ves Z2, Pa, Va)| <é 
whenever the points (%;, V;,2;, £;,7;) satisfy both 
py aa, NS. |AbSP,. igi) SQ, 
and 
|x, —%|<6,, |¥1—¥el <4, We 2a) =<. 0s,° | Py— Pal = 0.) (@y — Fa] = O.. 


In view of the uniform continuity of the functions u(x, y), f(x, y) and q(x, ¥) 
on R, there is another number 6*>0 (which for convenience will be chosen such 


386 J. B. Diaz: 


that ¢>6,>6*>0) such that 
|“ (Em) — # (Fo, %2)| <O<e, 
| (E1.m) — P G2, 12)| << O<e, 
|7(&,m) — 9(E2,%2)| <O.<e, 
whenever the points (€,, 1%), (€, 42) of R satisfy the inequalities 
\é,—&|< 6%,  |m— | < oF. 
Further, there is a positive integer N, such that 
|~(E,1) — Um,n,(E,)| < OF < 6,, 
[PAE i) Pree ees Oe 
[PE — an ae en a Oe Oey 
and (cf. Section 3 for the definitions of «,, and B,,) also 


Cn, és <0, fs Pnccoate yee ah 
whenever 


mn, > IN) ot SING, 


and (&, 7) is any point of the rectangle R. 


Let (x, y) be a point of R and m, and n, be positive integers such that mm, > N, 
and n,>N,. These positive integers m, and , and the numbers e, 6,, 6* will 
be supposed fixed during the immediate discussion. The notation of Section 2 
(for example, writing x, instead of x,,,,) will be used in the next computation 
for simplicity in writing. There are integers k and 1, with o<k<m,—1 and 
OSlSn,—1, such that 4<*%<x,,, and y,Sy<yj,.4, i.e., such that the point 
(x, y) being examined lies in the closed subrectangle R?y”". Seca that o (0) = (0) 


and consider the difference 


Un, PCS y) 


*o Vo 


which may be written 


time m%, Y) — [0 (2) + (9) = 6 (0) + 


+E LF HEn wn), Pen). alm) aban + 
+3 a Leona u(E.m), bE. m), a(E.m) dE dy + 
i ; f #(&0,«(E,n), 6 Em) (Em) a8 dn + 


t=1 xy yy 


eA GL es w(E,n), b(&,n), ¢(E, n)) 4& dn]. 


Xk VI 


a(x) +r(y) SI CURIGUWIGU NIGH ELE 
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Recalling the definition of u,,,,,(*, y) from Section 2, and the fact that 


fi-1,;-1= F(a Vis ore ae Mas Vi—1)> Pang tog %i—1 Vi—1) > Trig tp 41 » Vj—1) 


from Sections 2 and 3, one has, for example, that 


lfi-1j-1 — F(,n, u(E, 0), (En), 9(E,n))| 
= |f (4-4; Vi—1 > Voie n: %—15 Vj—1) 5 Pm ne X19 Via) > Yura 13 Vi) = 
PY (Mp2 Vas WM Vj—1)> P (X41, Vj)» V(%—-1 ¥j-1)) | + 
ae Coe Vj—1» W(%; 1, Ya), P(% 1, Vj) 7 (4-1, Vie 1) i 
— (én, u(€,n), b(é,n), gE ))| <2e, 


whenever x,_,Sé<¥%, and y;_;<7<y,;. Consequently, the absolute value of the 
difference u,,4,—[...| is less an (see Section 2 for the definition of the con- 
stants C and D) 


|o(x) —o(%)| +|t(y) — ti) + C-(*— x) +D-(y—y) + 


+ 28 (%, — Xo) (i — Yo) + 2e(% xs) ( 2 (9) — yj-a)) + 


j=1 


+ 2e(y—y)( 3B (%—-2)) + 20 — m9) 9 — 9) 


whenever m,>N, and n,>WN,, and hence 


u(x, y) = a(x) + t(y) - HI ING y, W(E,), B (En), 9(E,n)) aE an 


for any (x, y) in R. From this last equality it follows that du/éx, du/dy, and 
ulox dy (=u/dy Ox) exist and are continuous throughout the rectangle R. 
As a matter of fact 


2H (x,y) =0'(x) + f Hem u(xsn), (en), gem) dn, 


FeO =O) + JIE Hula), BEY). aE 9)) a8, 


while Cu 


— (%, y) = 


255 (x, 9) = f(x,y, u(x, 9), B(%y), a(x), 


55 i 


for any (x,y) of R. 

The proofs of Theorems 2 and 3 will be complete once it is shown that 
ujox=p and du/dy=gq. It suffices to consider only du/dx. Let e>0 be given, 
and the numbers ¢>6,>6,>0 and m,>N,, n,>WN, be as in the argument 
just carried out. Let (x, y) be a point of Rk. There are two cases to consider: 
either x,S*<%,4, and y,S yy, for suitable integers & and /, with OSkS 
m,—2 and 0S/Sn,—1 OF %m, 15% S%m,=%ta, and ySySy,4, with 
0</l<n,—1. Consider the difference 


Prvnd 9) —[o'o) +S Hon e(en), BC) ge 0) ar), 
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which may be written (in either of the two cases mentioned, with k=m,—1 in 
the second case) 


ants») —[o'e) +O F fs ees), boa), a(e1)) dn 


$F (xn wlssn)s Pen), gem) dn) 


Recalling the definition of Py,»,(x, y) from Section 5, and the fact that from 
Sections 2 and 3 


fiaj-1= Coma V7 Urry nA Xi—1 Vi=1)) Din ph Ay Y;-4)) Om bian y;-1)) , 


one has again, for example, that 


Pi een f(é,, w(&,n), p(é,m), q(é,n))| =< 2ZE- 


whenever %,_;S&<x, and y,_;<7Sy,. Besides, the mean value theorem of 
the differential calculus and the definition of the constant C of Section 2 imply 
that 

i r1,0— Uh, 0 


o' (x*) — o'(x)| S C|x* —x| = Ca,. 
*k+1— ¥k 


— 0'(x) 


Consequently, the absolute value of the difference #,,,,,,(%, y) —[...] is less than 


> 


R 


eC + 2e( D (%—%-1)) + 2e(x —) < e(C + 2a), 


4=1 


whenever m,>WN, and n,>N,, and hence 


p(x, ¥) =o'(x) + J f(x, u(x, 9), P(x, ), a(x, 9)) dy 


for any (x, y) in R. Since the right hand side of the last equation is already 
U ou 


known to be equal to ae (x, vy), 1t follows that sao =p, as desired. By sym- 


metry one has also that 
g(x y= ry) + SHE y, w(E), 09), 96,9) dé 


for any (x, y) in R, from which it follows that 6u/@y=gq, and the proof is complete. 


Under the hypotheses of Theorem 3, the preceding argument shows that any 
singly infinite subsequence {%p,»,(X, y)}, where lim m,=co and limn,=oo, 


1T> ow 


contains a subsequence which converges uniformly on R to a solution. On the 
other hand, under the hypotheses of Theorem 2 (in which case there is but one 


solution) the preceding argument implies that the whole double sequence {thin n(%s y)} 
converges to the solution, 7.e. that 


lim u#,,(%, ¥) 


m—> 
n—> 0O 


is the solution, the convergence being uniform on R. 


The Polygon Method 389 


References 


Asco_t, G.: Le curve limite di una varieta data di curve. Memorie della Reale 
Accademia dei Lincei, Ser. 3 18, 521—586. 

PEANO, G.: Sur le théoréme générale relatif a l’existence des intégrales des 
équations différentielles ordinaires. Nouvelles Annales de Mathématiques 11, 
79—82 (1892). 

Brancut, L.: Applicazioni geometriche del metodo delle approssimazioni succes- 
sive di Picard. Rendiconti della Reale Accademia dei Lincei, Ser. 5 3, 143 —150 
(1894). 

ARZELA, C.: Sull’esistenza degli integrali nelle equazioni differenziali ordinarie. 
Memorie della Reale Accademia delle Scienze dell’Istituto di Bologna, Ser. 5 
6, 131—140 (1896). 

Pricarp, E.: Sur les méthodes d’approximations successives dans la théorie des 
équations différentielles. (Note I to vol. 4 of G. Darspoux, Lecons sur la 
Théorie Générale des Surfaces, Paris, 1896, pp. 353— 367.) 

FusIn1, G.: Alcuni nuovi problemi, che si presentono nella teoria delle equazioni 
alle derivate parziali. Atti della Reale Accademia di Torino 40, 616—631 
(1905). 


[7] ArzeELA, C.: Esistenza degli integrali nelle equazioni a derivate parziali. Memorie 


della Accademia delle Scienze dell’Istituto di Bologna, Ser. 6 3, 117—141 
(1906). 


[8] Monret, P.: Sur les suites infinies de fonctions. Annales Scientifiques de |’Ecole 


[16] 


Normale Supérieure 24, 264—268 (1907). 
Briss, G. A.: Fundamental existence theorems. American Mathematical Society 
Colloquium Publications 3 (1913). 


| Gronwat1, T.H.: Note on the derivatives with respect to a parameter of the 


solutions of a system of differential equations. Annals of Mathematics, Ser. 2 
20, 292— 296 (1919). 


| Tone ttt, L.: Fondamenti di Calcolo delle Variazioni, vol. 1. Bologna 1922. 
| Ince, E. L.: Ordinary Differential Equations. London 1927. 
|] Tone ct, L.: Sulle equazioni integrali di Volterra, Memorie della Reale Acca- 


demia delle Scienze dell’Istituto di Bologna, Cl. Sci. Fis., Ser. 8 5, 59—64 
(1927/28). 


] Lewy, H.: Uber das Anfangswertproblem einer hyperbolischen nichtlinearen 


partiellen Differentialgleichung zweiter Ordnung mit zwei unabhangigen Ver- 
anderlichen. Math. Annalen 98, 179—190 (1928). 


] Tonettt, L.: Sulle equazioni funzionali del tipo di Volterra. Bulletin of the 


Calcutta Mathematical Society 20, 31—48 (1928). 

KAMKE, E.: Differentialgleichungen reeller Funktionen. Leipzig: Akademische 
Verlagsgesellschaft 1930. 

CaRATHEODORY, C.: Variationsrechnung und partielle Differentialgleichungen 
erster Ordnung. Berlin: B. G. Teubner 1935. 


] Bretot, M.: Les Principes Mathématiques de la Mécanique Classique. Paris: 


B. Arthaud 1945. 


] Gruxrano, L.: Generalizzazione di un lemma de GRoNWALL e di una diseguagli- 


anza di PEANO. Rendiconti dell’Accademia Nazionale dei Lincei, Cl. Sci. 
Fis. Mat. Nat., Ser. 8 1, 1264—1271 (1946). 


] Farpo, S.: Su un teorema di esistenza di calcolo delle variazioni e una pro- 


posizione generale di calcolo funzionale. Annali della Scuola Normale Supe- 
riore di Pisa, Ser. 2 12, 119—133 (1947). 


] Sansone, G.: Equazioni Differenziali nel Campo Reale, Seconda edizione, Parte 


prima, 1948, Parte seconda, 1949; Bologna: N. Zanichelli. 


2] Brecxert, H.: Existenz- und Eindeutigkeitsbeweise fiir das Differenzenverfahren 


zur Losung des Anfangswertproblems, des gemischten Anfangs-Randwert-, 
und des charakteristischen Problems einer hyperbolischen Differentialgleichung 
zweiter Ordnung mit zwei unabhangigen Variablen. Ber. Verh. Sachs. Akad. 
Wiss. Leipzig, Math.-naturwiss. Kl. 97, H. 4 (1950). 


- Aw ‘ eel ‘ pe a ee 
‘ 4 ee Meas eee at ‘* A, Mi, Na esis ae Ba 
1 vu Ty oie nae oes Ace y 
v " "¥ q mie 
p - $i oes - i ae 
a", cA, 


* r ; ‘, ie 
J. B. Diaz: The Polygon Method 


LrEHEY, P.: On the existence of not necessarily unique solutions of classical 
hyperbolic boundary value problems for non-linear second order partial diffe- 
rential equations in two independent variables. Ph. D. thesis, Brown Uni-_ 
versity, June 1950. . 

ZWIRNER, G.: Sull’approssimazione degli integrali del sistema differenziale 
02 2/ax Cy =f (%, ¥,2), 2(%o, Vv) =Y(V), 2(¥, Vo) = v(%). Atti dell’Istituto Veneto ¥ 
di Scienze, Lettere ed Arti, Cl. Sci. Fis. Mat. Nat. 109, 219—231 (1950/51). 

Contt, R.: Sul problema iniziale per i sistemi di equazioni alle derivate parziali 
della forma 2 = f (%, y;2™,..., 2; 2), Rendiconti dell’Accademia Nazio- 
nale dei Lincei, Cl. Sci. Fis. Mat. Nat., Ser. 8 12, NotaI, pp.61—65, Nota II, 
pp. 151—155 (1952). 

HARTMAN, P., & A. WINTNER: On hyperbolic partial differential equations. 
Amer. Journal of Mathematics 74, 834—864 (1952). 

Conti, R.: Sul problema di Darboux ‘per l’equazione 2,,=f/(%, V, 2, ¥x, 2y). 
Annali dell’Universita di Ferrara, Nouvo ser., Sez. 7, Sc. Mat. 2, 129—140 
(1953). 

Coppineton, E. A., & N. Levinson: Theory of Ordinary Differential Equa- 
tions. New York: McGraw-Hill 1955. 

Moors, R. H.: Proof of an existence and uniqueness theorem of Picarp for 
a non-linear hyperbolic partial differential equation, M. A. thesis, University 
of Maryland, June 1955. 


ALEXIEWIcZ, A., & W. OrLIcz: Some remarks on the existence and uniqueness ~ 


of solutions of the hyperbolic equation 6?z/dx dy =f (x, y, 2, @z/Ox, 0z/dy). 


Studia Mathematica 15, 201—215 (1956). ; 


Institute for Fluid Dynamics and Applied Mathematics 
University of Maryland 
College Park, Maryland 


(Received February 4, 1958) 


Applied Mathemaetess ’ 
and Statistics Laboratory 
Stanford University, California 


R. BERKER 


Technical University 
Istanbul 


L. CESARI 


Purdue University 
Lafayette, Indiana 


PuCOL TAG 


Institut fiir Angewandte Mathematik 
Universitat Hamburg 


J. L. ERICKSEN 


The Johns Hopkins University 
Baltimore, Maryland 


A. ERDELYI 


California Institute of Technology 
Pasadena, California 


IRMGARD FLUGGE-LOTZ 


Stanford University 
California 


HILDA GEIRINGER 


Harvard University 
Cambridge, Massachusetts 


H. GORTLER 


Institut fiir Angewandte Mathematik 
Universitat Freiburg i. Br. 


D. GRAFFI 


Istituto Matematico ,, Salvatore Pincherle“‘ 


Universita di Bologna 


Av. GREEN 
King’s College 
Newcastle-upon-Tyne 


J. HADAMARD 


Institut de France 
Paris 


H. LEWY 


University of California 
Berkeley, California 


A. LICHNEROWICZ 


Collége de France 
Paris 


CaCeEEN 
Mathematics Research Center-U. S. Army 
University of Wisconsin 
Madison, Wisconsin 


~ 


“ED ITO RIAL BOARD 
ce -# = i a 
S. BERGMAN athe. 


C. LOEWNER 


Applied Mathematics 
and Statistics Laboratory 
Stanford University, California 


W. MAGNUS 


Institute of Mathematical Sciences 
New York University 
New York City 


G. C. McVITTIE 


University of Illinois Observatory 
Urbana, Illinois 


J. MEIXNER 


Institut fiir Theoretische Physik 
Technische Hochschule Aachen | 


C. MIRANDA 


Istituto di Matematica 
Universita di Napoli 


C. MULLER 


Mathematisches Institut 
Technische Hochschule Aachen 


A. OSTROWSKI 


Mathematische Anstalt der Universitat 
Basel 


R. 3S) RIVEIN 


Division of Applied Mathematics 
Brown University 
Providence, Rhode Island 


M. M. SCHIFFER 


Mathematics Department 
of the University 
Amsterdam 


E. STERNBERG 


Division of Applied Mathematics 
Brown University 
Providence, Rliode Island 


W. TOLLMIEN 


Max-Planck-Institut 
fiir Str6mungsforschung 
Gottingen 


A. TONOLO 


Seminario Matematico 
Universita di Padova 


C. TRUESDELL 


Mathematics Research Center-U. 5S. Army 
University of Wisconsin 
Madison, Wisconsin 


He VIL eAgT 


47, bd. A. Blanqui 
Paris XIII 


CONTENTS 


EricksENn, J. L., & C. TRUESDELL, Exact Theory of Stress and Strain 


in Rods’and ‘Shells. .2 4. .- Sp pure los, eee ZO 
Copson, Ey T., On the Riemann-Green, Function 2-75.40 eee: 


Corson, E. T., On a Singular Boundary Value Problem for an Equation 


of Hyperbolic Ky pers... ks —. ac. ee ok ee 


Diaz, J. B., On an Analogue of the Euler-Cauchy Polygon Method for 


the Numerical Solution of uzy=/(%, y, u, Uig Uy) > de G~ ae oe ee Se 


Druck der Universitiatsdruckerei H. Stiirtz AG., Wiirzburg 


